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ABSTRACT 


Although numerous previous measurements of the radioactivity of terrestrial 
materials have been made, most of these researches are of qualitative value only in 
their applications to geology, geophysics, and cosmology because of inadequate recog- 
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nition of the analytical care necessary in order to obtain reliable results. In the 
present study of the radioactivity of terrestrial materials, a systematic program of 
standardization, calibration, and interchecking has been followed throughout. As 
part of an international intercalibration among the various workers in this field and 
a general program of helium age research, several hundred radioactivity measure- 
ments have been made. These results represent the most reliable collection of 
radioactivity determinations which have yet been made within the range of concen- 
trations involved. By combining these newer measurements with the limited number 
of well-authenticated earlier analyses available, average values have been obtained 
as follows: 1.37 + 0.17 X 10 g Ra/g for 43 acidic igneous rocks, 0.51 + 0.05 x 10 
g Ra/g for 7 intermediate igneous rocks, 0.38 + 0.03 X 10™ g Ra/g for 54 basic 
igneous rocks, and 0.57 + 0.08 & 10” g Ra/g for 28 sedimentary rocks; 3.0 + 0.3 x 10° 
g U/g, 13 + 2.0 X 10° g Th/g, and a Th/U ratio of 5.0 for 26 acidic igneous rocks; 
14 +02 10° g U/g, 44412 10° g Th/g, and a Th/U ratio of 26 for 6 inter- 
mediate igneous rocks; 0.96 + 0.11 « 10° g U/g, 39 + 06 X 10° g Th/g, and a 
Th/U ratio of 4.0 for 34 basic igneous rocks. These values are substantially lower 
than those obtained by Jeffreys in a compilation of most of the measurements re- 
ported prior to 1936. The present averages show a more marked decrease of radio- 
activity with increasing basicity, the Th/U ratios are considerably greater than 
those compiled by Jeffreys and are in better agreement with those to be expected 
from geochemical considerations. Two ultrabasic rocks were found to have radio- 
activities comparable to the low values for iron meteorites. Specific inaccuracies 
in earlier investigations have been discovered. Estimates are made of the rate of 
production of heat by radioactive decay based on the above average values for 
the different rock types. 


INTRODUCTION 


Knowledge of the radioactive content of terrestrial materials is of 
fundamental importance in many phases of geology, geophysics, and 
cosmology. The problem of the internal heat of the earth, and its relation 
to the source of magmatic energy, involves accurate information con- 
cerning the rate at which radioactive heat is being generated in the vari- 
ous types of rocks. The variation of radioactivity with rock type is also 
of interest in the study of the differentiation of magmas. The distri- 
bution of the radioactive elements among different minerals within a 
given rock may have important implications in geochemistry. Geological 
age measurements by radioactive methods depend upon accurate determi- 
nation of the content of the long-lived parent radioactive elements in 
rocks and minerals. The radioactivity of sea water and marine sedi- 
ments is illuminating in certain aspects of oceanography. The in- 
creasing applications of radioactive methods in geophysical prospecting 
focus attention on the need for further information on the radioactivity 
of rocks, particularly those of sedimentary origin. Measurements of 
the radioactivity of crude oils and associated sediments suggest that the 
chemical effects of alpha radiation may be important in the genesis 
of petroleum. 

While it is true that the determination of relative magnitudes of radio- 
activity has some value, the truly significant interpretations in geology, 
and in other fields of science, will be possible only after the radioactivi- 
ties of terrestrial materials have been established on an absolute basis 
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and have been evaluated in terms of the accurately known physical con- 
stants of the radioactive elements. Toward this goal, several hundred 
radioactivity measurements have been made on a wide variety of 
terrestrial specimens. Precautions have been taken to eliminate as 
many of the usual uncertainties as possible. A systematic program of 
standardization, calibration, and interchecking has been followed through- 
out this research. As the result of this care, it can be said with con- 
siderable confidence that these data represent the most reliable collection 
of radioactivity measurements which have yet been made within the 
range of concentrations involved. 
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RADIOACTIVITY MEASUREMENTS OF ROCKS 
GENERAL CONSIDERATIONS 


Most of the radioactivity measurements of rocks reported in this 
paper were made as part of the general program of helium age research. 
Other determinations resulted from investigations on the genesis of 
petroleum, the standardization of radioactivity measurements, and the 
international intercalibration program. Practicality dictated that each 
measurement have many-fold applications if possible. Consequently, 
the rock specimens show a somewhat less systematic geologic distribution 
than would have been selected had these measurements been made solely 
for the purpose of establishing the averages indicated for the different 
lithologic types. 

RADIUM MEASUREMENTS 

Techniques——The direct-fusion method of releasing radon from rocks 
as developed by Evans (1935, p. 102) has been employed, with essentially 
no modification, in all the radium measurements. Two techniques were 
used in determining the ionization effects of the radon and its immediate 
decay products. Total ionization measurements were made with the 
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double ionization-chamber, string-electrometer, photographic recording 
instrument (Evans, 1935, p. 104). Other determinations employed a 
more recently developed apparatus that counts the individual alpha rays, 
emitted in a single cylindrical ionization-chamber, by the use of the 
usual type of vacuum-tube electrometer (Penick, 1935, p. 115). Photo- 


200 T if T 
° 
° 
100 3 
a 
E 4 
° 
+ 
|__CURIES OF RADON x 107 : 


Ke) 2.0 3.0 
Ficure 1—Calibration of alpha-counting radon apparatus 


graphic recording affords continuous observation up to 20 hours on a 
single sheet of bromide paper 6 by 37 inches. This apparatus is a 
liberally modified version of the one described by Poole (1938, p. 595). 


Calibration.—The calibration of the double ionization chamber appa- 
ratus has been discussed in detail elsewhere (Evans et al., 1939, p. 937). 
The alpha-counting radium apparatus has also been calibrated with 
standard solutions prepared from the Lind-Evans radium stand- 
ard. The variation in observed alpha-count with radon content is 
linear over the lower range of values found in rocks (Fig. 1). On a 
basis of 5 grams of sample, the quantity normally used, the upper end 
of the graph corresponds to a radium content of about 0.6 x 10-7 g Ra/g 
of rock. With increase in radium content above this value, the reliability 
of the alpha-count method decreases rapidly due to the lack of resolution 
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of the individual counts when using a vacuum-tube electrometer and 
ballistic galvanometer. For this reason the alpha-count apparatus has 
a narrower useful range than the ionization apparatus, but this limita- 
tion is compensated by the somewhat greater accuracy of the former 
for low radium values. A sufficient number of check measurements 
have been made in the intermediate region of activities to assure good 
agreement between the two apparatuses. 


Measurements.—The observed radium contents of a number of rocks, 
expressed in the customary units of 10-’*g Ra/g, are arranged ac- 
cording to general lithologic types in five tables (Tables 1 to 5). Tables 
1, 2, and 3 contain the results on acidic, intermediate, and basic rocks; 
Tables 4 and 5, metamorphic and sedimentary rocks. Two samples, 30B 
and 35B included in Table 3, are of particular interest not only in re- 
flecting the unusually low radium content of ultrabasic rocks, but they 
also serve to indicate that the analytical techniques do not include ap- 
preciable contamination effects. Brief geologic notes for each rock speci- 
men are given in a later section. 


Accuracy.—An estimate of the over-all analytical probable error of an 
individual radium measurement can be obtained from a consideration of 
repeat measurements on uniform samples and on standard solutions used 
for calibration purposes. Using standard radium solutions for reference 
purposes, values of the calibration constant, Kp, obtained over a 2-year 
period, show a maximum spread in distribution of only 2.3 per cent from 
the mean (Evans et al., 1939, Table II, p. 937). These observations indi- 
cate: (1) the degree of reliability of the ionization measurements and 
(2) that the electrical properties of the detection apparatus have not 
changed significantly during the measurements. Duplicate analyses 
of homogeneous rock samples afford a check on: (1) the reliability of 
the direct-fusion method of releasing radon and (2) the efficacy of the 
transfer of radon from the furnace to the ionization chamber for measure- 
ment. The results for multiple analyses of rock samples 3A in Table 1, 
5B, 18B, 19B, and 28B in Table 3, and 9S in Table 5 are in agreement 
within 8 per cent of the mean in each case. Since these measurements 
were made on different portions of the same powdered sample of each 
rock, the interpretation given is that the parent element, uranium, and 
hence also its daughter element, radium, is fairly uniformly distributed 
throughout the hand specimens. After a rock is ground, different 5 to 
10 gram portions appear to contain equal radium contents within the ex- 
perimental uncertainty of measurement, which is estimated to be 8 per 
cent. At very low radium contents, the statistical fluctuations in the 
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TaBLeE 1.—Radium measurements of acidic igneous rocks 
Ra content in 10-! g/g rock 
Sample 
No. Individual Mean 
measurements value 
*4.1 +0.2 
3A Missourian granite.................. *3.17 +0.03 3.1 +0.1 
*3.14 +0.05 
*2.92 +0.03 
3.18 +0.03 
4A, Chelmsford granite................. 4.50 +0.05 2.5 +0.1 
4A, Chelmsford granite................. *2.48 +0.05 
*2.04 +0.03 
*3.10 +0.04 
*2.40 +0.03 
5A Rapakiwi .42 +0.0. +0. 
*2.20 +0.04 
1.96 +0.06 
2.81 +0.10 
SA Rockport xenolith.......:...:....0. 1.80 +0.06 1.80 +0.36 
d 9A 1.34 +0.03 1.34 +0.27 
Creighton 1.64 +0.04 1.25 +0.25 
10A, Creighton granite................... 1.20 +0.05 
10A; Creighton granite................... 0.95 +0.05 
0.80 +0.03 
11A Granite, Franklin, Me............... 1.24 +0.05 1.24 +0.25 
0.94 +0.03 0.94 +0.19 
*1.0 +0.05 
Fitchburg 0.88 +0.03 0.88 +0.18 
16A Granite, Kawatana, H.S............ 0.87 +0.01 0.87 +0.18 
17A | 0.41 +0.04 0.36 +0.05 
0.31 +0.03 
*0.21 +0.01 
*0.19 +0.04 0.20 +0.03 
*0.20 +0.03 
22A Rhyolite du Rosskopf............... 1.80 +0.07 1.80 +0.36 
0.93 +0.03 
24A Golden Gate rhyolite............... 1.23 +0.05 1.23 +0.25 
25A Porphyritic rhyolite................. 0.88 +0.04 0.88 +0.18 
26A Creighton aplite phase.............. 3.13 +0.07 3.13 +0.63 
1.36 +0.04 1.36 +0.27 
29A Sudbury micropegmatite............ 1.13 +0.02 1.01 +0.11 
0.88 +0.01 
Average for 29 acidic rocks:........ 1.64 +0.20 


* Sample stored for 30 days or longer before radon was removed for measurement. 
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emission of alpha-particles dictate the precision of the measurements 
(Evans and Neher, 1934, p. 144). 

Because of the indicated reliability for moderate radium values 
(0.1 10-** g/g or above), variations between duplicate analyses in 
excess of 10 per cent are attributed to an inhomogeneous distribution 


TaBLE 2.—Radium measurements of igneous rocks intermediate between acidic and basic 


Ra content in 10-! g/g rock 


No. ee Individual Mean 
measurements value 
1I Crescent Hill trachyte.............. 0.90 +0.03 0.90 +0.18 
21 Bunsen Peak dacite................. 0.60 +0.02 0.60 +0.12 
3I 0.52 +0.04 0.57 +0.04 
0.53 +0.04 
0.66 +0.05 
41 La Motte granodiorite.............. 0.57 40.04 0.57 +0.11 
5I Tholeiite, Cleveland dike............ 0.44 +0.03 0.44 +0.09 
61 Hornblende granodiorite............. 0.24 +0.04 0.29 +0.04 
0.33 +0.05 
71 Quartz dolerite, Whin Sill........... 0.17 +0.01 0.17 +0.04 
Average value of 7 intermediates:... ..............-. 0.51 +0.05 


of radium, and hence of its parent element, uranium, in the sample. 
Such inhomogeneity appears to be usual rather than exceptional for 
rocks and introduces a sampling probable error which may greatly 
exceed the analytical uncertainty. In the tabulation of data, the ana- 
lytical probable error, based on the residuals of hourly observations, is 
given after each individual measurement. The probable error of the 
mean of a number of individual measurements includes an observed, or 
estimated, field sampling uncertainty. In all single measurements, a 
sampling probable error of 20 per cent has been arbitrarily assigned when 
the observed value is listed as the mean for the particular rock body. 
For the few rock bodies that have been studied sufficiently to allow an 
estimate of their inhomogeneity, 20 per cent appears to represent the 
median variation in radioactive content. 

The probable errors of the average radioactivities indicated for the 
various rock types were calculated by the usual statistical methods from 
the residuals of the mean values. It is understood that average values, 
obtained by repeating the measurements on the same number of rock 
specimens selected at random, will have a 50 per cent chance of lying 
within the range set by the probable error of the present average values. 
In attempting to estimate the average radioactivity of the earth’s crust, 
much larger probable errors must be assigned to the various rock types, 
since such a small portion of the crust has been sampled and measured. 
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TABLE 3.—Radium measurements of basic igneous rocks 


Ra content in 10-! g/g rock 


Sample ‘ 
No. Individual Mean 
measurements value 
2B Gravenoire basalt... 0.64 +0.04 0.64 +0.13 
3B Crescent Hill basalt ................ 0.39 +0.01 0.39 +0.08 
4B 0.39 +0.03 0.39 +0.08 
5B Columbia River basalt.............. *0.325+0.009 0.31 +0.02 
*0.303+0.012 
0.305+0.007 
0.30 +0.03 
0.25 +0.01 0.21 +0.04 
0.16 +0.02 
7B Geode Creek basalt................. 0.20 +0.02 0.20 +0.04 
*0.15 +0.01 
0.20 +0.01 
OB, Watchung basalt.................0 0.18 +0.01 0.17 +0.03 
OB. Watchung basalt................... 0.16 +0.01 
10B 0.17 +0.02 0.17 +0.03 
11B 0.18 +0.02 0.15 +0.03 
*0.11 +0.01 
12B 0.14 +0.02 0.14 +0.03 
13B Oxbow Creek basalt................ 0.13 +0.02 0.13 +0.03 
Average value for 13 basalts:...... 0.30 +0.04 
14B Diabase (Horne Mine).............. 0.46 +0.03 0.46 +0.09 
‘ 15B, Sudbury olivine diabase............. 0.39 +0.04 0.34 +0.08 
15B, Sudbury olivine diabase............. 0.54 +0.03 
0.21 +0.02 
15B; Sudbury olivine diabase............. 0.26 +0.01 
0.24 +0.01 0.19 +0.02 
0.15 +0.01 
17B Somerset County diabase............ 0.21 +0.01 0.19 +0.02 
*0.17 +0.01 
18B Virginia Triassic diabase............ *0.17 +0.03 0.18 +0.01 
*0.17 +0.005 
0.186+0.004 
Average value for 5 diabase rocks:.. ..............-. 0.27 +0.03 
q 19B Worthington norite................. *0.85 +0.04 0.76 +0.11 
*0.70 +0.04 
*0.83 +0.05 
0.65 +0.04 
0.76 +0.05 
0.55 +0.02 0.55 +0.07 
0.54 +0.02 
21B er 0.25 +0.02 0.25 +0.05 
22B 0.14 +0.02 0.13 +0.01 
0.11 +0.03 
Avernpe value for 4 norites:. 0.42 +0.08 
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TaBLe 3—Radium measurements of basic igneous rocks—Continued 


Sample 


Ra content in 10—!? g/g rock 


No. Individual Mean 
measurements value 
0.82 +0.05 
0.47 +0.04 0.47 +0.08 
0.21 +0.01 
26B Riasg Buidhe monchiquite........... 0.68 +0.03 0.68 +0.14 
27B Kilchattan monchiquite............. 0.32 +0.02 0.32 +0.06 
28B Blueberry Mt. gabbro-diorite........ *0.168+0.005 0.16 +0.01 
*0.15 +0.007 
0.17 +0.01 
30B Dante, Weveter, N.C. *0.011+40.003 0.012+0.004 
*0.008+0 .002 
0.013+0.004 
*0.010+0.002 
*0.016+0.003 
31B Kimberlite, 5. Africa................ *0.62 +0.01 0.54 +0.06 
*0.54 +0.01 
*0.48 +0.02 
0.54 +0.01 
*0.53 +0.01 
32B Dalhousie andesite.................. 0.37 +0.03 0.37 +0.07 
34B Camptonite, Burlington, Vt.......... 0.085+0.008 0.077+0.018 
0.104+0.01 
0.041+0.007 
35B Wutchumna diorite................. 0.000 +0 .005 6.00 +0.01 
Average value for 35 basic rocks:... ................ 0.31 +0.03 
* Sample was stored for 30 days or longer before radon was removed for measurement. 
TABLE 4.—Radium measurements of metamorphic rocks 
Ra content in 10-! g/g rock 
Sample Rock 
No. Individual Mean 
measurements value 
Taroko-kyo gneiss.................. 0.77 +0.02 0.77 +0.15 
0.71 +0.03 0.71 +0.15 
3M __ Saastal eclogite..................... 0.077+0.007 0.07 +0.01 
0.072+0.006 


THORIUM MEASUREMENTS 

Techniques.—In the past, the thorium content of terrestrial materials 
has usually been determined from thoron measurements (for example: 
Evans, 1935, p. 110; Urry, 1936, p. 34). However, measurement of the 
total radioactive content from observations of the rate of alpha-particle 
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TABLE 5.—Radium measurements of sedimentary rocks 


Ra content in 10-1! g/g rock 


Sample 


No. Rock specimen Individual Mean 
measurements value 
18 Silurian bituminous limestone........ *2.46 +0.09 2.3 +0.3 
*2.20 +0.10 
3S:  Shaly limestone and sandstone....... 0.31 +0.03 0.30 +0.04 
3S,  Shaly limestone and sandstone....... 0.29 +0.03 
4S 0.29 +0.03 0.29 +0.06 
5S Carthage limestone................. *0.13 +0.01 0.11 +0.01 
*0.21 +0.02 
*0.08 +0.01 
0.12 +0.01 
*0.11 +0.01 
0.08 +0.01 
*0.08 +0.01 
*0.08 +0.01 
*0.13 +0.01 
Average value for 5 limestones:..... ..............-. 0.80 +0.20 
1.17 +0.06 
0.58 +0.05 
8S 0.67 +0.05 0.67 +0.15 
*0.23 +0.01 
*0.205+0.005 
*0.23 +0.008 
*0.21 +0.005 
Average value for 5 sandstones:.... 0.73 +0.12 
*0.047+0.003 
*0 .068+0.004 
0.084+0.005 
: Average value for 13 sedimentary rocks:............. 0.70 +0.12 


* Sample was stored for 30 days or longer before radon was removed for measurement. 


emission affords a simplified method of determining the thorium content 
of rocks, when the radium content is known. The equation relating 
these quantities is simply: 


‘ 
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Th = (Activity Index—104 (1) 
where Activity Index = the number of alpha particles (helium nuclei) emitted per 
hour per milligram of rock. 
Th= thorium content in 10° g/g of rock. 
Ra = radium content in 10™ g/g of rock. 


The derivation of Equation (1) and a detailed description of the pro- 
cedure used in the measurement of the activity index are included in the 
paper by Goodman and Evans (1941, p. 507). 

In the application of Equation (1), it is tacitly assumed that radio- 
active equilibrium exists. Since 1 million years suffices to establish 
equilibrium in the uranium series, and a shorter time is sufficient for 
the actino-uranium and thorium series, from a time standpoint this as- 
sumption is applicable to all but very young rocks. Recent alteration 
and weathering may preferentially remove certain series members and 
upset equilibrium, but this effect appears rather unlikely for nearly 
all the specimens measured, particularly the igneous rocks. 


Calibration—The alpha activity measurements are not dependent 
upon radioactive standards. (Goodman and Evans, 1941, p. 507). 
Hence, frequent recalibration, as is routine in radium measurements, is 
not customary in this technique. Some internal standardization meas- 
urements have been made (Finney and Evans, 1935, p. 511) using rock 
powders containing added known amounts of chemically analyzed radio- 
active ores. The results of these observations are in good agreement 
with the theoretical predictions for the alpha-counting method. How- 
ever, the most convincing evidence of reliability is the concurrence of 
activity index measurements based on separate radium and thorium 
analyses of powdered rock samples also measured for alpha activity 
(Evans et al., 1939, Table VI, p. 939). A number of similar measure- 
ments have subsequently been made on the Palisade diabase, the results 
of which are included in Table 8. 

These intercheck thorium analyses are based on measurements of 
thoron by the streaming method? and depend upon calibration with 
standard thorium solutions. No official thorium standards have ever 
been established.2 In the past, the customary procedure has been to 
prepare reference standards from chemically analyzed thorium minerals, 
usually a monazite, thorianite, or thorite. As a temporary standard, 
the authors have adopted a solution prepared from an analyzed thorian- 


1 The method used is similar to that described by Urry (1936, p. 35) but with a number of changes 
in the detection apparatus and improvements in the chemical techniques. 

2The National Research Council Committee on Standards of Radioactivity is endeavorirg to fill 
this need. Through the co-operation of Drs. A. F. Kovarik and J. P. Marble, work is now in 
progress on the preparation of a primary thorium standard, which it is hoped will later be adopted 
as the primary international thorium standard. 
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ite furnished by the late Professor Herman Schlundt. This solution is 
considered to contain 4.18 X 10-* g Th/ce and 2.1 X 10-* g U/ce, with 
0.7 per cent HNO, for protective purposes. Comparison with thorium 
standards prepared from other analyzed thorium minerals (Evans et al., 
1939, p. 939) indicate agreement within the experimental uncertainty 
of a few per cent. 

Measurements.—In Tables 6, 7, and 8 are listed the individual measure- 
ments and resulting mean values for activity index determinations of 
a number of rock specimens,* most of which also have been analyzed for 
radium as reported in Tables 1, 2, and 3. The uranium content of a 
given specimen is calculated from the mean radium value using the 
conversion: U = 2.84 10° Ra. The calculated thorium content results 
from substitution of the mean activity index and mean radium value in 
Equation (1). In the last column of these tables, the ratios of thorium 
to uranium have been listed. The Th/U ratio is an index of the relative 
activity of the two series in the production of helium, radiogenic lead, 
and heat in rocks as discussed in a later section of this paper and by 
Goodman and Evans (1941, p. 506). 

Accuracy.—The remarks, concerning inhomogeneity, field sampling 
uncertainty, and calculation of probable errors, included in the discus- 
sion of radium measurements, are equally applicable to thorium measure- 


ments. In addition, the small amount of sample used (less than 100 milli- 
grams) in an individual alpha activity determination increases the ana- 
lytical uncertainty to about 20 per cent for most samples. In order to 
decrease this uncertainty to the point where it is comparable to that of the 
radium analyses, several alpha-count measurements are usually made on 
each sample. 


DISCUSSION OF RESULTS 

General statement.—In considering the radioactivity measurements of 
rocks presented in Tables 1 to 8, the reader may justly ask: Do these 
data merely add to the several hundred other results on similar materials 
obtained by Strutt, Joly, Poole, Piggot, Fletcher, Urry, Evans, and the 
numerous other investigators of the past 30 years, and hence merely define 
our knowledge of terrestrial radioactivity a little more sharply? The 
answer to this inquiry is that the present measurements have consider- 
able more significance in that they constitute the most complete set of 
well-authenticated results which are now available. An ever increasing 
mass of evidence indicates that few of the earlier measurements are 
sufficiently reliable to allow other than the broadest of geological deduc- 
tions. Whereas individual determinations have been reported as having 


3 All the activity indices were determined from alpha-counting measurements, except for some 
separate radium and thorium measurements made on samples 16B and 18B as reported in Table 8. 
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TaBLE 6.—Radioactivity measurements of acidic ig 
Activity index in alphas/hr 


Rock specimen 


Rockport xenolith 


Mt. Manitou granite 
Creighton granite 


eee 
“28 


Creighton ite 
Granite, Glen Cove, Colorado 


4. 
3 
3 
2. 
3. 
2 
2 
2 
2 
3. 
3 
2. 
2. 
2 
2. 
2 
2 
2 
2 
1 
0 
0 


Number of 


Golden Gate rhyolite 
Obsidian Cliff 
Pegmatite, Ruggles Mine 


* Radium measurements by Evans and Raitt (1935, p. 172), who have also previously reported 


| 
measurements 
Granite; Freakin, Maina... 1.24 
21A Granite, Winnipeg, Manitoba.................... 0.36 
0.48 
0.56 
0.59 


rements of acidic igneous rocks 


ivity index in alphas /hremg 


Uin Thin Th/U 
/g rock 10-°g/g rock ratio 
+0.09 +0.5 6.5 +0.9 17, +4. 2.6 
t +0.08 
+0.1 3.2 40.4 +0.5 20. +4. 5.4 
+0.3 2.8 +0.3 3.5 +0.7 17. +4. 4.8 
+0.1 
+0.2 
+0.06 2.8 +0.3 6.5 +1.4 4.3 +1.0 0.7 
+0.07 
+0.1 2.8 +0.3 2.5 +0.5 21. +85. 8.4 
+0.2 
+0.2 
+0.1 
+0.06 
+0.05 
+0.06 
+0.06 
8B +£0.98 2.1 5.1 +1.0 2.6 +0.5 0.5 
+0.06 
+0.10 1.7 +0.3 2.7 +0.5 8.1 +2.0 3.0 
+0.10 
+0.08 
+0.07 
ee 1.24 +0.16 3.5 +0.7 0.0 +0.5 0.0 
ae ed 0.36 +0.07 0.34 +0.06 2.7 +0.8 7.8 
2 a 2.3 +0.2 3.8 +0.4 10.3 +1.7 3.7 
(11) (9) (9) 
+0.1 3.6 +0.8 4.6 +0.9 22 +5. 4.7 
+0.1 3.5 +0.8 3.5 +0.7 25 +5. 7.2 
+0.1 3.3 +0.7 3.9 +0.8 21 +5. 5.4 
+0.03 
+0.04 
D +0.04 
2.4 +0.2 3.9 +0.3 13.4 +1.6 4.2 
2) ee (15) (12) (12) (12) 


reviously reported radium measurements on 6A and 21A, 


AS 

23% 


of acidic igneous rocks 


index in alphas/hremg 


U in Th in Th/U 
10-*g/g rock 10-*g/g rock ratio 
+£0.09 +05 6.5 +0.9 17, 2.6 
+0.08 
+0.1 3.2 40.4 %3.8 +0.5 20. +4. 5.4 
+0.3 2.8 +0.3 3.5 +0.7 4.8 
+0.1 
+0.2 
+£0.06 2.8 +0.3 6.5 +1.4 4.3 +1.0 0.7 
+£0.07 
+0.1 2.8 +0.3 2.5 40.5 21. x85. 8.4 
0.2 
+0.2 
+0.1 
+£0.06 
B +0.05 
+0.06 
+£0.06 
B +£0.08 2.1 5.1 +1.0 2.6 +05 0.5 
y +0.06 
41.7 2.7 8.1 +2.0 3.0 
+0.10 
+0.08 
+£0.07 
1.24 40.16 3.5 +0.7 0.0 +0.5 0.0 
0.36 40.07 0.34 +0.06 2.7 +0.8 7.8 
2.3 +0.2 3.8 +0.4 10.38 +1.7 3.7 
(11) (9) (9) (9) 
40.1 3.6 408 46 40.9 22. <5. 4.7 
+0.1 3.5 +0.8 3.5 +0.7 25. +65. 7.2 
+0.1 3.3 +0.7 3.9 +0.8 21 +5. 5.4 
B +0.03 
B 0.04 
+0.04 
sense 2.4 +0.2 3.9 +0.3 13.4 +1.6 4.2 
ee (15) (12) (12) (12) 


reviously reported radium measurements on 6A and 21A. 
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an analytical uncertainty of 5 per cent or less, intercomparative analyses 
between different laboratories may differ by factors of 2 to 100. Itisa 
painful process to discard the results of years of research, and we are 
naturally reluctant to do so in the light of anything other than the most 
convincing proof of unreliability. However, in view of the extensive 
efforts that have been made to avoid all known sources of error in the 
present work, it seems justifiable to question results by other laboratories 
when a difference is found to exist, particularly if it is known that im- 
portant analytical precautions have not been taken in these laboratories. 


Comparison with former measurements——A direct comparison of ter- 
restrial radioactivity measurements between laboratories can only be 
made on identical samples, preferably powdered samples in order to 
reduce the uncertainty due to inhomogeneity. If a difference is observed, 
a comparison of standards, in the case of radium and thorium measure- 
ments, will aid in isolating the difficulty. An indirect and much less satis- 
factory comparison between laboratories can be made from a consideration 
of the average values obtained for the various rock types. This method 
requires a large number of measurements, and its accuracy is limited by 
the broad range of radioactivities observed for a given rock type, as is 
evident from the spread of values given in the tables. For evaluating 
many of the earlier measurements reported in the literature, only the in- 
direct method is possible, since the samples studied are no longer avail- 
able. Among investigators still active in the field, a systematic program 
of interchecking has been in progress for several years, and a number of 
definite results can be stated at this time. 

A co-operative investigation, with two independent methods of meas- 
urement, has shown (Evans e¢ al., 1939, p. 931) that all radium deter- 
mination published by Urry between July 1934 and January 1939* are 
low by a factor of over two. Between 30 and 40 per cent of this error is 
explicable by unreliable radium standards, but the remaining source of 
discrepancy is of an unexplained origin. Because of this uncertainty it 
does not appear possible to correct Urry’s radium measurements by a 
constant factor. Urry’s thorium measurements are also open to slight 
question, because of analytical uncertainties. Until further work is done 
to validate these results, it does not appear sound to include them in tabu- 
lations of better authenticated radioactivity measurements. For these 
reasons, extensive research, which otherwise would be an important sup- 
plement to our knowledge of terrestrial radioactivity, must be discarded, 
or at least held in abeyance. 


*These publications include Lane and Urry (1935, p. 1101-1120), Lane (1935-1938) and Urry (1936a, 
p. 40-48; 1936b, p. 1218). 
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An exchange of five powdered rock samples (2A, 14A, 20A, 36A, and 
18) with F. Behounek ° has clearly established the unreliability of his 
radium measurements on rocks. Unsystematic errors, ranging up to a 
factor of 30, were found by this intercheck. Evidence has previously 
been presented (Evans, 1934, p. 36, 40) which strongly suggests that all 
of Graven’s (1930, p. 181) and Graven and Kirsch’s (1932, p. 521) terres- 
trial radioactivity measurements are without foundation. It can also be 
shown by a simple calculation that the thorium technique described by 
Hirschi (1920, p. 545) cannot have a sensitivity sufficient to give signifi- 
cant values. Five hand specimens (5A, 2B, 11B, 17B, and 3M), previously 
analyzed by Joly and Poole, very kindly have been furnished by Professor 
Poole for intercheck measurements. The radium measurements on these 
rocks, reported in Tables 1, 3, and 4, are invariably lower than the earlier 
values of Joly and Poole, given in the geologic notes. It will be seen by a 
comparison of the two sets of values for these five specimens that no fixed 
ratio obtains. However, Joly and Poole’s results show a greater ratio to 
the present values for the specimens with lower radium content. This 
Observation suggests the possibility of a systematic error in the earlier 
work which may be due to radioactive contamination in the chemicals 
used by Joly and Poole. The recent independent agreement on the radio- 
activity of sea water between Fgyn and his coworkers (1939, p. 10) and 
Evans et al. (1938, p. 243) also furnishes additional evidence that Joly’s 
radium values, on all types of specimens, are systematically too high. 
It appears that the measurements of a number of other workers (Evans 
et al., 1938, p. 248, Table I) are also subject to question as the result of 
these researches. The present international situation has interfered 
with further interchecking with European laboratories, but it is hoped 
that in the not too distant future these differences can be studied in greater 
detail. 

An exchange of a number of rock specimens has also been made with 
the Geophysical Laboratory in Washington, D. C. Radium measure- 
ments on four rock specimens (2A, 14A, 20A, and 36B) show ratios be- 
tween the earlier values of Piggot, given in the geologic notes, and the 
present values, reported in Tables 1 and 3, of 1.4 to 2.8, with an average 
f of very close to 2. Intercheck measurements on the standard solutions 
used for calibration are in good agreement, indicating that these solutions 
are not the source of the differences. More recent radium measurements 
by Urry, working in Piggot’s laboratory, using the 12 standard rock 
samples (see 3A in geologic notes) for comparative measurements, are in 
quite good agreement with the values reported for these samples in Tables 


5¥Formerly Director of the Radiological Institute of Czechoslovakia. 
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1, 3, and 5. Further work is in progress to resolve the disagreement on 
the earlier specimens. 

On the other hand sufficient intercheck measurements have been made 
on the same rock samples among Evans, Clarke, Finney, Goodman, Keevil, 
and Raitt to establish the mutual agreement of the results summarized in 
Table 9° and those presented in the preceding tables. These data are 
combined in Table 10 to give more representative average values for the 
different rock types. 

A number of summaries of the radioactivity of terrestrial materials have 
previously appeared,’ but the emphasis has invariably been on the geologi- 
cal implications with little or no consideration of the probable reliability 
of the measurements themselves. As a result many of the geological con- 
clusions suffer from the lack of critical evaluation of the data on which 
they are based. 

Jeffreys’ (1936, p. 149-170) compilation of the radioactivities of igneous 
rocks is the most complete summary available. The measurements of 
Evans and his coworkers, given in Table 9, and omitted from Jeffreys’ 
summary, would lower Jeffreys’ averages to a slight extent. Other data 
not included are those of Urry, which were omitted because they appeared 
anomalously low compared to other observations. As indicated above, 
we now know that these results are erroneously low. Intercheck measure- 
ments with other workers whose results are included in Jeffreys’ averages 
suggest that where other errors exist they are very likely in the direction 
to give unreliably high values. For these reasons a comparison of the 
weighted average values of Table 10 with Jeffreys’ average values gives 
a rather reliable indication of the downward revision that appears neces- 
sary in estimates of the radioactive content of igneous rocks (Table 11). 
Too few specimens of intermediate rocks have been measured to give 
conclusive results. However, since a definite decrease in the average 
radioactivities is indicated for both the acidic and basic rocks, a similar 
suggested change in the intermediates lends support to the significance 
of the observed differences. A decrease of 40 per cent in the estimated 
average radium content of acidics and a 60 per cent decrease for basics 
are indicated in Table 11. Decreases of 25 and 50 per cent in the esti- 
mated average thorium contents are also indicated for acidics and basics 
respectively. The decreases in both the radium and thorium values are 


®6In addition to the results included in Table 9, Keevil (1938a, p. 146; 1938b, p. 690-691) has pub- 
lished a number of other measurements. These analyses were made on rock speci also m 
by the authors and included in Tables 1 to 8. Because inclusion of Keevil’s duplicate analyses in 
Table 9 would give undue weight to certain specimens in the averages of Table 10, his results have 
been omitted. 

7 Among others are: Holmes (1914, p. 12; 1915, p. 2; 1926, p. 1055; 1931, p. 559; 1937, p. 763), 
Hoffmann (1939, p. 189), Hatuda (1936, p. 1), Jeffreys (1936, p. 149), Poole and Joly (1924, p. 819), 
Piggot (1931, p. 28), and Urry (1933, p. 125). 
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hardly surprising since Jeffreys’ averages are strongly dependent upon 
the measurements of Joly and Poole, Graven and Kirsch, and Hirschi, all 
of which appear to be on the high side. 


Possible sources of error —One interpretation of the larger percentage 
decrease indicated for basic rocks as compared with acidic rocks is that 
the mean errors in the earlier radium and thorium determinations were 
fairly constant. Under these conditions, the error would contribute pro- 
portionately less to the more radioactive (and, therefore, more readily 
measurable) acidic rocks than to the less radioactive basic rocks. Con- 
tamination errors are of this nature. If a reliable correction factor could 
be found, it might be possible to salvage a significant fraction of the 
earlier work in this field. 

The indicated increase in Th/U ratio may be the result of proportion- 
ately larger contamination errors in the earlier radium measurements 
than in those of thorium. As yet, there is no definite proof of this point, 
but Fgyn et al. (1939, p. 40) have observed a greater likelihood of radium 
contamination than uranium and thorium in ordinary chemical reagents. 


General discussion.—It was observed quite early (Strutt, 1905, p. 472) 
that the concentration of radioactive elements in rocks increases with 
silica content. A number of other investigations in addition to the present 
work have amply confirmed this general conclusion. A correlation of 
radium with potassium content has been reported (Evans and Williams, 
1935, p. 447; Hatuda, 1936, p. 15). A notable exception to this latter ob- 
servation is the lack of correlation with potassium but the apparent cor- 
relation with biotite content observed by Barth (1938, p. 240) in the 
Finnish granites measured by Piggot (1938, p. 227). The principal diffi- 
culty with most of these investigations is that the radioactivity measure- 
ments have been made on the composite rock instead of on the separate 
mineral constituents. In addition, two aspects of this problem seem to 
have escaped most workers in this field. The first, and rather obvious 
consideration, is that, while radium is the substance usually measured, 
the parent element, uranium, must be the one concentrated in the rock 
at the time of its solidification in order for radium and the other elements 
in the series to be present today, except for very recently formed rock 
bodies. The second point is that there is no particular reason to expect 
the same correlation to apply to rocks of totally different genesis, even 
though they may be of the same general chemical composition. Certainly 
the conditions of crystallization were entirely different for the Lassen 
lavas (Evans and Williams, 1935, p. 441) as compared with those that 
obtained when the Finnish granites (Barth, 1938, p. 231) were formed. 
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Since it now appears necessary to base helium age measurements on 
selected mineral constituents of rocks rather than on the rocks as a whole 
(Hurley and Goodman, 1941, p. 553), an increasing number of radio- 
activity measurements on separate minerals will probably appear as an 
important by-product of age researches in the future. These results 
should be of considerable value in solving the geochemical problems of 
the distribution of radioactive elements in rocks. 


RADIOACTIVE HEAT 
GENERAL STATEMENT 


The radiations emitted by the radioactive elements contained in the 
earth are absorbed by the rocks, and the energy of disintegration appears 
in the form of heat. Since radioactivity is a nuclear property of atoms, 
the rate of decay would not be expected to be affected by changes of pres- 
sure, temperature, or chemical composition—even such high temperatures 
and pressures as must exist deep in the earth. There is considerable ex- 
perimental evidence to support this conclusion. 

The energies of the three types of radiation emitted by the members 
of the uranium, actino-uranium, and thorium series have been accurately 
measured. From these data and the disintegration constants of the parent 
radioactive elements, the specific rate of generation of heat for each 
radioactive series (in equilibrium) can be computed as shown below. 


DERIVATION OF CONVERSION FACTORS 


The energy accompanying the emission of alpha radiation constitutes 
the source of more than 85 per cent of the heat generated by the naturally 
occurring radioactive elements. The alpha-ray energies for the UI and 
Th series, including the correction for the recoil nucleus, are summarized 
in Table 12. Combining the total values for these series with the follow- 
ing accepted physical constants.*® 
Aur = 4.87 & 107 sec™, Arn = 1.58 sec™* 

Atomic weights: U = 238.14, Th = 232.12 


Avogadro’s number: N = 6.02 10” atoms per gram equivalent 
1 calorie = 4.19 & 10° ergs, 1 MEV =1.60 X 10° ergs, 1 year = 3.15 & 10° sec. 


the heating effect for the radiations can be calculated as follows: 


Ratio (Rutherford et al., 1930, p. 161) a: 8 : y= 1248 :63:9.4 


hence: a + 6B + oo 140.5 


1048 = 1.125 


S\ur from Kovarik and Adams (1941) including Nier’s (1939, p. 150) value for the actino- 
uranium /UI activity ratio. Arn from Kovarik and Adams (1938, p. 413). 
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Therefore, heat per gram UI per year = 
43.34 & 1.125 & 1.60 x 10° X 3.15 & 10° N Av — 0.708 calori 


Heat per gram U per year = 0.741 calories’ 


Similarly, heat per gram Th per year = 
36.82 X 1.125 X 1.60 X 10° X 3.15 X10°NdAm _ 0.203 calori 
4.19 X 10° Th 


TABLE 12.—Uranium and thorium series alpha-ray energies 


Mean range in cm Disintegration energy 
Element In air at 15°C in million electron 

760 mm Hg volts 
3.09 4.66 
4.05 5.58 
4.66 6.11 
3.84 5.41 
3.95 5.47 
5.64 6.90 
35% Th C (av.)..... 4.73 2.16 


The heating effect of potassium is estimated from Miilhoff’s (1930, 
p. 205) determination of the specific activity as 283(+ 15 per cent) beta 
rays per second per gram of K and values of 0.7 MEV for the maximum 
beta ray continuous spectrum (Anderson and Neddermeyer, 1934, p. 653) 
and 2 MEV for the gamma ray energies (Miilhoff, 1930, p. 205; Gray and 
Tarrant, 1934, p. 694). The latter investigators have shown that three 
gamma ray quanta are emitted per 100 disintegration electrons. Combin- 
ing these data, the heat per gram of K per year = 4.6 X 10~-° calories. 
Taking into account the uncertainties in the physical quantities, the 
rate of heat generation for the three most abundant naturally occurring 
radioactive substances can be summarized as: 
(UI + AcU) Series: 0.74 + 0.03 cal/yr/gm U 


Th Series: 0.20 + 0.01 cal/yr/gm Th 
Potassium: 5 + 2 & 10° cal/yr/gm K 


° Addition of 4.6 per cent for AcU series (Nier, 1939, p. 150). 
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HEAT GENERATION IN ROCKS 


Using the above conversion factors and the average radio-active con- 
tent of rocks from Table 10, an estimate of the rate of generation of 
heat in the different types of rocks can be made. The resulting values, 
including the relatively small contribution from potassium,’® are com- 
pared with a similar earlier estimate by Holmes (1930, p. 136; 1931, p. 
559) in Table 13. 


TABLE 13.—Estimates of radioactive heat generation in rocks 


Heat production in 10-* calorie per gram rock per year 


Acidic Intermediate Basic 


igneous igneous igneous Sedimentary 
(Holmes’ values*) 

15.9 13.3 5.7 
(Present values) 

2.2 +0.2 1.0 +0.15 0.7 +0.08 0.9 +0.09 
2.6 +0.4 0.9 +0.2 0.8 +0.01 0.7 +0.1 
_ 0.14+0.02 0.10+0.02 0.07+40.01 0.10+0.02 

| ee 5.0 +0.6 2.0 +0.4 1.6 +0.1 1.7 +0.2 


* Holmes’ classification of igneous rocks into the three broad types differs somewhat from that used 
in the present discussion. Holmes’ values have been rearranged according to the present grouping of 
rocks as indicated in Tables 1 to 8. 


The extent to which Holmes’ values have been generally adopted is re- 
flected by their inclusion, in part or as a whole, in numerous publica- 
tions (Daly, 1933, p. 70; Grout, 1932, p. 146; Ellsworth, 1932, p. 65; 
Hevesy and Paneth, 1938, p. 274). The estimates based on the present 
values for the radioactive content and the conversion factors are con- 
siderably lower than those of Holmes. These differences are largely 
attributable to (1) the lower estimates for the radioactivities and (2) a 
decrease in the estimated geothermal importance of potassium as the 
result of more precise studies of the radioactive properties of this element 
in recent years. 

CONCLUSIONS 


Several hundred radioactivity measurements have been made on rocks 
using unusual analytical care. The results of these measurements have 
been combined with the limited number of well-authenticated results 


10The potassium content has been estimated from average chemical analyses reported by Daly 
(1933, p. 9-28) and Grout (1932, p. 127) to be: 2.8+ 0.5 X 10-2g/g for acidic igneous rocks, 2.0 + 0.4 X 
10-°2/g for intermediate igneous rocks, 1.4 + 0.3 X 10-°g/g for basic igneous rocks, and 2.0+ 0.4 X 
10-*g/g for sedimentary rocks. 
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previously reported by other investigators to obtain estimated average 
values for broad lithologic types. When compared with previous aver- 
ages, the present values indicate that a substantial downward revision 
is necessary in the present conception of the radioactive content of 
igneous rocks. Geological considerations based on the earlier values will 
very likely need revision. In particular Holmes’ estimates of the rate 
of generation of radioactive heat within the earth may need to be sub- 
stantially decreased. 


GEOLOGIC NOTES ON ROCK SPECIMENS REFERRED TO IN TABLES 
Aciwic Rocks 


1A Brifosse granite obtained together with 22A from E. Rothe for intercheck 
radioactivity measurements. Gamma-ray activity of these rocks reported by 
Rothe and Hopcewicz (1937, p. 165) and by Rothe and Hee (1937, p. 1835). 


2A Kumlinge granite, Aland Archipelago, submitted as a powder by C. 8S. Piggot, 
one of four specimens (others: 14A, 20A, and 36B) for intercheck measure- 
ments. C. S. P.’s mean value. 5.51 & 10° g Ra/g, and geologic notes previously 
reported (Piggot, 1938, p. 228) using designation P-75. 


3A Missourian granite 50-pound specimen collected by U. 8. Geological Survey 
from the Sheahan Quarry, in about Section 14, T. 34.N., R. 3E., Graniteville, 
Missouri. One of 12 rocks (others: 4As, 19A, 5B, 8B, 18B, 28B, 30B, 31B, 5S, 
9S, and 118) being used as standard samples for the international intercalibra- 
tion of radioactivity measurements (Goodman and Evans, 1940, p. 557). This 
specimen was erroneously designated (Curtiss et al., 1940, p. 602) as “Carthage” 
granite. 


4A, Chelmsford granite collected by N. B. Keevil from H. E. Fletcher Co. quarry, 
N. Chelmsford, Mass. Age: possibly Mississippian (Emerson, 1917). 

4A, Chelmsford granite, 50-pound specimen collected by C. Goodman from same 
location as 4A:. One of standard samples; see 3A. 


5A Rapakiwi granite from Viborg, Finland, one of five hand specimens (others: 
2B, 11B, 17B, and 3M) submitted for intercheck radium measurements by 
J. H. J. Poole. These specimens are particularly valuable since they not only 
furnish an intercheck with the recent work of Poole (1938, p. 595) but also 
afford a means of evaluating the earlier results of the late J. Joly, since some 
of the rocks were measured by the latter. Poole (personal communication, 
1938) reports Joly’s value as 5.77 X 10°” g Ra/g for 5A. 


6A Granite collected by R. A. Millikan from boulder in Glen Cove, Colorado. 
See Evans and Raitt (1935, p. 172) and Raitt (1935) for other radioactivity 
measurements, designation C-13. 


7A Rockport granite, Lat. 42° 40’ N, Long. 70° 35’ W., collected by N. B. Keevil 

from Upper Pigeon Hill Quarry, Cape Ann, Mass. Intrudes Salem gabbro 

diorite, which in turn cuts pre-Cambrian rocks and is possibly of Devonian 

age. It is cut by diabase dikes interbedded with Triassic sedimentary beds 

= Deerfield and Mt. Holyoke, Mass. Age: Probably pre-Triassic and post- 
ambrian. 


8A Rockport xenolith, fine-grained and cognate, from the same large hand speci- 
men as 7A. 


9A Murray granite collected by H. C. Cooke, 1939, from railway cut, 4 to 5 feet 
from surface, lot 10, range V, McKim Township, Sudbury, Ont. 


10A: Creighton granite collected by H. C. Cooke, 1939, from road cut 5 feet from 
surface, SW quarter, lot 9, concession 1, Snider Township, Sudbury, Ont. 
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Creighton granite collected by E. A. Collins from surface rock SE of Mud Lake, 
Copper Cliff, Ont. Considered to be outcrop of 10As. 

Creighton granite collected by E. A. Collins from the Creighton Mine, 40 
level, 14 X-cut, Copper Cliff, Ont. 

Granite collected by N. B. Keevil from quarry in Franklin, Me.; see Dale 
(1923, p. 219). Age: Middle Paleozoic. 

Quincy granite collected by N. B. Keevil from J. S. Swingle quarry, Quincy, 
Mass. Age: Possibly Mississippian, probably late Paleozoic. 

Quincy granite drill core from over 600 feet deep, same location as 12Ax, ob- 
tained from R. A. Daly. 

Granite from collection of G. D. Louderback. Late Jurassic from the Sierra 
Nevada, Mokelumne River region, quarry for dam. 

Granite collected by J. J. Sederholm and T. F. W. Barth in summer of 1932 
from Ava, north of Brando, Aland Archipelago. Submitted as a small hand 
sample by C. S. Piggot for intercheck measurements; see 2A. Piggot’s mean 
value, 1.21 10™ g Ra/g, and geologic notes previously reported (Piggot, 
1938, p. 229) using designation P-83. 

Fitchburg granite collected by N. B. Keevil from Rollstone Hill quarry, Fitch- 
burg, Mass. Age: probably Carboniferous. 

Granite from Kawatana, H. S., Japan submitted by Z. Hatuda together with 
three other specimens (18A, 27A, and 1M) for intercheck radioactivity 
measurements. 

Granite collected by N. B. Keevil from Barryfield quarries, 2 miles east of 
Kingston, Ont. Age: post-Grenville and pre-Keweenawan. 

Granite from Azi, collection of Z. Hatuda; see 16A. 


Granite, 50-pound specimen obtained by C. Goodman indirectly from H. E. 
Fletcher Co. quarry in Milford, Mass. One of standard rock samples; see 3A. 
Compact, massive rock, somewhat above medium grain and of light color. 
Its two especially characteristic constituents constantly present are blue quartz 
and a microcline-microperthite. Age: post-Cambrian and pre-Carboniferous. 
Granite from Ounasjok in Rovaniemi, end of the railway in northern Finland, 
collected by J. J. Sederholm and T. F. W. Barth in the summer of 1932. Sub- 
mitted as a small hand specimen by C. S. Piggot for intercheck measurements; 
see 2A. Piggot’s mean value. 0.43 & 10°" g Ra/g, and geologic notes previously 
reported (Piggot, 1938, p. 229) using designation P-86. 

Granite collected by R. A. Millikan from outcrop in Red River, Winnipeg, 
Manitoba. See Evans and Raitt (1935, p. 172) and Raitt (1935) for other 
radioactivity measurements, designation C-8. 

Rhyolite du Rosskopf obtained from E. Rothe for intercheck radioactivity 
measurements; see 1A. 

Rhyolite collected by A. D. Howard 100 feet north of Tower Creek Bridge, 
Yellowstone Park. 

Rhyolite collected by C. Max Bauer from Golden Gate, Yellowstone Park. 
Porphyritic rhyolite collected by C. Max Bauer from near 8-mile Bridge, 
Yellowstone Park. 


Late aplite phase of Creighton granite, 10A, collected by H. C. Cooke, 1939, 
from 4 feet below original surface, lot 1, concession III, Snider Township, Sud- 


bury, Ont. 
Pitchstone from Aso volcano, collection of Z. Hatuda, see 16A. 
Obsidian collected by C. Max Bauer from Obsidian Cliff, Yellowstone Park. 


Micropegmatite collected by H. C. Cooke, 1939, from railway cut 6 feet 
below surface, lot 1, concession III, Wisner Township, Sudbury, Ont. 
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Granite from Mt. Manitou. Colo.. collection of R. A. Millikan. See Evans 
and Raitt (1935, p. 172) and Raitt (1935) for other radioactivity measurements, 
designation C-9. 


Lebanon granite collected by E. P. Kaiser (1938, p. 107) near Hanover, N. H. 
Granite collected by N. B. Keevil of unknown location, designation K-G-1. 
Pegmatite sample collected by A. C. Lane, 1937, from Ruggles Mine, Grafton, 


aN. 
Basic Icngzous Rocks 


Porphyritic basalt collected by A. D. Howard from SE base of Cresent Hill, 
Yellowstone Park. 


Vesicular basalt from Gravenoire volcano, Clermont, France, submitted as a 
small hand specimen for intercheck radioactivity measurements by J. H. J. 
Poole who reports 1.19 X 10°” g Ra/g as Joly’s value (personal communication, 
1938); see 5A. 


Cresent Hill basalt collected by A. D. Howard % mile west of crossing of 
Geode Creek, Yellowstone Park. 


Basalt from the well-known Giant’s Causeway in Ireland, submitted as a 
large hand specimen for intercheck measurements by Arthur Holmes (1937, p. 
155); also see 5I and 71. 


Columbia River basalt, 50-pound specimen secured through A. W. Fahrenwald. 
One of standard samples; see 3A. Composed of calcic labradorite, augite, 
olivine, and a little magnetite. Age: Middle Miocene. 


Basalt from Lat. 40° 40’ N., Long. 40° 45’ W., collected from a flow by A. C. 
Hawkins, 1936, 50 feet from the surface in a working quarry at Oldwick, N. J. 
The top of the flow has been removed by erosion, and the bottom is not 
visible. Similar flows in the area, 16B and 17B, of slightly different periods 
of eruptivity have been identified as being of Triassic age from fossils in 
adjoining sedimentary beds (red shales and sandstones). Age: Probably 
Triassic. 

Basalt collected by A. D. Howard *4 mile E of crossing of Geode Creek, Yellow- 
stone Park. 


Deccan Trap Plateau basalt, 50-pound specimen obtained from near Bombay, 
India, through the courtesy of the U.S. Geological Survey and A. S. Kalapesi, 
St. Xaviers College. One of standard samples; see 3A. Composed of plagio- 
clase (calcic andesine), augite (roughly 50 per cent), less magnetite, and a 
small amount of secondary minerals including chloro-phalite-delessite. 


Second Watchung trap collected Dec. 7, 1939, by P. M. Hurley from old quarry 
on Murray Hill short distance SW of Summit, N. J., on top of Second 
Watchung Mtn. Specimen from 10 feet above base of a flow at least 30 feet 
thick and taken about 25 feet below original surface. Rock has fine diabasic 
texture. 

First Watchung trap collected Dec. 7, 1939, by P. M. Hurley from fresh blast 
in operating quarry on First Watchung Mtn. opposite Summit, N. J., about 1 
mile east of Summit. Specimen from block of fresh broken rock probably 
coming from a point in the face between 30 and 50 feet below the original 
surface. No sign of limonite on any of the material. 


Basalt collected by A. D. Howard from Overhanging Cliff, Yellowstone Park. 
Basaltic lava from Hawaiian Islands submitted as a small hand specimen for 
intercheck radioactivity measurements by J. H. J. Poole who reports 
1.25 < 10°" g Ra/g as Joly’s value (personal communication, 1938); see 5A. 
Basalt collected by C. Max Bauer from near 8-mile Bridge, Gardiner, Yellow- 
stone Park. 

Oxbow Creek basalt collected by A. D. Howard % mile west of crossing of 
Oxbow Creek, Yellowstone Park. 
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Diabase from Lat. 43° 15’ N., Long. 79° 2’ W. collected by P. Price 3 feet from 
chilled contact of a N-S dike at the 17th level (1975 ft.) of the Horne Mine, 
Norando, Quebec. In a previous publication (Evans et al., 1939, p. 946) this 
rock was indicated as post-Cobalt in age. Subsequently, a personal communi- 
cation from H. C. Cooke informs us that this dating is very uncertain and 
the most that can be said is that this rock is post-Timiskaming. (See 
Geological Survey of Canada Memoir 166, p. 143-144.) 


Olivine diabase (coarse-grained phase) collected by E. A. Collins from hanging 
wall of pit at fault, Creighton Mine, Copper Cliff, Ont. Considered to be 
outcrop of 15Bz. 


Olivine diabase collected by E. A. Collins from the Creighton Mine, 20 level, 
Main No. Dr., Copper Cliff, Ont. 


Olivine diabase collected by H. C. Cooke, 1939, from about 2 feet below original 
surface, west side of lot 10, concession V, McKim Township, Sudbury, Ont. 


Palisade diabase collected Dec. 8, 1939, by P. M. Hurley from old quarry on 
Forest Ave., Staten Island. Specimen from 8 feet below original surface— 
fresh hard diabase, surface glaciated, and no weathering. 


Palisade diabase collected by P. M. Hurley, Dec. 8, 1939, from road cut on 
River Road at corner of 16th St., North Bergen, N. J. Specimen is a fine- 
grained trap, 20 feet below original surface. This point is the south end of 
the main Palisades. 

Palisade diabase collected by P. M. Hurley, Dec. 8, 1939, from large old quarry 
by side of River Road facing the Hudson, 1 mile north of 16Be, near Gutten- 
berg, N. J. Specimen of dense trap, 150 feet below original surface, possibly 
close to bottom of sill. 


Palisade diabase collected, Dec. 8, 1939, by P. M. Hurley from road cut on 
State highway 5 (N.J.) as it winds up the Palisade bluff above the 125th St. 
Ferry, Edgewater. Specimen from 70 feet below top of bluff and probably 
at least 30 feet from original surface. 


Palisade diabase collected Dec. 6, 1939, by P. M. Hurley from Kingston Trap 
Quarry 1% miles from Kingston, N. J., on road toward Rocky Hill, from a 
fresh blast, about 20 feet original surface; quarry operating. 


Diabase from Somerset Co., N. J., submitted for intercheck measurements by 
J. H. J. Poole, who reports (personal correspondence, 1938) Joly’s value as 
1.18 X 10°" g Ra/g; see 5A. 


Triassic diabase collected through the courtesy of the U. S. Geological Survey 
from Centerville, Va. Composed of calcic labradorite, about 40 per cent augite, 
a very little magnetite, and sparse chloritic and serpentine alteration products. 
Rock very fresh; grain size about 0.5 to 1.0 mm. 


Norite collected by E. M. Kindle from nickel-norite offset dike rock, Worthing- 
ton Mine, railway cut about 100 feet from Worthington Station, Ont. 


Norite collected by H. C. Cooke, 1939, from railway cut 15 feet below surface, 
lot 1, concession III, Wisner Township, Ont. 


Norite collected by E. A. Collins from hanging wall of pit at fault, Creighton 
Mine, Copper Cliff, Ont. Considered to be outcrop of 21B. 


Norite collected by E. A. Collins from Creighton Mine, 40 level, 14 X-cut, 
Copper Cliff, Ont. 


Norite, Lat. 45° 1’ N., Long. 109° 25’ W., collected July 1937 by A. C. Lane 
from the metamorphosed Stillwater complex of Quad Creek, Mont., 1 mile 
north of State line. The Stillwater complex is bordered by a Cambrian con- 
glomerate at an erosional contact and is cut by granite and by basic dikes. 
Age: Probably early pre-Cambrian. 


Trap collected by E. A. Collins from the Creighton Mine, 40 level, 14 X-cut, 
Copper Cliff, Ont. 
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Trap collected by E. A. Collins from the hanging wall of pit at fault, Creighton 
Mine, Copper Cliff, Ont. Considered to be outcrop of 23B: 


Trap rock, Lat. 44° 5644’ N., Long. 109° 29%4’ W., elev. 9800 feet, collected from 
a small working quarry by A. C. Lane, July 1937, 2 feet from the margin of a 
100-foot dike at Long Lake, Wyoming. The dikes of this series are intru- 
sive into the early pre-Cambrian complex (Stillwater series) and appear to be 
older than the middle Cambrian rocks of Beartooth Butte, but the exact 
relations are obliterated by overburden for a considerable distance. Age: 
possibly pre-Middle Cambrian, and post early pre-Cambrian. 


Trap rock, ditto 24B, except taken from the center of the 100-foot dike. 


oe rock, Lat. 44° 58’ N., Long. 109° 29’ W., elev. 10,600 feet, collected by 
A. C. Lane, July 1937, as part of Beartooth Range suite, from center of 25-foot 
dike, dip 45°. strike NNW. 


Monchiquite dike, Lat. 55° 45’ N., Long. 5° 3’ W., submitted by Arthur Holmes 
in 1936 from Riasg Buidhe on the Island of Colonsay, Scotland. Geologic re- 
lations nowhere provide evidence of age beyond the fact that it cuts forma- 
tions regarded as Torridonian. Age: Post-Devonian. 


Monchiquite dike, ditto 26B but from Kilchattan instead of Riasg Buidhe. 


Gabbro-diorite, 50-pound specimen collected by C. Goodman from General 
Crush Stone Quarry, Woburn, Mass. One of standard samples; see 3A. 


Sudbury gabbro collected by H. C. Cooke, 1939, from 10 feet below original 
surface in road cut, lot 3, concession III, Graham Township, Ont. 


Dunite, 50-pound specimen collected through the courtesy of the U.S. Geological 
Survey, from near Balsam Gap, N. Carolina. A part of the Webster dunite 
body. Composed dominantly of high magnesium olivine; a few per cent of 
chlorite, vermiculite, and chromite are present. Age is in question; some hold 
it to be pre-Cambrian; others that it is probably Ordovician. One of the 
standard samples; see 3A. 


Kimberlite, 50-pound specimen obtained through the courtesy of U. S. 
Geological Survey and S. J. Shand, from Kimberly, S. Africa. One of standard 
samples; see 3A. Fragmental rock composed of olivine, its alteration products 
consisting largely of serpentine. The only important primary mineral, other 
than olivine, is a very pale-colored type of phlogopite. Fragments of several 
invaded rocks are also present. Some of the olivine is completely altered, but 
other grains up to several mm. in diameter are largely fresh with only veinlets 
of serpentine. Age: lower Cretaceous. 


ecg collected by K. G. Bell (1939, p. 10) from Dalhousie Mtn. in Nova 

cotia. 

Andesite of unknown origin. 

Camptonite dike collected by R. J. Lougee from quarry in Burlington, Vt., 

visited by A-1 excursion of 1933 International Geological Congress. Cuts 

Cambrian quartzite; age possibly Triassic. 

Diorite collected by G. D. Louderback from Wutchumna Hill, north side, 

Lindsay-Strathmore canal, % mile NE of Lemon Grove, Calif., Lemon Grove 

quadrangle. Holocrystalline, granitic texture, medium grain, hornblende- 

hypersthene diorite. For petrographic and other radioactivity measurements 

see Finney (1934, p. 49). 

Pumiceous lava with sparse phenocrysts of olivine belonging to the flow from 

Halemaumau of March, 1921, submitted as a powdered specimen by C. S. 

Piggot for intercheck measurements; Piggot’s mean value, 0.94 X 10 g Ra/g, 

— geologic notes previously reported (Piggot, 1931, p. 3) using designation 
-10. 


Tower Falls basalt collected by A. C. Lane, 1937, in Yellowstone Park at about 
14 of Guide Book (Bighorn Basin Yellowstone Valley Tectonics Field Con- 
ference, Aug. 3-5, 1937), 134 feet vertically and 30 feet horizontally from the 
outcrop surface. Age: Recent. 
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Trap, Lat. 44° 5614’ N., Long. 109° 2914’ W., elev. 9800 feet, collected by A. C. 
Lane, as part of Beartooth Range suite, from altered glassy margin of Long 
Lake Dike (see also 24B:) by new road from Red Lodge to Yellowstone, near 
Long Lake, Montana, at about 35 of Guide Book (see 37B). Age: Possibly 
latest pre-Cambrian. 

Trap rock from Gogebic range sill, Mikado Mine, Wakefield, Mich. From 
drill hole No. C-18 at 105 to 110 feet below collar; hole vertical and showing 
about 20 feet of drift above ledge. Location: SE-NW, Sec. 18-47-45, Mich. 
Specimen submitted by A. C. Lane. 

Trap rock from Gogebic range sill, Mikado Mine, Wakefield, Mich. From 
drill hole No. G-48 at 855 to 860 feet below collar; hole vertical and showing 
about 20 feet of drift above ledge. This specimen is a more acidic phase than 
39B. Location: NE corner of SE-NW, Sec. 18-47-45, Mich. Specimen sub- 
mitted by A. C. Lane. 

Trap rock from Gogebic range, Mich., collected by L. M. Scofield, 1936. See 
Lane (1938, p. 66) for additional geologic notes and other radioactivity measure- 
ments. Lane’s designation: No. 31. 

Trap rock from Gogebic range, Mich., same suite as 41B. Lane’s designation: 
No. 27. 

Trap rock from Gogebic range, Mich., collected by L. M. Scofield, Jan. 16, 
1936. “Normal” from the 26th level of the Brotherton mine, cuts the upper 
phase of the Palms (Huronian). Lane’s designation: No. 37; see 41B and 42B. 


IcNeous Rocks INTERMEDIATE BETWEEN ACIDIC AND Basic 


Trachyte collected by A. D. Howard from SE base of Cresent Hill, Yellow- 
stone Park. 

Dacite porphyry collected by A. D. Howard from summit of Bunsen Park, 
Yellowstone Park. 

Gray tuff collected by W. D. Urry and W. L. Whitehead, Oct, 1935, 1 foot 
below glaciated surface and 2 feet above fossiliferous shale bed at fresh road 
cut, Newburyport Turnpike and Crowley Street, Rowley, Mass. Age: lower 
Devonian or upper Silurian, provisionally placed on dividing line. 
Granodiorite from La Motte, Quebec collection of R. A. Daly. Other 
measurements by Keevil (1938, p. 145). 

Tholeiite from Cleveland Dyke, Bolam County, Durham, England, submitted 
as a small hand specimen for intercheck measurements by Arthur Holmes 
(1937, p. 118); see 4B. 

Hornblende granodiorite, Lat. 48°14’N., Long. 77°54’ W., collected by 
W. C. Gussow in 1936 from a specimen freshly blasted from a road cut 10 
feet from the surface. The granodiorite is intrusive into a sedimentary 
complex of Keewatin or possibly earlier age and is thought to be older than 
the diabase dikes in the area and the Cobalt series to the southwest. Age: 
Probably pre-Cobalt and post-Keewatin. 

Quartz dolerite from Whin Sill, Scordale Beck, Westmoreland, England, sub- 
mitted as a small hand specimen for intercheck measurements by Arthur 
Holmes (1937, p. 118); see 4B. 


MeEtTAMoRPHIC Rocks 


Gneiss from Taroko-kyo, collection of Z. Hatuda; see 16A. 


Gneiss collected by T. F. W. Barth in Finland, designated as “augen Gneis” 
and classified as a migmatite. Lead has been quantitatively separated from 
a large sample of this rock by E. Gleditsch and coworkers and submitted to 
A. O. Nier for isotopic analysis in an attempt to determine the age by the 
lead method. 

Eclogite from Saastal, Switzerland, submitted as a small hand specimen for 
intercheck radioactivity measurements by J. H. J. Poole, who obtained 
0.42 « 10°" g Ra/g (personal communication, 1938); see 5A. 
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SEDIMENTARY Rocks 


Limestone submitted as a powdered specimen by F. Behounek, Director of 
State Radiological Institute of Czechoslovakia in 1936 for intercheck measure- 
ments. Behounek (personal communication, 1937) reports 24 « 10™ g Ra/g 
and designates as “Silurian Bituminous Limestone”. We observed the decrease 
in weight upon heating at 1000° C in air for 1 hour was 4.0 and 5.5 per cent 
for two portions of specimen. 

Shaly limestone core from 3780-00 feet, Moore No. 7 Schauers, Fitts Pool, 
Oklahoma. Age: Ordovician. 

Shaly limestone core from 3810-35 feet, Moore 1A Harper, Fitts Pool, Okla- 
homa. Age: Ordovician. 

Limestone, shale, quartz core from 4265-70 feet, Moore No. 7 Schauers, Fitts 
Pool, Oklahoma. Age: Ordovician. 

Shaly limestone, sandstone core from 4245-81 feet, Moore 1A Harper, Fitts 
Pool, Oklahoma. Age: Ordovician. 

Dolomite core from 4370-75 feet, Moore No. 7 Schauers, Fitts Pool, Oklahoma. 
Age: Ordovician. 

ae, Limestone, 50-pound specimen obtained through the courtesy of 
the U. S. Geological Survey from Carthage, Missouri. One of standard 
samples, see 3A. Composed exclusively of calcite with obvious fossil frag- 
ments of which crinoid remains are dominant. So completely cemented that 
no pore space is visible under the microscope; grains about 0.5 mm. 

Quartz sandstone core from 5108-16 feet, McKinney No. 5, Sun Field, Starr 
Co., Texas. Age: Oligocene. 

Quartz sandstone core from 4645-50 feet, Olivares No. 5, Sun Field, Starr Co., 
Texas. Age: Oligocene. 

Quartz sandstone core from 4657-8 feet, McKinney No. 6, Sun Field, Starr 
Co., Texas. Age: Oligocene. 

Sandstone core from 494 feet, Nowata Co., Oklahoma. Age: Pennsylvanian. 
Ditto 78: from 505 feet. 

Ditto 7S: from 512 feet. 


Shaly sandstone core from 2510-40 feet, Moore No. 7 Schauers, Fitts Pool, 
Oklahoma. Age: Pennsylvanian. 

Berea sandstone, 50-pound specimen obtained through the courtesy of the 
U. S. Geological Survey from Cleveland, Ohio. One of standard samples; see 
3A. Composed of subangular sand grains, a little chertlike quartz, and less 
calcite and feldspar. Age: Mississippian. 

Quartz sandstone core from 5889-94 feet, No. 1 Dailey, Navarro Crossing, 
Texas. Age: Cretaceous. 

Quartz standstone core from 5934-35 feet, No. 1 C. C. Hill, Navarro Crossing, 
Texas. Age: Cretaceous. 

Quartzite, 50-pound specimen obtained through the courtesy of the U. S. 
Geological Survey from Bull Run Mountain, Va. One of standard samples; 

see 3A. Composed of quartz in interlocking grains (thoroughly recrystallized), 
a little muscovite in parallel arrangement, and a few grains of magnetite. 

Clay shale core from 5096-7 feet, McKinney No. 5, Sun Field, Starr Co., 

Texas. Age: Oligocene. 

Clay shale core from 4649-53 feet, McKinney No. 6, Sun Field, Starr Co., 
Texas. Age: Oligocene. 

Clay shale core from 4636-7 feet, Olivares No. 5, Sun Field, Starr Co., Texas. 
Age: Oligocene. 

Shale core from 5851-67 feet, No. 1 Dailey, Navarro Crossing, Texas. Age: 
Cretaceous. 

Shale core from 5916-17 feet, No. 1 C. C. Hill, Navarro Crossing, Texas. Age: 
Cretaceous. 
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ABSTRACT 


The radioactive methods of geologic age determination involve the four basic 
requirements: (1) a known, systematic rate of disintegration of the radioactive 
elements; (2) accurate measurement and sampling; (3) absence of disintegration 
product as a primary constituent; (4) no addition or subtraction of the disintegra- 
tion product or its source during the history of the material. A critical review is 
presented of our knowledge of the fulfillment of these requirements by the various 
radioactive methods. Emphasis is placed on the helium method because of its 
possible wide applicability, because no adequate review has heretofore appeared, 
—— considerable confusion exists concerning the present status of this 
method. 

Following a discussion of the necessity for completely discarding all of Urry’s 
helium age measurements, a summary of the remaining data shows that (1) with 
only a few exceptions the helium age ratios lie in the proper geologic sequence when 
the results for acidic rocks and for basic rocks are considered separately, (2) the 
age ratios for acidic rocks are consistently lower than the age ratios for correspond- 
ing basic rocks, and (3) the helium age ratios, even for basic rocks, are substantially 
lower than corresponding lead ages on radioactive minerals. These general incon- 
sistencies stimulated the present reinvestigation of fundamental requirements under- 
lying the radioactive methods. The retentivity of rocks for helium has been found 
to be the major source of uncertainty in the present application of this method. A 
review of the lead method indicates that only comparatively few ages can be con- 
sidered as possible fixed points in the lead time scale. The uncertainty in the isotopic 
constitution of ordinary lead is seldom an important source of error, but the need 
for isotopic analysis in all lead age measurements is emphasized by the relative in- 
sensitivity of atomic weight determinations. Further work needs to be done on the 
effects of leaching and weathering of radioactive minerals. Brief reviews of the iso- 
topic ratio, rubidium-strontium, and possible potassium methods of age determination 
are included for the sake of completeness. 
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INTRODUCTION 


By the united application of the three main subdivisions of Historical 
Geology—Stratigraphy, Lithology, and Structural Geology—the relative 
sequence of many of the geological formations of the lithosphere have 
been established. 

The evidence from Paleontology, strictly speaking a branch of Stratig- 
raphy, is so distinctive that our extensive post-Proterozoic time scale is 
based on the characteristic forms of life which existed during each geologic 
period. However, the igneous rocks, which constitute more than 95 
per cent of the earth’s crust, contain no fossils and can be accurately 
placed in the geologic time scale only if there is a clear-cut structural 
relation to fossil-bearing sedimentary rocks. The vast expanse of pre- 
Cambrian time is represented by practically no fossils. 

The structural and paleontological methods, however, are inherently 
unable to give more than a qualitative indication of age—that is, of 
elapsed time in thousands or millions of years. A number of important 
results accrue from the establishment of a quantitative time scale. It 
provides a measure of the time required for the biological changes observed 
in the study of paleontology. A knowledge of the length of time involved 
in decomposition, disintegration, and deposition is of value in the study 
of Sedimentation. Similarly, the quantitative knowledge of ages serves 
to measure the rate of diastrophic processes, which is of interest to the 
structural geologist and petrologist. A more exact span can be assigned to 
certain of the periods and epochs than can be obtained from paleontological 
evidence, and a truer estimate of their relative lengths thus obtained. But 
the most valuable contribution to geology will come from the application 
of the quantitative method to the host of rocks that cannot otherwise be 
dated or even correlated. This need is especially acute for pre-Cam- 
brian rocks. 

The systematic disintegration of the geologically long-lived radioactive 
elements provides the most reliable method of measuring geologic time. 
Although we shall consider briefly methods based on the disintegration 
of some of the lighter elements, our discussion will mainly concern itself 
with age measurements based on the disintegration of the two parent 
elements, uranium and thorium, and their daughter elements, represented 
schematically as shown in Figure 1. 

It is seen that the element uranium is composed of two main isotopes,’ 
uranium I (UI) and actino-uranium (AcU), with atomic weights of 238 
and 235 respectively. The most accurate measurement of the relative 


lIsotopes of a given element have the same atomic number, hence possess the same chemical 
properties, but have different atomic weights (indicated by superscripts). 
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abundance of these two isotopes (Nier, 1939, p. 150) indicates 99.29 per 
cent UI and 0.71 per cent AcU. The transmutation of one atom of 
UI, through the uranium series, into a stable lead isotope, Pb°*, produces 
as a by-product 8 atoms of helium, He. Similarly, an atom of AcU 


THE RADIOACTIVE SERIES 


MOWo 


o 


Ficure 1—The radioactive series 


ends as Pb**? and produces 7 He, and an atom of Th ends as Pb®* and 
produces 6 He. 

This systematic generation of lead isotope end-products and of helium 
as a by-product are the properties of these parent radioactive elements 
which make possible their use as indicators of geologic time. The age 
equation can be expressed in the fundamental form: 


__ Amount of disintegration product (1) 


Age - 
r Rate of production of product 


in which the numerator represents the amount of radiogenic lead or 
helium formed and the denominator measures the rate, expressed in the 
proper units, at which the disintegration product is being generated in 
the material whose age is being measured. The Age refers to the time 
that has elapsed since final solidification or crystallization of the rock or 
mineral containing the radioactive elements. The application of this 
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INTRODUCTION 495 
equation to terrestrial materials involves the four basic requirements: 
(1) a known, systematic rate of disintegration of the radioactive elements; 
(2) accurate measurement and sampling; (3) absence of disintegration 
product as a primary constituent; (4) no addition or subtraction of the 
disintegration product or its source during the history of the material. 

In the light of present-day knowledge of nuclear physics, the first re- 
quirement is clearly well established. All efforts to alter the rate of 
radioactive decay have invariably led to negative results from which 
it is concluded that radioactive decay proceeds at the same rate inde- 
pendent of physical environment within the earth. Evidence from 
pleochroic haloes proves that the rates of decay of radioactive elements 
have been constant over geological periods of time. 

Evaluation of the extent to which the other three requirements are 
fulfilled by different geological materials has been the subject of con- 
siderable research and constitutes the major discussion of this paper. The 
problem divides itself into two parts: (1) the helium method as applied 
to ordinary rocks and minerals, and (2) the lead method as applied to 
radioactive minerals.? The reason for this general division of applicability 
can be understood from Figure 2. On a vertical logarithmic scale, repre- 
sentative total uranium (U)—. e., UI -+- AeU—and thorium (Th) con- 
tents for radioactive minerals, igneous rocks,® and a hypothetical rock, 
to be discussed later, are indicated in each case by the two solid lines on 
the left. The amounts of radiogenic lead (Pb*) and helium (He) that would 
be generated in these materials in 100 million years and the amounts of 
original ordinary lead (Pb) to be expected are indicated on the right. 
The lead content of such a radioactive mineral can quite readily be sep- 
arated quantitatively and weighed by the usual chemical techniques. 
However, the minute lead content found in ordinary rocks is considerably 
beyond the range of chemical methods, as indicated by the schematic 
arrangement of Figure 2. Similarly the uranium and thorium content of 
radioactive minerals can be determined gravimetrically, but the content 
of these elements in ordinary rocks can be determined only by physical 
methods using sensitive electrical detection apparatus. 

The helium content, corresponding to chemically separable quantities 
of radiogenic lead, would be several ce per gram.* Hence, it is hardly 
surprising that an appreciable fraction of the helium generated in even 
relatively young radioactive minerals has escaped, since the internal 


* Throughout this discussion the term ‘‘radioactive mineral’ is intended to designate terrestrial 
materials that contain uranium or thorium as important constituents, present to the extent of 0.01 
per cent or more, as distinct from rocks or minerals that contain only traces of these elements. 

%The U and Th values for the igneous rock lie midway between the averages for acidic and basic 
rocks given by Evans and Goodman (1941a, Table 10). 

‘In age measurements, helium volumes are referred to 0° C. and 760 mm. of Hg pressure. 
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pressure often reaches several atmospheres. It is this widely different 
nature of the two decay products, lead and helium, that divides the ap- 
plicability of the two methods of geological age determination. Thus, in 


1 
‘ 
; Pb = ordinary lead 
Pots radiogenic lead 
H 
-1 
Cc 
o 
be 
ro) ‘ 
= 
= tis 
H 
at = 4 ' 
A 
L 
a 5 ; ; 
tig 
7 
1 
' 
U Th PO'Pb He U Th PD Pb He U Th Pb'Pb He 
RADIOACTIVE IGNEOUS HYPOTHETICAL 
MINERAL ROCK ROCK 


Ficure 2.—-Schematic representation of radioactive content and 
decay products 


general, the lead method is limited to radioactive minerals, and the 
helium method to ordinary rocks and minerals. In the following sections, 
an attempt is made to evaluate the extent to which the fundamental re- 
quirements are fulfilled by the two methods. 
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By fortuitous circumstances a material might be found with an inter- 
mediate composition, such as that shown for the hypothetical rock in 
Fig. 2. In this case both the helium and lead methods might be ap- 
plicable to the same specimen. Thus far a direct comparison of this 
type has not been made. 


HELIUM METHOD 
GENERAL CONSIDERATIONS 


Since the discovery by Strutt (1908, p. 272) that the helium content of 
radioactive minerals could be correlated with their probable geological 
age, considerable work has been devoted to the application of this prin- 
ciple to the measurement of geologic time. It was shown quite early 
(Strutt, 1909, p. 166) that some helium escapes from radioactive minerals, 
even at room temperature and atmospheric pressure. This loss of helium 
lowers the numerator of our age equation and, hence, the resulting ratio 
is lower than the true age of the radioactive mineral. For this reason, 
age ratios on such minerals usually only indicate minimal values for the 
age.” As we have seen from Figure 2 and have mentioned in the intro- 
duction, the internal pressure of helium associated with these high con- 
centrations of uranium and thorium is so great that there is little wonder 
some helium escapes. Rather, the large quantities retained under these 
conditions over geologic time attest to the inherently high retention of 
helium in these minerals. The early investigators did not appreciate this 
viewpoint and did not have sufficiently sensitive apparatus to extend the 
helium method to igneous rocks. As a consequence, the helium method 
lay in disrepute for nearly 20 years following Strutt’s classical work. 

While radioactive minerals undoubtedly do not fulfill the fourth funda- 
mental requirement concerning helium, it does not necessarily follow that 
the same is true of igneous rocks, where the concentrations of helium, 
uranium, and thorium, as we have seen, are only about one-millionth of 
those found in radioactive minerals. With this in mind, Dubey and 
Holmes (1929, p. 794) made the first helium age determinations on 
igneous rocks. The results were so encouraging that a re-exploration of 
the long neglected helium method was begun. During the interval fol- 
lowing Strutt’s initial researches, the analytical techniques required for 
measuring the minute quantities of helium and radioactivity in rocks 
had been developed. Paneth and his coworkers (see Paneth and Urry, 


5 Throughout this paper and the one by Hurley and Goodman (1941), we are faced with the 
problem of differentiating between the age indicated by experimental measurements and the true 
numerical age of the geological material, which, we shall see, may be quite different. Accordingly, 
the term, age, represents the true age and the term, age ratio, the calculated quantity. If the four 
fundamental requirements are fulfilled for a given material, the age and the age ratio will, of 
course, be equal. 
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1930, p. 110, for references) had greatly improved the technique for the 
isolation and measurement of small quantities of helium. Joly, Poole, 
Piggot, Paneth, and Evans had developed methods of measuring the 
minute radioactivities involved. Thus, rapid strides in the application 
of the helium method were possible because of the availability of the 
necessary detection instruments. 

Under the guidance of the National Research Council Committee on 
the Determination of Geologic Time by Atomic Disintegration, A. C. 
Lane, Chairman, and supported by the American Philosophical Society 
and the Geological Society of America, an intensive program of helium 
age measurements was begun. Between 1932 and 1936 an elaborate 
helium time scale evolved largely as the result of measurements by 
W.D. Urry. (See Urry, 1936, p. 1219, for references.) Some investiga- 
tions of the basic geological assumptions underlying the helium method 
were made during this time. In general, however, the establishment of 
the validity of the method rested upon the accumulation of a large 
number of apparently concordant measurements. The most convincing 
evidence of reliability came from the apparent agreement with cor- 
responding ages based on lead ratios in radioactive minerals. Thus it 
appeared that a reliable time scale had been established which could be 
used as a reference for all geological age measurements. 

In the spring of 1936 the researches herein reported were begun. The 
direct-fusion furnace, which had proved so satisfactory in the liberation 
of radon from rocks (Evans, 1935, p. 102) appeared to offer similar ad- 
vantages in the determination of helium. Preliminary alpha-activity 
measurements on rocks (Evans, 1934 a and b, p. 29-42; Finney, 1934, p. 
1-98; Finney and Evans, 1935, p. 503-511; Raitt, 1935, p. 1-78) suggested 
the feasibility of this method of determining the total rate of helium 
production—~. e., the denominator of the age equation. The combination 
of these two techniques offered a direct, physical method of age de- 
terminations which would be independent of radioactive standards and 
much more rapid than the radon-thoron-helium method used by Urry 
(1936, p. 34-48). Since mass production of age results seemed in order, 
the development of assembly-line methods appeared appropriate. 

Initial measurements by the alpha-helium method, on rocks previously 
analyzed by Urry, gave age ratios lower by a factor of about two. Further 
study indicated that this difference was primarily due to the larger 
activity index ® obtained by the alpha-counting method as compared with 
that resulting from Urry’s radium and thorium measurements. At this 


®8In the alpha-helium method the rate of production of helium is observed as the number of 
alpha particles (helium nuclei) emitted per hour per milligram of rock, which we define as the 
“activity index.’’ This quantity is determined indirectly in the radon-thoron-helium method from 
the measured radon and thoron contents. 
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time the most likely explanation for the differences appeared to be the 
presence of a previously undiscovered radioactive series in terrestrial 
materials. This exciting possibility stimulated an investigation into the 
exact nature of the analytical differences. An extensive program of 
calibration and intercomparison of measurements on the same rock 
samples was carried out between the two methods. The results of this 
investigation have already been published by Evans, Goodman, Keevil, 
Lane, and Urry (1939, p. 931-946). It was established that there are 
no unknown strong alpha-emitting radioactive elements in ordinary ter- 
restrial materials.’ This conclusion follows from the good agreement 
between activity-index measurements made by direct alpha counting and 
those based on separate radon and thoron measurements of the same 
rock samples. This agreement appeared only after Urry’s radium ap- 
paratus had been recalibrated, the earlier calibrations having been 
found to be in error by a factor of over two. As previously mentioned 
(Evans and Goodman, 1941a, p. 470-471) part of this discrepancy can be 
explained by the use of unreliable radium standards, but the remaining 
difference has not yet been resolved. 

With a constant Th/U ratio, the use of a constant radium correction 
factor would lower evenly all the earlier helium ages. However, as is 
evident from the Th/U values given by Evans and Goodman (1941a, 
Tables 6-9), this ratio varies considerably, and a correction, applied to 
the radium measurements alone, would cause an uneven decrease in 
Urry’s helium time scale. Another difficulty is that when the earlier 
radium measurements are recomputed, using the corrected calibration 
constant, the ages based on these values are also higher than more recent 
measurements (by both the alpha-helium and radon-thoron-helium meth- 
ods) on the same hand samples (Evans et al., 1939, Table XI, p. 948). 

The apparent agreement between ages based on the lead method as 
applied to radioactive minerals and Urry’s helium ages on fine-grained 
igneous rocks was cited (Urry, 1936, p. 1230; Holmes, 1937, p. 179) as 
the most cogent evidence for the validity of both methods. With the 
downward revision of helium age ratios necessitated by the newer studies, 
this apparent agreement, which we now know was merely fortuitous, 
disappeared. In its stead we were left with only a limited number of 
reliable values constituting a series of helium age ratios that appeared 
to be in a sequence consistent with geologic evidence, but in magnitude 
at variance with the corresponding lead ages. At this juncture two pos- 
sibilities appeared likely—either the lead method involved some previ- 
ously undisclosed error, or one or more of the fundamental requirements 


7Evans and Raitt (1935, p. 174) have previously shown that there are no strong unknown 
gamma-ray emitters in terrestrial materials. 
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of the helium method had not been fulfilled by the rock samples measured. 
An inquiry into the lead method, most of which is included in a later 
section of this paper, showed that while considerable uncertainty exists 
in many of the lead age ratios, some are unusually well authenticated. 
Attention was then turned to a critical evaluation of the various uncer- 
tainties underlying the helium method, a summary of which was pre- 
sented before the Geological Society (Mead, Evans, and Goodman, 1939, 
p. 1921), and the details of which are given in subsequent sections of 
the present paper. In brief, it was concluded that with the proper care 
the second requirement of accurate measurement and sampling can be 
fulfilled. The possible presence of primary helium can be avoided by a 
not too restrictive selection of material, thus satisfying the third require- 
ment. The fourth requirement, particularly as it pertains to the loss 
of helium, remained as the major source of uncertainty in the application 
of the helium method to rocks. Accordingly, a systematic program of 
research was undertaken, in collaboration with Dr. P. M. Hurley, which 
was specifically designed to test the fulfillment of the fourth requirement 
by various rocks and minerals. The details of this work are presented 
by Hurley and Goodman (1941, p. 549-556) but, for the sake of com- 
pleteness in the present consideration, are also summarized here as follows: 


The retentivity of rocks for helium has been found to be extremely variable and 
is the major source of uncertainty in the application of this method. Researches 
on separate mineral components indicate that, with only a few exceptions, helium 
ages obtained on rocks as a whole, even fine-grained, unaltered mafic rocks, can 
represent only minimal values. Certain minerals, notably the feldspars, appear to 
lose a large fraction of their radiogenic helium. On the other hand, certain of the 
mafic minerals, pyroxene and magnetite in particular, have a much higher retentivity 
and show helium ages which approach the values to be expected from the lead time 
scale. Thus, all helium ages made on rocks as a whole represent mean values de- 
pendent upon the proportion and helium retentivity of the mineral components. It 
is anticipated that further research on selected rock minerals will yield helium ages 
suitable for the establishment of a new helium time scale in substantial agreement 
with the present lead time scale. 


DETERMINATION OF HELIUM 


General statement —We have seen that the fundamental age equation 
involves knowledge of the helium content, in the numerator, and the rate 
of production of helium, in the denominator. In this section the tech- 
niques employed for the measurement of the helium will be briefly 
reviewed. 

Most rocks contain between 10~° and 10~* ec of helium per gram. The 
determination of these minute quantities involves: (1) liberation, (2) 
purification, and (3) measurement. There are a number of ways in which 
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these three steps can be performed, but the techniques summarized below 
appear to offer the optimum combination of reliability and convenience. 


Liberation of gases—The development of the direct-fusion furnace and 
its application to radioactivity measurements have been described previ- 
ously by Evans (1935, p. 102). In principle the same type of furnace 
is used in the release of helium from terrestrial materials. However, 
because of the significant helium content of the atmosphere (see footnote 
16), a vacuum-tight design is necessary. For this reason, the earlier 
methods of sealing have been superseded by a double-gasketed face plate 
and special inset electrode seals. These and other features can be seen 
in the photographs (Pl. 1). In Figure 1 of Plate 1, the furnace is opened, 
and the upper graphite radiation screen removed to show the mounting of 
the graphite crucible, the glass vacuum leads, and the rubber hose water 
connections for cooling. The swinging bed allows convenient removal 
of the sublimate that lines the furnace after each fusion of rock. In 
Figure 2 of Plate 1, the furnace is closed and ready for action; the mercury 
bubbler trap is seen just above the peep-hole in the top of the furnace. 

A large volume of gas, of which the helium constitutes only a minute 
fraction, is released during the fusion of a rock sample. This gas acts as 
a flushing medium in quantitatively sweeping all helium out of the fur- 
nace and through the mercury bubbler trap which prevents back diffusion. 

Experimental evidence for the quantitative release of helium from rocks 
by direct fusion (at 1500° to 2000° C.) is afforded by the absence of helium 
in repeat measurements on the rock melt at elevated temperatures (2000° 
to 2500° C.). Studies of the fractional liberation as a function of 
temperature indicate that 99 per cent or more of the helium is released 
from most rocks within 5 minutes at temperatures only slightly above 
1500° C. It is routine practice to heat the sample for at least 5 minutes at 
a temperature of 2000° C. or above. Comparative measurements on quar- 
tered portions of a number of granulated rock samples show good agree- 
ment between the carbonate-fusion and direct-fusion techniques (Evans, 
et al., 1939, Table VIII, p. 941). 


Purification—Claude (1914, p. 861) was the first to conduct extensive 
research on the purification of the noble gases. He showed that helium 
and neon are absorbed to a negligible amount on charcoal at liquid air 
temperatures when compared with the other inert gases and hydrogen. 
Paneth and Urry (1930, p. 110) applied this method to the separation of 
pure helium from the other gases released when meteorites were dissolved 
in acid or fused with alkali carbonates. 

The most troublesome contaminant is hydrogen, since it is only par- 
tially absorbed on the charcoal. The method previously used (Paneth 
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and Urry, 1930, p. 118; Urry, 1933, p. 340) for removing the last traces 
of hydrogen was the introduction of pure, electrolytic oxygen followed by 
circulation of the gases over heated palladium (approximately 300° C.). 
This procedure catalyzes the formation of water vapor which is subse- 
quently absorbed on charcoal together with any excess oxygen. 

A bulb containing liquid oxygen has been found to be a somewhat more 
convenient and equally pure source of oxygen. The general arrangement 
of the apparatus is represented by the left-hand side of Figure 3. The 
liquid oxygen is poured into a flash equipped with two outlet tubes, one 
leading to a bubbler through which the excess gas can escape, the other 
leading into the main body of the system through a pyrex tube con- 
taining stick cupric oxide. The cupric oxide, when heated, serves the 
same purpose as the palladium but is not susceptible to poisoning and 
hence does not decrease in efficiency with use. 

The charcoal tubes (Fig. 3, C1-C5) dipped in liquid air are used in 
successive order for absorption of all the gases except helium, which re- 
mains distributed throughout the vacuum-tight, pyrex system. 


Measurement.—After final purification the volume of helium is meas- 
ured in the modified McLeod gauge equipped with a ground capillary 
(Rosenberg, 1938, p. 258). This type of capillary has been found to 
have a more uniform correction (for capillarity) and hence yields more 
reproducible results than the smooth-walled variety used by other 
workers. Ten or more measurements of the length of the gas column in 
the capillary are made at various pressures, observed by the differential 
height of the mercury in the side-arm tube. These observations, after 
the capillary correction is applied, usually show a maximum variation of 
less than 5 per cent from the mean value of the pressure-volume product. 
By the use of the perfect gas law, involving the observed temperature 
and calibration constant of the apparatus,® the volume of helium con- 
tained in the known weight of rock sample can be readily calculated. 

An accurate estimate of the experimental uncertainty in helium 
measurements on rocks is difficult to make, largely because of the some- 
what inhomogeneous distribution of helium in rocks. Repeat measure- 
ments on separate portions of the same hand specimen usually agree well 
within a probable error of 10 per cent of the mean, and on quartered 
portions of the same sample within 5 per cent of the mean. Some rocks, 
however, show variations in helium distribution of several times this 
amount. The relatively narrow range of values observed for the calibra- 
tion constant (Evans, et al., 1939, Table VII, p. 941) indicates that the 


* The calibration constant is defined as the ratio of the effective total volume of the apparatus 
to the volume of the McLeod gauge. The capillary is calibrated by weighing measured lengths 
of mercury at a known temperature. The diameter of the bore generally used is about 0.25 mm. 
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purification procedure is quantitative to within 2 per cent for moderate 
volumes of helium, around 10-* ec. Thus for most measurements the 
sampling error is predominant. For volumes of 10~° cc of helium, the esti- 
mated probable error amounts to about 30 per cent. 


DETERMINATION OF RADIOACTIVITY 


Methods.—Two methods can be used to determine the rate of produc- 
tion of helium in terrestrial materials. The older method involves sep- 
arate measurements of the uranium and thorium and the use of their 
decay constants to calculate the corresponding rate at which helium is 
being produced. In actual practice the radon and thoron contents are 
determined, and the proportion of the parent elements calculated, assum- 
ing radioactive equilibrium. Since modern detection instruments make it 
readily possible to count individual alpha rays, a direct, physical method 
of determining the rate at which helium is being produced in terrestrial 
materials has been developed. This newer method possesses the distinct 
advantages of being free from chemical manipulation and having only 
a second-order dependence upon standard radium or thorium solutions 
and upon the decay constants of uranium, thorium, and actino-uranium, 
which enter only in a small correction term for very old rocks. Both 
these methods have been used for the age analyses reported in this paper 
and in the paper by Hurley and Goodman (1941). Consequently, a brief 
review of the two methods is included. 


Alpha-count method.—Since the only known natural source of helium 
is from the alpha particles emitted by terrestrial radioactive sources,’ the 
true rate of production of helium is given directly by the total alpha 
emission within a rock. However, it has already been shown by Finney 
and Evans (1935, p. 503) that such measurements are complicated by 
the small range of alphas in these materials and the variation of range 
for alphas from the different radioactive elements. For example, the 
alpha rays of uranium I have a mean range of only about 14 micra in 
rock, while the normal alphas from thorium C’ penetrate only about 46 
micra. Except for the rare, long-range alphas from the C’ bodies and 
the very short range ones from samarium,’® the alphas from the 20 other 
emitters lie between these limits. For sources thick to these alpha rays, 
the conversion of the observed count to the true emission depends directly 
upon an accurate knowledge of the stopping power of the rock, the pro- 
portion and ranges of the various alpha emitters, and an experimental 


® The recently discovered (Alvarez and Cornog, 1939, p. 613) stable isotope, He®, in ordinary helium 
has a relative abundance of less than 10-7, and its origin is not yet known. 

10In practically all rocks, the helium contributed by samarium is negligible compared with that 
from the other radioactive elements. 
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factor involving the minimum detectable range of alpha rays in the 
counting apparatus (Finney and Evans, 1935, p. 504). While each of 
these quantities is either known or can be measured within a probable 
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Ficure 4—Alpha-count conversion curves 


error of about 10 per cent, the combined application results in consid- 
erable uncertainty of measurement for thick sources, which usually off- 
sets the advantages to be derived from the alpha-counting method. Al- 
though rock sources sufficiently thin to cause only a fractional absorption 
of even the shortest range alphas give a somewhat lower counting rate 
than can be obtained from thick sources, the conversion of the observed 
count to the true emission only involves the above factors in terms 
of second-order importance. Consequently, considerable improvement in 
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accuracy is obtained by the use of thin rock sources, and the present 
discussion is limited exclusively to such deposits. 

Following the general equations derived by Finney and Evans (1935, 
p. 503-511) which relate the observed counting rate to the true rate of 
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Ficure 5—Graphical relationship between F and Th/UI 


emission of alpha particles within the source, the graphical relationships 
shown in Figure 4 can be established." The ratio of the true to the observed 
rate of alpha emission, (8Ny + 7N,4 + 6Nr)/(ny + ng + np), is obtained 
for a given thickness of source t plus absorber a, by interpolation between 
the curves for F = 0 to F=1.0. The quantity F represents that fraction 
of the total alpha rays produced in the source due to the uranium and 
actino-uranium series. Consequently, values of F from 0 to 1.0 cover the 
range of all possible values and are related to the thorium-uranium ratio, 
Th/UI, as shown in Figure 5. At a Th/UI ratio’? of 4.13—that is, 
F = 0.5—the rate of production of helium by the thorium series is equal 
to the rate of production of helium by the uranium plus actino-uranium 
series. 

Referring again to Figure 4, it is seen that for a thickness of rock source 
plus absorber of t+ 2a=1.0 air-cm., which corresponds to slightly 
more than 5 micra, a maximum uncertainty of less than 2 per cent is 


1 The details of this derivation and a general description of the alpha-helium technique are 
being published in physical journals (Evans and Goodman, 1941b). 

12 Average values of 5.0 and 4.0 were obtained for the Th/U ratio (essentially the same as 
the Th/UI ratio) of acidic and basic rocks respectively by Evans and Goodman (194la, p. 474). 
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introduced by assuming F = 0.5. Hence, for most rocks, accurate knowl- 
edge of the proportion of Th and U is unnecessary in determining the 
denominator of the age equation by direct alpha counting. 

To obtain the approximate age ratio, to, in billions of years (10°) from 
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Ficure 6.—Old age correction curves 


the observed helium content, He, and activity index,!* I, the general ex- 
pression, Equation (1) is rewritten as follows: 

Helium in 10-*ce/g of rock 

Activity Index in a/hr-mg 


Age =t, =0.0308 


In all rocks the present-day rate of production of helium is less than 
in the past due to the wearing out of the very slowly decaying parent ele- 
ments. This effect is not significant for rocks younger than about 500 
million years but introduces an increasingly important correction as the 
age increases above this value. Thanks to the reliable researches of 
Kovarik and Adams and of Nier, we now have an accurate knowledge 
of the decay constants of uranium, thorium, and actino-uranium, which 
allows this old-age correction to be made with considerable certainty. 
The true age ratio, t, is obtained from the approximate age ratio, t,, by 
reference to Figure 6. It is seen that, for ages < 1 billion years (t= 1), 


In the alpha-helium method, I = (8Nu + 7Na + 6Nr) and is obtained by multiplying 
the observed rate of alpha emission, (mv + ms. + nr), by the conversion factor obtained from 
Fig. 5. In the radon-thoron-helium method, I = (1.04 Ra + 0.0886 Th) where Ra = radium 
content in 10-!2g/g of rock and Th =thorium content in 10-®g/g of rock. 
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a maximum error of 5 per cent results from the assumption that F = 0.5. 


For greater ages, an accurate value of F can always be obtained from a 


radium determination, because: F = bl 


Radon-thoron method—The direct-fusion technique of measuring the 
radium, and hence the uranium, content of rocks has been discussed in 
detail by Evans and Goodman (1941la). The determination of thorium 
from the thoron content has likewise been mentioned insofar as it pro- 
vides additional evidence for the validity of the thorium values obtained 
from alpha-count measurements. Instead of repeating this material, the 
present discussion supplements these previous remarks with details that 
particularly pertain to the age determinations reported in a later section 
and in the paper by Hurley and Goodman (1941). 

We have seen in the preceding section that a negligible experimental 
error in the determination of the rate of helium production is introduced 
by the lack of knowledge of F for alpha counts on rock deposits of about 
5 micra in thickness. Since the area is about 40 cm?, such sources only 
represent a mass of 60 mg of rock. Except for quite homogeneous rocks, a 
single measurement on this small quantity of sample cannot be expected 
to give a very representative average value for the radioactive content. 
Consequently, it is customary to grind between 1 and 2 grams of sample 
from which several deposits are made for alpha counts. The average 
value is then used for the age calculation. This procedure gives quite 
satisfactory results unless the activity index of the sample is unusually 
low—i. e., less than 0.4. Below this activity, the radon-thoron method 
is definitely superior to the alpha-count method. For this reason both 
methods are in active operation at the Massachusetts Institute of Tech- 
nology laboratories and supplement each other in covering the entire 
range of activities found in terrestrial materials. 

The thoron method being used at present is similar to that described 
by Urry (1936, p. 35) but embodies a number of changes in technique 
that appear to improve the reproducibility of the results. Since Urry 
(1933, p. 335) employed a sodium carbonate fusion for releasing helium 
from the rock sample, he used the same flux for radium and thorium as 
well. For inhomogeneous rocks, this procedure offers a distinct advantage 
in using the same portion of sample for all three measurements. A 
pyrex, bell-jar furnace suitable for carbonate-fusion helium measure- 
ments is shown in Figure 3. Instead of the thimble-shaped, nickel 
crucible used by Urry, a hemispherical nickel crucible, 5.0 em in diameter, 
is used, because it allows more rapid solution of the melt following 
fluxing. The use of sintered-glass filter plates with suction greatly 
accelerates the otherwise very laborious chemical manipulations and 
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reduces the uncertainty concerning the time of separation of the thorium 
X from its parent, radiothorium. It has also been found that nitric acid 
solutions of thorium X have a higher and more reproducible emanating 
power than corresponding hydrochloric acid solutions. For this reason, 
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Ficure 7.—Diagram of thorium apparatus 


the procedure of dissolving the final residue, after removal of silica, in 
50 ce of dilute nitric acid has been followed. 

This solution is then placed in a 75 cc flask and the thoron measurement 
begun using the apparatus represented in Figure 7. Radioactively inert 
nitrogen ** is bubbled at a constant rate, as indicated by the tapered glass 
flow meter, through the solution and continuously sweeps the thoron into 
the ionization chamber. The gases are then exhausted through an upper 
outlet tube, but a definite fraction of the thoron atoms disintegrate in the 
chamber (Evans, 1935, p. 109-112). The ionization created by each 
alpha ray from thoron and its decay products is collected by the electro- 
static field between the center wire and the shell of the chamber. The 
resulting pulses are amplified by the electrometer tube (FP-54) and fed 
to a sensitive galvanometer. The position of the galvanometer is followed 
photographically. Since the 6” by 37” bromide paper used for recording 
moves at a known rate, the alpha particle emission rate in alphas per 


4 Ordinary tanked nitrogen, stores for 30 days or longer to allow any radon present to com- 
pletely decay, has been found quite satisfactory. 
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hour is obtained directly from a visual count of the developed record. 
When the emission rate is plotted against the time, a histogram similar to 
Figure 8 is obtained. It is seen that the rate increases abruptly above 
the normal background activity upon introduction, at the 20 hour mark, 
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Ficure 8.—Thorium calibration and background runs 


of a solution containing thorium X, in this case 15 ce of our temporary 
thorium standard. (See Evans and Goodman, 1941a, p. 469-470). The 
average activity for the 8 hours during which streaming continued is 
found to be 190 alphas per hour above a background of 22 alphas per 
hour. Other calibrations over the range of activities observed for most 
rocks give the curve shown in Figure 9. 

It should be noted from Figure 8 that the background following the 
rock run does not fall immediately to the normal value. Instead there is 
clearly a gradual decrease over a period of 50 to 60 hours after removal 
of the thorium X solution. This activity is attributable to the active 
deposit of thorium B* collected on the central electrode of the ionization 
chamber. This element decays with a half-period of 10.6 hours, which is 
much longer than that of the other daughter products of thorium X. 
Hence, the rate of decay of the residual deposit activity is dictated by 


1% The genetic relation between thorium B and the other members of the thorium series can 
be recalled by reference to Figure 1. 
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thorium B. In practice it is quite inefficient to wait 50 to 60 hours be- 
tween runs, yet accurate knowledge of the background is essential for 
reliable measurements. Accordingly, the procedure is to make extended 
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Ficure 9—Calibration curve of thorium apparatus 


determinations of the normal background every 5 or 6 runs and observe 
the decaying background for 16 hours after each run. Knowing the decay 
rate, the true average background for each rock measurement can be 
computed from the combined observations. This schedule makes possible 
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one thorium measurement or calibration run every 24 hours during the 
week and an extended normal background measurement over the week- 
end. The most important desideratum is to make the calibration, rock 
measurement, and background runs all under as close to identical con- 
ditions as possible. Toward this aim the above schedule is followed for 
both calibration and rock measurements, and the same procedure fol- 
lowed for backgrounds with distilled water substituted for the thorium X 
solution. Under ordinary conditions, the moisture picked up from the 
solution in the flask would condense to some extent in the ionization 
chamber and cause electrical leakage across the amber insulation of the 
center electrode. A few turns of resistance wire placed around the 
chamber and dissipating 5 to 10 watts raises the temperature in the 
chamber sufficiently to prevent this condensation. 

After completing the thorium measurement, the same solution can be 
used for a radium determination using the procedure for liquid solutions 
described by Evans (1935, p. 103). In brief, the procedure is: (1) trans- 
fer the solution to another flask, (2) seal this flask to one of the reflux 
condensers shown in Figure 3, (3) de-emanate the solution by heating, 
while bubbling through radioactively inert nitrogen, (4) shut off for a 
known period of time (usually several days) during which a known frac- 
tion of the total radon is produced by the radium in the solution, and (5) 
transfer this radon to the ionization chamber and measure the rate of 
alpha emission by means of the apparatus shown in Figure 3. Calibration 
with standard radium solutions (see Evans and Goodman, 1941a, Figure 
1) following the same procedure, gives the conversion to equivalent 
radium content. The main difficulty with the procedure is the uncertainty 
due to radium contamination from the chemicals and glassware used in 
the fluxing, purification, and dilution of the rock sample. Even with very 
elaborate precautions, it has been found very difficult to avoid picking up 
between 0.2 <X 10-"* and 0.3 X 10-1” grams of radium in the over-all 
procedure, as indicated by blank measurements. This correction, which 
is not always reproducible, becomes very significant for rock samples of 
low radium content and may even exceed the radium content of 5 grams 
of sample. Accordingly, it is preferable to limit this method to rocks with 
higher than average radium contents and make direct-fusion measure- 
ments on separate portions of sample for the less active rocks. The direct- 
fusion furnace technique avoids all these sources of contamination in 
radium measurements. 

PRIMARY HELIUM 

Rayleigh (1938, p. 410) and Shepherd (1938, p. 311) have established 
the presence of the noble gases, notably argon, in rocks. Although there 
is a possibility that some argon may be produced from the disintegration 
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of potassium, as shown in a later section, the major proportion of these 
inert gases is very likely of primary origin. Since magmas are known to 
contain uranium and thorium, they must also contain some helium due 
to radioactive disintegration. If other inert gases in rocks are derived 
from the magma, it is only reasonable to suppose that at least a part of 
the helium content of some igneous rocks must be of a similar origin. If 
primary helium constitutes a significant fraction of the observed helium 
content of a rock, the third requirement is not fulfilled and the age cal- 
culated from Equation (1) will be in error on the high side. 

At first sight, it would seem that helium measurements on recent effusive 
rocks would provide a direct means of ascertaining the amount of primary 
helium. However, the texture of these rocks vitiates their usefulness. 
As the result of their mode of origin, effusives are generally either glassy 
or porous. Since it is well established (Paneth and Peters, 1928, p. 253; 
Urry, 1932, p. 3887; 1933, p. 3242; Rayleigh, 1936, p. 350) that many 
glasses are selectively permeable to helium, and since the atmosphere 
contains helium '*, glassy, effusive rocks are likely to contain sufficient 
atmospheric helium to mask any contribution from their magmatic 
sources. The helium content of a number of recent Vesuvian lavas, as 
observed by Strutt (1907, p. 436), is probably of atmospheric derivation. 
Other more porous effusives have been found by the authors to contain 
neon in the proportion of about three parts to one part of helium by 
volume, which suggests that air is either absorbed or entrapped in the 
rock. 

The information obtained from studies of recent trap rocks is of con- 
siderably more value in studying the question of primary helium. The 
absence of air from the environment of the rock at the time of its in- 
trusion eliminates the possibility of absorption of atmospheric helium. 
However, helium measurements alone are not sufficient to evaluate the 
presence of primary helium, for it is not possible to obtain specimens of 
these intrusive rocks that are geologically very recent. In addition, the 
helium content of rocks, even of the same age, varies over a wide range 
owing to the considerable variation in the radioactive content, and hence 
of the rate of production of helium, as shown by Evans and Goodman 
(1941a). For this reason, it is necessary to measure the radioactivity of 
the rock as well as its helium content. These data combine, according 
to Equation (1) to give the age of the rock. Thus we see that helium 
ages themselves provide the most direct means of studying the question 
of primary helium. But the term, age, can justifiably be applied only 
when the four fundamental requirements are rigorously fulfilled. Hence, 


In a recent compilation, Paneth (1937, p. 436) gives the following percentages for the rare 
gas constituents of air: Ne: 1.8 X 10-8, He: 5.3 x 10-*, Kr: 1 & 10-*, He: 5 X 10-5 and Xe: 8 X 10-4. 
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in this circuitous method of evaluation, in which helium ages are used to 
test one of the requirements on which helium ages are based, one must 
assiduously avoid the pitfall of reasoning in a circle. 

First, we must assume that we have some qualitative notion of geo- 
logical time. Suppose we regard 50 million years as a relatively short 
span in geological time. (There is considerable independent evidence 
supporting this generalization.) Secondly, we must assume that by 
selection of samples we are able to evaluate one uncertainty at a time. 
This is necessary, since if primary helium is present and helium leakage 
has occurred subsequent to solidification of the rock, the two effects are 
mutually compensating, and it is not possible to evaluate either on a 
given specimen. Suppose for the moment we assume that fresh, fine- 
grained, basic igneous rocks do not allow helium to diffuse either into or 
out of themselves. Then age ratios on recently formed intrusives of this 
type should give some indication of the presence or absence of primary 
helium. 

Measurements have been made on a number of such Tertiary rocks, 
and the results are presented in Table 1. It is seen that the observed 
helium contents are quite small. Even if all this helium were of primary 
origin, it would be of negligible amount compared to the radiogenic 
helium that would be produced in these rocks within a few hundred 
million years by the radioactive contents observed. Hence, corresponding 
amounts of primary helium in much older rocks would not significantly 
affect the measured helium age, except in cases of exceptionally low 
radioactivity. Since we have every reason to believe that the specimens 
in Table 1 are typical, they furnish fairly definite evidence that for trap 
rocks at least the amount of primary helium, if any, is usually negligible 
in helium age measurements. It remains for future research to evaluate 
this question for the more plutonic rocks.'* 


MIGRATION OF HELIUM AND OF RADIOACTIVE ELEMENTS 


General considerations—We now consider the fourth and final re- 
quirement of the helium method—namely: has there been any addition 
or subtraction of helium or of the radioactive elements during the history 
of the material whose age we are attempting to measure? Since the 
preceding three requirements can be fulfilled simultaneously for certain 
selected rock samples, the fate of the helium method appears to hinge on 
the fulfillment of this fourth requirement. 

Following the establishment of the reliability of the analytical tech- 


17 There is some evidence that primary helium occurs in certain of the Sudbury rocks, which may 
explain to some extent the anomalous helium age ratios observed (Hurley and Goodman, 1941, 
p. 546). However, the data are too limited to draw any definite conclusions, and in addition these 
rocks are not typical of plutonics in general. 
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niques, it became increasingly apparent that only after a thorough study 
of the fourth requirement had been made would a true evaluation of the 
helium method be possible. Accordingly, in the Fall of 1939, a systematic 
investigation of the problem of helium retention in rocks and minerals 
was initiated in collaboration with Dr. P. M. Hurley of the Massachusetts 
Institute of Technology Geology Department. Although the problem of 
the migration of the radioactive elements in rocks is far from settled, the 
question of possible helium migration seemed to be the more uncertain of 
the two aspects of the fourth requirement. For this reason, study of the 
latter has been undertaken first. It is in the sampling of the proper 
materials by the geologist for analysis by the physicist that the co-opera- 
tive aspects of such research becomes of paramount importance. Our 
present knowledge of the fourth requirement is the sum of the following 
inventory of information and the results of the recent work reported in 
the paper by Hurley and Goodman (1941). 


Helium migration—A number of laboratory studies have been made 
of the diffusion of helium in rocks and minerals. Urry (1933, p. 3242) 
and Rayleigh (1936, p. 350) each measured the permeability to helium of 
one rock specimen and a number of minerals. Because of the minute 
quantities of transmitted helium, relatively large pressures were required. 
The general nature of the results indicates a very low permeability. 
However, extrapolation to geological conditions involves so many uncer- 
tainties that these observations only afford qualitative evidence for the 
retention of helium over geologic expanses of time. 

The fractional loss of helium at elevated temperatures has been studied 
for a very limited number of meteorites and rocks and for a more 
extensive, but for our purpose less important, number of radioactive 
minerals. Paneth and Urry (1930, p. 140) have demonstrated that the 
loss of helium, from three iron meteorites maintained at 1000°C for 
several hours, was equivalent to the helium contained in only a thin 
surface layer. Such meteorites show the greatest retention of helium of 
all mineral matter studied, with the possible exception of certain of the 
native metals. 

The escape of helium as a function of grain size and time has been 
studied for a few rock samples. One 8-gram sample of each of two fine- 
grained specimens, a basalt and a felsite, were found by Urry (1933, p. 
335) to lose less than 0.05 and 0.03 per cent respectively of their helium 
per 100 years, based on measurements made after being sealed in air-free 
oxygen at atmospheric pressure for 3 months. Urry also found no dif- 
ference in the helium content of a Keweenawan basalt, whether taken in 
one piece, several small pieces, or finely powdered, after standing for 24 
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hours between grinding and fluxing. In contrast to these results, Holmes 
and Paneth (1936, p. 385) observed a 60 per cent loss of helium from 
finely powdered kimberlite kept in a glass tube with a cork stopper at 
atmospheric pressure for 4 years. Too much emphasis should not be 
placed on these measurements since the helium leakage from powders is 
only of real significance in evaluating the analytical techniques. It is to 
be expected that considerable helium might escape from a rock powder 
prepared from a rock which in its native state would retain its helium 
completely. Truly significant measurements of this type should be made 
on specimens of at least hand size. 

The most elaborate measurements thus far made on the helium reten- 
tion of rocks and rock-forming minerals are those reported by Holmes 
and Paneth (1936, p. 385). Although these observations are quite 
complete from an analytical standpoint, it is to be regretted that this 
excellent work had to be done on specimens which are so atypical of 
average igneous materials. The excessive hydrothermal alteration of 
the various rocks composing the Kimberley diamond pipes indicates 
clearly that such rocks cannot be expected to yield anything approaching 
reliable helium ages. The relatively high helium contents of certain of 
the basic appinites and eclogites are significant only in that they indicate 
retention of a large fraction of their helium under the most adverse 
geological conditions. The phlogopite glimmerite and basic granulite are 
so highly altered that helium loss from these specimens can only represent 
extreme conditions compared with those to be expected for typical igneous 
rocks. 

For the sake of completeness, some mention should be made of the 
loss of helium from radioactive minerals, which in general contain about a 
millionfold more helium than ordinary igneous rocks (Fig. 2). In his 
early work on the age of such minerals, Strutt showed that there is an 
appreciable loss of helium even at room temperature. Kano and 
Yamaguchi (1928, p. 244) have studied the yield of helium from samars- 
kite and ishikawaite, with and without the addition of flux, when held 
under vacuum for 2 hours at temperatures from 200 to 1000°C. These 
investigations also indicate that more than 90 per cent of the helium is 
lost from ishikawaite during the first 30 minutes when held at 700°C for 
8 hours. Time intervals shorter than 30 minutes were not used. Other 
researches (Chlopin, Herling, and Joffe, 1934, p. 28) have shown that the 
release of helium from radioactive minerals is greatly accelerated by the 
presence of hydrogen in the surrounding gases. Since hydrogen very 
seldom occurs in nature, its influence would be rare, although this fact 
may have an important bearing on the genesis of helium-bearing natural 
gases. 


on 
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A piece of work '® particularly interesting since it was made only one 
year after the discovery of helium on the earth was carried out by 
Ramsay and Travers (1897, p. 325). They found that fergusonite, when 
heated to a temperature not exceeding 500°C, becomes suddenly in- 
candescent and evolves much of its helium, while its density decreases. 
The investigators thought at that time that the effect was due to “break- 
down of a helium compound” and they called fergusonite an “endothermic 
mineral.” X-ray study in recent years has shown that fergusonite 
normally has not a crystalline atomic structure but on heating to 400°C 
it becomes crystalline, tetragonal. This suggests the possibility that a 
temperature condition, which causes a state of dynamic disorder either 
among part or among all of the constituents of a crystal structure, might 
release part of the helium held in the structure, even if the disorder is 
only momentary. 

With modern knowledge of crystal structure, it should be possible to 
make a theoretical estimate of the probabilities of diffusion of helium 
through minerals. The most important assumption involved in such 
considerations is that of assigning the distribution of the radioactive 
elements which generate the helium. Practically no definite information 
is available concerning this interesting problem. In the absence of 
experimental evidence it is generally assumed that the uranium and 
thorium atoms are randomly distributed throughout the crystal structure 
of the rock mineral. On this basis Keevil (1940, p. 311) has recently 
shown theoretically that a surprisingly large retention of helium should 
exist for igneous rocks which have not been subjected to hydrothermic or 
dynamic metamorphism. Feldspar, in particular, is considered in detail, 
and Keevil arrives at the conclusion that this mineral, when fresh, should 
not lose an appreciable fraction of its helium during a length of time 
comparable to the age of the earth itself. Hurley and Goodman (1941) 
present experimental evidence which indicates quite the contrary in the 
case of feldspar. Accordingly it must be concluded that this theoretical 
treatment does not adequately include all factors involved. 

As is so frequently the case, laboratory facilities cannot duplicate 
geological conditions, and in order to evaluate the fourth requirement 
reliably we must base most of our conclusions on field evidence. By the 
judicious selection of certain critical specimens, it is possible from age 
measurements to make important deductions concerning the retention of 
helium in rocks. If there is any appreciable loss or migration of helium 
either in the rock body as a whole or in the individual mineral consti- 
tuents of the rock, different values for the helium age ratio should be 
obtained, since the helium would no longer be in immediate association 


18 The authors are indebted to Dr. P. M. Hurley for calling their attention to this early work. 
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with the corresponding amount of the parent radioactive elements that 
produced it. In the course of his extensive age measurements of the 
Keweenawan, Urry (Lane and Urry, 1935, p. 1101) selected a few rocks 
for their bearing on this question of helium retention. These measure- 
ments appeared to afford convincing evidence of the lack of migration of 
helium. However, as has been mentioned previously, the radium 
measurements obtained in these studies were too low by a factor of more 
than two. Because of the variations in the thorium-uranium ratio of 
these specimens the downward revision of age ratios results in widely 
discordant values, with the result that these critical measurements on 
helium retention must be discarded entirely. 

The variation of the helium age ratio, laterally or with depth in the 
same igneous body, would suggest the migration of helium. Two of the 
early measurements by Dubey *® on acidic felsite of the Mt. Pawagarh 
plug indicated the same age ratio for a specimen from the summit and 
one from 100 feet below. Holmes has discounted the absolute value of 
these age ratios because of the cryptocrystalline nature of the rocks. In 
recent years Keevil (1938, p. 123) found the same age ratio over a depth 
of 600 feet in Quincy granite, determinations of ages having been made at 
five points over this range. Further age measurements, on this well- 
known rock, from quarries separated by a tenth of a mile and one mile, 
also gave the same age ratio, indicating a lack of differential migration 
of helium. Keevil also obtained the same age 2 feet from the margin and 
in the center of a Keweenawan trap dike. However, these latter results 
are not thoroughly convincing evidence of retention of helium since the 
radioactivities were essentially the same, and hence no gradient of 
helium existed. 

Some of the first age measurements made by Dubey (1930, p. 807) on 
the Gwalior series gave progressively lower age ratios with an increase in 
grain size for three rocks that were considered to represent different 
phases of the same basic intrusive. These discrepancies in age were 
attributed to the loss of helium in coarser-grained materials. It is inter- 
esting to note that in the original publication of these results the textures 
of the rocks were indicated merely as fine-grained or medium-grained. 
Subsequent recapitulation has given three gradations of grain size, in- 
creasing in the order of decreasing values of the reported age measure- 
ments. A number of recent determinations by Goodman and Keevil, 
included with other results in a later section, fail to confirm this variation 
of age ratio with texture for basic rocks. 

A number of metamorphic and sedimentary rock types are known to be 
unsuitable for helium age determinations because of the escape of helium. 


19 Personal communication to Arthur Holmes; see National Research Council, 1931, p. 413. 
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The major portion of the rocks reported by Holmes and Paneth (1936, 
p. 385) fall in the class of metamorphic rocks comprising hydrothermally 
altered igneous rocks. The age ratios obtained indicate extreme con- 
ditions of helium loss. 


TaBLE 2.—Age measurements on metamorphosed rocks 


Sample P Helium content Activity index in* Age ratio in 
No. Geologic age in 10-5 ce/g @ /hremg 106 years 

38B Pre-Cambrian......... *2.7 +0.2 **0.97+40.05 85+5 
4M Pre-Cambrian......... 2.2 +0.2 *0.82+0.05 8448 
5M Pre-Cambrian......... 2.4 +0.2 *2.0 +0.1 3743 
6M Pre-Cambrian......... *0.33+0.04 *0.18+0.01 5649 
7M Cambrian (7)......... *1.3 +0.1 0.88+0.04 46+5 

31A eee 3.4 +0.3 *2.8 +0.1 37+3 

18 Pre-Cambrian*f....... 0.43 0.84 16 

20 Pre-Cambrian*f....... 3.7 1.17 98 


+ Activity index measurements by alpha count except for samples 18 and 20, where radon-thoron 


method was used. 
+ Determinations by Holmes and Paneth (1936, p. 408). 
* Mean of two determinations. 

** Mean of four determinations. 


Two of Holmes and Paneth’s results on pre-Cambrian rocks together 
with measurements by the authors on six geologically old rocks are given 
in Table 2. These age ratios are much lower than those observed for 
unaltered rocks of comparable geologic age.*® The most likely explanation 
is that metamorphic changes cause at least a partial loss of helium from 
rocks. While the introduction of uranium and thorium would afford an 
alternative explanation, this process is much less likely, as indicated in 


the next section. 


Migration of radioactive elements—The migration of radioactive ele- 
ments commonly occurs as the result of movement of solutions (Hender- 
son, 1939, p. 2388) in which these elements are either suspended or 
dissolved. Some migration occurs in the gaseous phase. Diffusion in the 
solid state appears to be exceedingly rare. With only a few exceptions, 
the important uranium and thorium minerals are formed by high-tem- 
perature replacement of earlier pegmatitic minerals. The concentrated 
solutions of these parent elements, that effect the replacement, represent 
the uranium and thorium content of large volumes of magma. These 


2 There is some possibility that the period of metamorphism may be estimated by studies of 
this type, but the uncertainties involved in such an application of the helium method are even 
more numerous than those for unaltered rocks. 
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solutions escape from the molten rock as the result of the filter-pressing 
action of the crystallized fraction of the cooling rock mass. 

The migration of the parent radioactive elements into rocks, as the 
result of replacement, or out of rocks, as the result of leaching, will 
respectively decrease or increase the age ratio. The only value we can 
observe for the denominator of the general age equation is the present- 
day rate of production of helium, which, as we have seen, is proportional 
to the radioactive content. If any significant amount of replacement 
occurred in a common rock, there would be a number of indications seen 
upon petrographic examination. Rocks of this type should not be studied 
by the helium method until sufficient knowledge has been obtained re- 
garding the importance of this effect. Similarly, fresh, unweathered 
specimens should always be chosen, since leaching may cause erroneous 
results. These considerations emphasize the need for adequate petro- 
graphic study prior to helium age measurement. This policy has only 
been followed in the more recent age research. 

A study of the distribution of the radioactive elements over a regional 
scale provides a check on the possible migration of these elements since 
solidification of the rock mass. It is difficult to conceive of a process 
that would result in uniform addition or subtraction of both uranium and 
thorium throughout a large mass of rock. The results of a study of the 
Triassic intrusives of New Jersey, Pennsylvania and Virginia are given 
by Hurley and Goodman (1941). Good agreement is observed between 
age ratios for selected minerals separated from specimens of the same 
rock mass over a lateral range of 45 miles. These results indicate that 
in this rock body, at least, no significant migration of the radioactive 
elements has taken place since its solidification. 


AGE MEASUREMENTS 


General statement—We have seen that the reliability of the helium 
method can be evaluated to some extent on the basis of independent 
physical and geological knowledge, but in the final analysis this method, 
or any other geological tool, must be judged by its success in actual ap- 
plication. The road of research is seldom smooth, and this is particularly 
true in the development of radioactive methods of age determination. 
Considerable work has been done during the past 3 decades, but much of 
it has been discarded as newer developments disclosed earlier misconcep- 
tions and discrepancies. Our present position can best be appreciated by 
a review of the past progress. 


Early investigations —Following the suggestion by Rutherford in 1905 
that the accumulation of helium in radioactive minerals might provide a 
means of measuring geological time, Strutt made an extensive survey of 
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the helium, uranium, and thorium content of a number of the more 
radioactive terrestrial materials. He obtained helium age ratios for 
phosphatic deposits and hematite (1908, p. 272), other iron ores (1909, p. 
96), zircons (1909, p. 298), and sphenes and other minerals (1910, p. 194). 
From these results a broad correlation between the geological age of a 
mineral and the helium age ratio was established. However, Strutt (1909, 
p. 166) was the first to point out that, owing to leakage of helium from all 
these minerals, the age ratios observed could only represent minimal 
estimates of the true age. These helium age ratios were found to be only 
30 to 70 per cent of the values based on the accumulation of lead in 
corresponding radioactive minerals. For this reason, it has been necessary 
to discard most of the excellent helium age researches of Strutt. 

Waters (1909, p. 677; 1910, p. 903) made preliminary helium age 
studies of some common rock minerals. Unfortunately, the electrical 
detection methods used were too insensitive to measure other than the 
more radioactive mineral inclusions. Because these materials showed 
evidence of helium leakage, the work was discontinued. 

Following the development of methods of measuring minute quantities 
of helium, Paneth and his coworkers determined helium age ratios for a 
large number of meteorites (Paneth and Urry, 1930, p. 127). The magni- 
tude of the age ratios obtained, together with some laboratory tests on 
the materials, indicates that iron meteorites have a remarkably high 
retentivity for helium, even at elevated temperatures. Similar measure- 
ments were attempted on native metals, but their extremely low radio- 
activities prevented application of the method. Subsequently, Paneth 
has collaborated with Holmes (1936, p. 385) on an investigation of the 
rocks composing the Kimberley diamond pipes. These measurements 
are mainly of value in emphasizing the type of material not to be used 
for helium age studies. 

The first real stride in the application of the helium method to rocks 
can be said to have begun with the work of Dubey and his collaborators. 
These results, summarized in Table 3, illustrate a number of interesting 
points. The felsite (1) shows a much lower age ratio than the Cleveland 
Dike (2), although the two rocks are of approximately the same geologic 
age. The age ratios of (2 - 5) are in the proper sequence, but the actual 
values observed are somewhat lower than the corresponding lead ages. 
Specimens (5 - 7) were selected to test the effect of texture on the re- 
tention of helium as well as to estimate the geologic age of the Gwalior 
series of India. In the original publication Dubey (1930, p. 807) listed 
(5) as a fine-grained basalt and both (6) and (7) as medium-grained 
dolerite from a Morar-group sill. Holmes, (National Research Council, 
1931, Table LX XVI, p. 415), presumably on the basis of a more exact 
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description, subsequently indicates (6) as medium-grained and (7) as “of 
coarser grain than 6.” With this revision in the texture, it is seen that 
the age ratios for (6 - 8) decrease with increasing grain size. In the 
light of these observations, it seemed desirable at that time to concentrate 
on close-grained, crystalline basaltic or other basic rocks.”* 


TaBLeE 3.—Age measurements of igneous rocks 


Helium con- Activity Age 
Rock specimen Geologic age tent in index Rn-Th ratio in 
10-5 ce/g method 108 years 
2. Cleveland Dyke...... Oligocene or Miocene... 0.11 0.12 28 
3. Deccan basalt....... Early Eocene.......... 0.12 0.10 37 
Late Carboniferous. .... 0.36 0.057 196 
5. Gwalior basalt... .... Late pre-Cambrian..... 0.55 0.034 500 
6. Morar dolerite....... Late pre-Cambrian..... 0.42 0.033 389 
7. Morar dolerite....... Late pre-Cambrian..... 0.31 0.047 203 


1-7 Geologic notes, references and further discussion, National Research Council, 1931, p. 411-415. 
* Mean of two determinations. 


Following this work, researches were initiated in 1931 by Dr. W. D. 
Urry in collaboration with Professor A. C. Lane on an extensive suite of 
Keweenawan rocks and minerals, 


“with a view to co-ordinating age results with certain geological and mineralogical 
aspects pertaining to the specimens, in addition to the important question of the age 
and duration of the Keweenawan” (Lane and Urry, 1935, p. 1119). 


The magnitude and scope of this work were indeed admirable. As 
subsequent events have proved, however, it was premature. Age ratios 
were obtained for more than 30 rock specimens and minerals. With only 
a few exceptions, excellent agreement with geologic conditions was ob- 
served. The magnitude of the age ratios also appeared to be in accord 
with equivalent lead ages. On the strength of this encouragement age 
ratios were determined on 39 selected rock specimens ranging in geologic 
age from Cambrian to Pliocene. Urry (1936, p. 1226) reports that of 
these results 


“only five are incompatible with the correct geologic sequence as determined from 
field relations. All five specimens should be rejected as unsuitable on strict adherence 
to the standard conditions set for selecting material.” 


The helium time scale represented by the summation of these measure- 
ments was also found to agree quite well with the somewhat less extensive 


21 Holmes (National Research Council, 1931, p. 414) states as follows: ‘A systematic programme 
of research has already been started by Dr. Dubey in India and by the writer in collaboration 
with Prof. Dr. F. Paneth, and Dr. R. W. Lawson who have arranged to make helium determinations 
and to investigate quantitatively the escape of helium under various conditions.” Unfortunately ; 
these investigators have found it possible to complete only partially this admirable program. 
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lead time scale. This mutual compatibility of the two radioactive 
methods of age determination appeared as the final conclusive proof of 
the reliability of both methods when properly applied. While there were 
a number of gaps left in the structure and much work remained to be 
done in the pre-Cambrian, in general, the geologic time scale seemed a 
definitely established framework on which future research would build. 

Concurrent with these developments (1929-1936), Evans was engaged 
in the development and application of sensitive detection instruments for 
the measurement of terrestrial radioactivities. (See Evans, 1935, p. 99-112, 
for key references.) The logical extension of this work was in the 
direction of age measurements. The physical methods used offered a 
number of analytical advantages that should make analyses more accurate 
and more routine and, hence, more generally available for geologic 
problems. 


Intercalibration measurements—In February 1936 collaboration be- 
tween the authors began on helium age studies. The first step taken was 
to test the reliability of the alpha-counting method of determining 
radioactivity and the direct-fusion method of measuring helium, using 
portions of hand samples of rock previously analyzed by Urry. These 
preliminary determinations, on samples 19B, 34B, and 31, gave age ratios 
substantially lower than those obtained by Urry. The differences ap- 
peared to be primarily due to a lack of »greement on the radium contents. 
However, it also appeared possible that the discrepancies in radium 
results might be due to inhomogeneities in the samples. Accordingly, a 
set of Keweenawan rocks, 39B to 43B inclusive, were measured by the 
alpha-helium method. These age ratios were in even more drastic dis- 
agreement with Lane and Urry’s (1935, p. 1108) extensive list of ages for 
the Keweenawan. Thus it seemed advisable to make a number of direct 
comparisons between the alpha-helium and the radon-thoron helium 
methods. Portions of the same granulated rock samples and of the 
standard solutions used for calibration purposes were measured independ- 
ently in the two laboratories.2?. The results of this work already have 
been summarized in the preceding general discussion. For the purposes 
of the present review, we need only repeat that, because of apparently 
unsystematic errors in the radioactivity measurements, all of Urry’s 
helium age ratios had to be discarded. 


Summary of age data—The helium age data remaining after this 
drastic, but necessary, re-evaluation are summarized in Tables 4 and 5. 


2A detailed description of this intercalibration and comparison has been given by the par- 
ticipants, Evans, Goodman, Keevil, Lane, and Urry (1939, p. 931-946). 
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The values given include the earlier results of Dubey and Holmes from 
Table 3 and one measurement by Holmes and Paneth. The remaining 
age ratios include those obtained in connection with the interchecking 
work mentioned above, other concurrent measurements by Keevil (work- 
ing in Urry’s laboratory), and subsequent determinations by the authors. 
The helium age ratios fall into two groups, one for basic rocks (Table 4) 
and one for acidic rocks (Table 5). These have been placed in their ap- 
propriate positions on the left hand side of the time scale shown in Figure 
10. The field evidence, on which the geological age is based, is in some 
cases only very approximate, and in these instances the position assigned 
is the midpoint of the range within which the age can be defined. The 
lead ages shown on the right side are also divided into two groups, one for 
lead ratios and one for lead isotope ratios. These two methods yield 
results that are in general agreement, although there are notable excep- 
tions, for example the Permian pitchblende from Bohemia and the 
Cambrian Swedish Kolm. These lead age methods and the results ob- 
tained by their application are discussed in some detail in the next section. 
For our present purpose, the important generalizations to be drawn from 
Figure 10 are: (1) With only a few exceptions the helium age ratios lie 
in the proper geologic sequence when the results for acidic rocks and for 
basic rocks are considered separately, (2) the age ratios for acidic rocks 
are consistently lower than the age ratios for corresponding basic rocks, 
and (3) the helium age ratios, even for the basic rocks, are substantially 
lower than corresponding lead ages on radioactive minerals. These three 
general inconsistencies in helium age ratios indicated quite clearly that 
there was some fundamental failure in the helium method. A suggestion 
of the possible difficulty had come from Keevil’s (1938, p. 406-416) in- 
vestigation of the Quincy granite. Keevil observed no lateral or vertical 
variation of age ratio in this rock body, but he did find marked differences 
in the age ratios for the three main mineral components of this rock. A 
theoretical treatment of the probability of helium leakage from minerals 
(Keevil, 1940, p. 311) with particular consideration of feldspar, which 
was the Quincy rock mineral that gave the lowest age ratio, led Keevil ** 
to conclude that the helium retentivity of practically all igneous rock 
minerals is sufficient to ensure the retention of essentially all the helium 
produced in them during geologic time. On the strength of this conviction, 
Keevil (1938, p. 406) attributed the observed differences in age ratios to 
leakage of helium resulting from alteration of the feldspar (chiefly 
microperthite). He states (p. 412): “if the feldspar crystals were per- 
fect, no loss of helium could occur, and even with a normal amount of 


% This investigation was made in 1937 although its publication was delayed until 1940. 
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crystal imperfection, the permeability would be low.” However, Keevil 
did suggest that if further work on fresh minerals substantiated these 
observations (p. 416), “absolute helium age determinations will be ob- 
tained only from certain minerals and not from the rock as a whole.” 


Recent developments.—In view of these numerous uncertainties in the 
applicability of the helium method, it appeared essential that a really 
fundamental investigation be made of the basic requirements underlying 
this method. The preceding review of our present knowledge concerning 
the first three requirements constitutes a part of this work. The progress 
made to date in collaboration with Dr. P. M. Hurley on the evaluation 
of the fourth requirement is summarized in the paper by Hurley and 
Goodman. Definite differences in helium retentivity have been found 
for various common rock minerals. These differences offer a very logical 
explanation for the variations in age ratios observed for basic and acidic 
rocks, owing to the higher feldspar content of the latter. It can also be 
stated with considerable surety that all helium age ratios on rocks con- 
taining feldspar as an essential mineral are too low owing to the leakage 
of helium from this mineral. Since all the acidic rocks in Table 5 and 
most of the basic rocks in Table 4 contain appreciable proportions of 
feldspar, the large majority of these age ratios must be considered only 
as minimal values. Hence, these results are of little use in the establish- 
ment of a helium time scale or for correlation purposes. However, there 
is a possibility that the helium retentivity of even the more permeable 
rock minerals can be established as a fairly specific quantity. If this 
can be done, a knowledge of the mineral composition of some of the rocks 
previously studied may enable estimates to be made of the true helium 
age. Until this investigation has been made, it appears necessary to in- 
clude in the helium time scale only those measurements made on selected 
rock minerals. This restriction limits the acceptable age ratios to those 
reported by Hurley and Goodman. Fortunately, the suitable minerals 
appear to be rather widely distributed in igneous rocks. Magnetite, in 
particular, offers considerable promise. Age-ratio measurements made on 
magnetites from widely varying horizons appear to be in a sequence 
consistent with geologic knowledge and in quite good agreement with 
corresponding lead ages. By the continuance of this work, it is anticipated 
that a new helium time scale can be established on a sound scientific 
basis. Although complicated by the necessity for greater selection in the 
choice of samples, the new helium method should provide a valuable aid 
to geologists in the solution of the numerous important problems awaiting 
a reliable dating method for rocks. 
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LEAD METHOD 
APPLICABILITY 


By reference to Figure 2, the fundamental difference between the 
helium method and the lead method of measuring geologic time has been 
discussed. It is evident that with present analytical techniques the lead 
method must be restricted, almost exclusively, to the study of radioactive 
minerals. In the present section an attempt is made to evaluate the 
extent to which the second, third, and fourth basic requirements for reli- 
ability are fulfilled by radioactive minerals. These considerations seem 
particularly appropriate at a time when lead ages are being used as 
indubitable points of reference on the geologic time scale. 


ANALYTICAL TECHNIQUES AND ACCURACY 


The quantitative determination of Pb, U, and Th in radioactive 
minerals can usually be made with an uncertainty of only a few per cent. 
The precision of chemical separations and gravimetric measurements, of 
course, depends on the concentrations of these elements in the mineral. 
There is also some dependence on the chemical nature of the mineral; 
certain of the more siliceous minerals are particularly difficult. 

The presence of lead as an impurity in the reagents may be a significant 
source of error unless due precautions are taken, as Wells (see Lane, 1936, 
p. 66) has point out. Some of the earlier work on Pb/U and Pb/Th 
ratios undoubtedly suffers from this defect. However, a more important 
cause of inaccuracy in these measurements was the failure to determine 
the amount of ordinary lead which, as indicated in Figure 2, often occurs 
in appreciable quantities in radioactive minerals. It is necessary to be 
able to identify the two types of lead, ordinary, Pb, and radiogenic, Pb*, 
in order that the numerator of Equation (1) accurately represent only 
lead generated by radioactive disintegration. 

Three methods have been used for distinguishing between ordinary and 
radiogenic lead. These three methods—(1) atomic weight determina- 
tion **, (2) determination of the abundance of the lead isotopes mass 
spectrographically, and (3) a similar determination from the hyperfine 
structure of the lead spectrum—all depend upon the difference in atomic 
weight of the individual lead isotopes, Pb®°*, Pb?°*, Pb?°? and Pb?°*. Be- 
cause of its convenience, accuracy, and need for only small quantities of 
sample, together with the fact that in some cases it gives additional in- 
formation not otherwise obtainable, the mass spectographic method has 
largely superseded the other methods. Aston (1927, p. 224; 1929, p. 313) 

24Since the atomic weight of ordinary lead, 207.21 (chemical scale), is substantially different 
from that of uranium lead, approximately 206.03, and thorium lead, approximately 207.98, the 


proportions of the different leads in the mineral lead can often be estimated from an atomic 
weight determination. 
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made the first isotopic abundance measurements of ordinary and radio- 
genic leads. Nier (1938, p. 1571; 1939, p. 153) has recently greatly im- 
proved the accuracy of this technique and has made a systematic study of 
quite a number of ordinary and radiogenic leads whose atomic weights 
had previously been determined by chemical methods. From these re- 
sults it appears that the present analytical techniques are adequately 
accurate for the identification of lead in age measurements of the large 
majority of radioactive minerals. 
PRIMARY LEAD 

In the preceding discussion a distinction has been made between 
ordinary lead, Pb, and radiogenic lead, Pb*. Until recently it was 
assumed that these two leads represented the only types found in terres- 
trial material. This assumption was based on the observation (National 
Research Council, 1931, p. 218; Baxter, et al., 1935, p. 467; 1937, p. 702) 
that the atomic weight of ordinary lead from a wide variety of sources 
was apparently constant. Nier (1938, p. 1574) has recently observed 
significant variations in the relative abundance of the isotopes of a num- 
ber of ordinary leads. These variations, however, produce little or no 
effect on the chemical atomic weight of the lead, which is a very insensi- 
tive indication of isotopic composition where four isotopes are involved. 
These variations are much larger than can be explained by possible 
isotopic separation taking place in nature. Instead it seems likely that 
what was previously considered to be ordinary lead is itself a mixture of 
primal lead and varying amounts of radiogenic lead. In order to dis- 
tinguish between these several forms for purposes of discussion, the fol- 
lowing definitions are included: 


Primal lead.*®—From spectroscopic analysis of solar radiation, it appears that the 
atoms composing the sun are identical with those of our earth. It also appears prob- 
able that the stable atoms have remained essentially the same throughout the life- 
time of our planet. Accordingly, some of the lead found in rocks and minerals today 
must constitute part of the original inheritance of the earth. The term “primal lead” 
is used to represent this form of lead. Since the isotope Pb™ is not produced by 
any known radioactive disintegration, it is considered as indicative of the proportion 
of primal lead. However, the isotopic constitution of primal lead is shrouded in 
mystery; the best approximation for the abundance of its isotopes comes from a 
study of very old lead minerals. Thus far the lead from the Great Bear Lake galena 
of pre-Cambrian age (estimated to be 1300 million years) has shown the highest 
proportion of Pb™ and is tentatively selected as representing primal lead. The lead 
from the Joplin, Mo., galena of Late Carboniferous age (estimated to be 230 million 
years) has the lowest proportion of Pb®* of all the ordinary leads thus far studied. 
Nier’s (1938, p. 1574) values for the relative abundance of the isotopes of these leads 


are as follows: Atomic Wt. (Chem. Scale) 
Galena 204 206 207 208 Physical Chemical 
Method Method 
Great Bear 1.48 23.59 22.64 52.29 207.218 207.206 
Joplin, Mo. 1.26 27.31 20.00 51.43 207.178 207.22 


% The general adoption of this term was suggested by Professor A. C. Lane during the Conference 
on Applied Nuclear Physics, held at the Massachusetts Institute of Technology in October 1940. 
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A comparison of the atomic weights shows the relative insensitivity of the chemical 
method. 


Ordinary lead—By the term “ordinary lead” is meant the lead of commerce de- 
rived from such minerals as galena, cerussite, and wulfenite. Other terms that have 
been used for this purpose are: “ore lead”, “primary lead”, and “common lead”. 
Ordinary lead is composed predominantly of primal lead together with varying small 
amounts of radiogenic lead. It is the occurrence of ordinary lead in radioactive 
minerals that complicates the lead method of geologic age determination. 


Radiogenic lead—The lead resulting from the radioactive disintegration of uranium 
and thorium is designated as “radiogenic lead”. 


Primary lead—The lead present in the mineral at the time of its formation. This 
initial time is the zero point in the age of the mineral—that is, it represents the time 
t=0. While there is some uncertainty in the designation of this time, depending 
on the mode of formation of the mineral, this uncertainty is probably quite small 
in comparison to the geologic age of the mineral. 

From the determination of the abundance of Pb? in the lead separated 
from a radioactive mineral, it is thus possible to correct for the amount of 
primary lead. Nier’s results show that while this correction can generally 
be made with satisfactory accuracy it necessarily involves some uncer- 
tainty owing to the lack of specific knowledge of the isotopic constitution 
of this lead. Because of the extreme difference in the geochemistry of Pb 
and that of U and Th, the amount of ordinary lead in radioactive minerals, 
except those of recent origin, is generally much less than the radiogenic 
lead produced in them during their lifetimes. For this reason, the un- 
certainty due to variations in the isotopic abundance of ordinary lead is 
seldom significant, although it may dictate an appreciable correction for 
very young minerals. These considerations emphasize the need for 
isotopic analysis in all lead age measurements. 


ALTERATION OF RADIOACTIVE MINERALS 


The fourth basic requirement as applied to the lead method involves 
not only the question of migration of the disintegration product, lead, but 
also the possible variation in the concentrations of the parent elements 
as the result of hydrothermal alteration. Radioactive minerals are par- 
ticularly susceptible to leaching, replacement, oxidation, and reduction. 
Microscopic examination of radioactive minerals shows that they are in- 
variably intensely shattered. Ellsworth (1932, p. 88) has observed that 
the degree of shattering in radioactive minerals increases with the uranium 
and thorium content, and that the rock matrix surrounding radioactive 
minerals is likely to be penetrated by fractures which radiate in all 
directions from the uranium- or thorium-bearing minerals. These effects 
are hardly surprising in view of the continual bombardment of the min- 
eral and its surroundings by the energetic radioactive radiations and the 
atomic and molecular rearrangements which are necessitated by the finite 
recoil of the active atoms with each disintegration. The resulting heating 
effect and internal pressure of helium supplement in the development of 


; 
ae 


LEAD METHOD 533 


stresses which are sufficient to cause fractures throughout the mineral 
within a geologically short space of time. 

These fractures afford openings through which solutions can percolate 
in the rocks and minerals. The resulting alteration may profoundly 
affect the lead ratio of the mineral. There has been considerable specula- 
tion concerning the preferential leaching of lead from radioactive minerals 
(Holmes, 1937, p. 172), especially thorium minerals, but at the present 
time our knowledge concerning this possibility is exceedingly meager. 
Primary uranium minerals appear to be formed under reducing condi- 
tions,2* and it is generally assumed that all, or nearly all, the uranium 
was originally in the lowest stage of oxidation represented by UO.. As 
the result of subsequent oxidation, whether by autoxidation as advocated 
by Ellsworth (1932, p. 89) or the oxidizing effects of the meteoric solu- 
tions, many uranium minerals show definite layers containing high per- 
centages of UO;.?”7_ In the form of this higher valence state, uranium is 
probably more susceptible to leaching action than lead. Consequently, 
altered uranium minerals may show misleadingly high lead age ratios. 

Another possibility, that would result in quite the opposite effect, is the 
loss of intermediate members of the radioactive series from the mineral. 
For example, the radioactive gases radon, Rn, thoron, Tn, and actinon, 
An, might escape to some extent from the mineral by diffusion. If this 
should happen, the end-product lead would be deposited elsewhere. Two 
factors oppose this loss: (1) These gases are relatively short-lived °°, 
hence are only present in minute amounts at any instant and could not 
travel far without decaying, and (2) these gases are so heavy that their 
rates of diffusion are extremely small. However, it seemed worth while 
to test this possibility in at least one instance. 

The Swedish kolm is the only material of high radioactive content that 
is precisely dated by geologic evidence; distinctive fossils prove it to be 
of latest Cambrian age.*® Ironically enough, this unique material is 
particularly resistant to analysis. The lead age ratios observed average 
about 400 million years, whereas the lead isotope ratio for the kolm shows 
a value of 770 million years (Nier, 1939, p. 159). Since escape of radon 
would be in the direction to account for this difference, measurements 
were made to evaluate this factor. A 750-gram hand sample of kolm was 


% Carbon-hydrogen and hydrocarbon minerals are found associated with uraninites and other 
radioactive minerals in pegmatites in Canada, Sweden and Russia (Ellsworth, 1932, p. 88). 

7A method of estimating geologic age by the proportion of oxygen to that of uranium and 
lead, assuming that all the uranium was originally in the form of UOs, was suggested by Ellsworth 
and has been the subject of investigation by Lane and Khlopin. This method, however, requires 
study of the crystal structure and density of the mineral, as well as accurate determinations of 
oxygen and the ratio of UOz to UOs. 

%Trn = 3.82 days, Trn = 54.5 secs., and Tan = 3.92 secs. 

® For further details, including various analyses and geologic notes, see National Research Council, 
1931, p. 299. 
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found to lose 0.25 X 10-° curies of radon per 24 hours when sealed in 
nitrogen at atmospheric pressure. Gamma-ray measurements on the 
sample indicated a radium content of 0.8 & 10~—*g, which corresponds to 
about 0.3 per cent uranium. These results indicate an emanating power 


Tasie 6.—Lead age ratios of selected minerals 


Age ratio in 
Mineral Locality Geologic age 10° years 
Samarskite....... Glastonbury, Conn............. Pre-Triassic....... 270 
Pitchblende....... Jachymov, Bohemia............ Late-Paleozoic..... 220 
ae Giillhégen, Sweden............. Latest Cambrian. . . 400 
Bréggerite........ Karlhus, Raade, Norway........ Pre-Cambrian..... . 900 
Cleveite.......... Aust-Agder, Arendal, Norway.... Pre-Cambrian...... 1000 
Uraninite......... Keystone, S. Dakota............ Pre-Cambrian...... 1500 
Uraninite......... Sinyaya, Pala, Carelia, Russia... . Pre-Cambrian..... . 1850 


For geologic notes and references, see National Research Council, 1931, p. 7; Evans, et al., 1939, p. 
944; and Nier, 1939, p. 159. 


of only 0.15 per cent, which is hardly sufficient to account for any sig- 
nificant loss of radiogenic lead from the mineral due to diffusion of radon.*° 


LEAD TIME SCALE 


As we have seen in the introduction to this paper, the real benefit to 
geology of the radioactive methods of age determination results from the 
establishment of a quantitative time scale and its use for correlation 
purposes. In building up this time scale, the rocks or minerals measured 
must not only fulfill the four basic requirements but. in addition must be 
of well-defined geologic age. Unfortunately, only a very limited number 
of samples that have been measured satisfy these criteria for fixed points 
in the time scale. Essentially all the helium age ratios, prior to the recent 
work on magnetites, must be discarded for the reasons given. Most of 
the early lead age ratios can be considered only as maxima because of the 
failure to consider any correction for ordinary lead. Many other lead 
age ratios are of little or no value, because the measurements were made 
on badly altered material. Of the remaining radioactive minerals that 
apparently satisfy the four basic requirements, only the few shown in 
Table 6 can be considered as being even approximately well dated 
geologically. With the exception of the kolm, the span of possible horizons 
in these specimens is so great that they can hardly represent anything 
approaching fixed points in the time scale. These results clearly point out 


%° The possibility remains that the meteoric waters might remove a much higher fraction 
of the equilibrium amount of radon. 
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the need for more detailed study of the geological evidence concerning 
the age of radioactive minerals. For the present, therefore, the lead time 
scale derives its strength from the fact that a larger number of age ratios 
fall in the proper geologic sequence although individually they cannot be 
specifically assigned to definite horizons. Most of the minerals listed 
under the lead method in Figure 10 have been arranged on this basis. 
Although these values present a rather impressive array, covering a 
broad range, they cannot be used for comparative purposes with helium 
age ratios for rocks except in a very general way. 

If there were a clear-cut stratigraphic relationship between the geologic 
occurrence of a radioactive mineral, suitable for dating by the lead 
method, and a rock or mineral, suitable for helium age studies, an indirect 
but sharp comparison would be possible. Thus far these conditions have 
not been found. Instead, in general we are dependent upon separate 
geologic dating of the two materials. Hence, the comparison must be 
based upon the interdigitation of time scales containing numerous lead 
and helium age ratios covering a wide range of geologic time. It is the 
necessity for this type of comparison, plus the added complication that it 
is apparently quite difficult to assign geologic ages to radioactive minerals, 
that introduces so much uncertainty in the present status of radioactive 
age methods. 

ISOTOPE RATIO METHODS 


Because the two isotopes of uranium, UI and AcU, have different rates 
of decay, the proportions of the daughter elements and the end products, 
Pb**?/Pb?°*, have varied systematically during geologic time. This 
variation can be measured in two quite different ways which provide two 
other methods of estimating geologic ages. 

In discussing the lead method, the necessity for isotopic analysis of the 
separated lead was emphasized. Not only does such an analysis afford a 
measure of the amount of ordinary lead present, but it also provides a 
measure of the Pb?°?/Pb?°* ratio, which in turn gives an independent 
estimate of the geologic age of the specimen. Nier (1939, p. 159) has 
applied this method to a number of leads whose age ratios were also known 
from other data. These values are included in Figure 10. Among the 20 
measurements by Nier on leads from 13 different uranium-bearing min- 
erals, the Pb?°*/U?** ages are in fair agreement with those previously 
calculated from Pb/U ratios. The Pb?°?/Pb?°* ages on eight of these 
minerals are in good agreement with the Pb?°*/U?** ages, but serious dis- 
crepancies appear in the other five minerals. The differences are most 
likely due to alterations in the minerals not previously recognized as 
serious. 
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A distinct advantage of the Pb”°’/Pb?* ratio method is that it is much 
less sensitive to alteration of the radioactive mineral than the Pb/U ratio. 
Nier (1939, p. 162) has considered at length the possible effects of altera- 
tion and concludes “that with the exception of isolated cases any simple 
alteration will manifest itself in a difference in the ages computed by the 
two different methods.” Accordingly, where there is good agreement be- 
tween these two independent methods on the same specimen, it is very 
convincing evidence of their mutual reliability. 

The relative intensity of coloration of the Ac C and Ra C’ rings in cer- 
tain pleochroic haloes provides the other, more approximate, method of 
determining the ratio of past actinium series activity to uranium series 
activity and hence the age of the mica in which these haloes are found. 
Artificial coloration studies (Henderson and Bateson, 1934, p. 563) show 
that the blackening of the pleochroic haloes is proportional to the alpha- 
ray ionization in the early stages of blackening. Subsequently, the 
blackening increases more and more slowly with increasing exposure to 
ionization. For these reasons care must be taken in the selection of haloes 
to be studied, and in general it is necessary to use minerals whose age is 
greater than about 400 million years. These restrictions limit the ap- 
plicability of the method, and thus far only three micas have been studied 
quantitatively (Henderson, 1934, p. 591). Of these three specimens, only 
one, the Murray Bay biotite, has been dated by the lead method. The 
lead age ratio is 1070 million years, whereas Henderson obtained only 750 
million years. However, the latter value was based on an actinium : ura- 
nium series activity ratio of 0.04, Nier’s (1939, p. 150) more accurate 
value of 0.046 + 0.001 brings Henderson’s estimated age to 960 million 
years, in much better agreement with the lead value. 


RUBIDIUM-STRONTIUM METHOD 


Another method of geologic age determination has recently been sug- 
gested by Hahn and his coworkers (1937, p. 189). Ordinary rubidium is 
composed of two isotopes, Rb’? and Rb**, with abundances of 28 and 72 
per cent respectively (Nier, 1936, p. 1041). Only the heavier and less 
abundant isotope has been shown to be radioactive (Hemmendinger and 
Smythe, 1937, p. 1052). Rb** is beta-ray active and decays, with a half 
period of about 10"! years, to stable Sr’. In the application of the general 
age equation (1) to this method, the amount of radiogenic Sr*’ is the 
numerator, while the denominator is directly proportional to the rubidium 
content of the material. Ordinary strontium also contains Sr*’, but in an 
abundance of only 7.5 per cent. 

A handicap to the general applicability of this method is the rarity of 
rubidium. It occurs in small amounts in certain scarce lithium materials, 
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in microcline, and in rhodizite, a rare borate of Al, K, Cs, and Rb (Ells- 
worth, 1932, p. 12). The richest source of rubidium is lepidolite which in 
some cases contains as much as 3 per cent Rb.O. An additional limita- 
tion is the quantitative determination of the minute amount of radiogenic 
strontium present, even in the oldest rubidium-bearing minerals. Hahn 
and his associates separated 147 mg of 99.7 per cent pure Sr’? (Mattauch, 
1937, p. 189) from 1.011 kg of a lepidolite analyzed as containing 2.6 per 
cent 

At present the most serious uncertainty in the method is in the knowl- 
edge of the disintegration constant of rubidium. Evans (1940, p. 284- 
287) has recently reviewed the various measurements of \ rp» and selects 
3.5 X 10—?* sec—! as most nearly correct, although even this figure may be 
in error by 30 to 40 per cent. This value is based on Hahn’s measurements 
on the lepidolite from Silver Leaf Mine, near Winnipeg, Canada (not 
from Great Bear Lake, as reported), assumed to be the same age (1800 
million years from Pb/U) as the near-by Huron claim pegmatite (Ells- 
worth, 1932, p. 105). Thus the situation is one of circuitous relationship. 
Lead ages are used to estimate the decay constant of rubidium which in 
turn is used to measure other Rb-Sr ages. An accurate experimental de- 
termination of \ rp, while very desirable, is subject to large uncertainties 
owing to internal absorption and scattering of the soft beta-rays, E max 
= 0.13 Mev (Libby and Lee, 1939, p. 251), within even relatively thin 
sources of rubidium. 

The logical extension of the Rb-Sr method in the future is its applica- 
tion to a number of suitable minerals of widely varying geologic age, the 
interpretation of the results being one of over-all agreement with other 
age measurements. 


POSSIBLE POTASSIUM METHODS 


Potassium, like its heavier chemical analog, rubidium, is also beta 
radioactive. However, whereas no gamma rays have been observed for 
rubidium, potassium is known to emit a homogeneous gamma radiation 
of quantum energy 2 Mev, the ratio of the number of quanta to the num- 
ber of disintegration electrons being about 3 : 100. 

The beta-active isotope has been shown (Hemmendinger and Smythe, 
1937, p. 178-182) to be K*°, the least abundant of the three main isotopes, 
K**, K*°, and K*', present in fractional abundances of 0.932, 0.00011, and 
0.068 respectively (Nier, 1935, p. 283). Accordingly, with the emission 
of a negative beta-ray, K*° changes to Ca*®. Unlike rubidium, however, 


8t These values are from a personal communication from Professor O. Hahn to Professor A. C. 
Lane and are slightly different from those quoted by Hahn et al. (1937, p. 189). Since there 
is a possibility that the recovery of strontium was not complete, the 147 mg. represents a lower 
limit which Hahn thinks may be 10 per cent too small. 
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the Ca*® decay product of potassium is the predominantly abundant 
isotope (96.76 per cent) of ordinary calcium. 

Because potassium is so generously distributed in terrestrial materials 
and has a long half life (Tx*® = 1.06 + 0.3 X 10° years), it would be 
particularly suited for geologic age measurements. However, because 
ordinary calcium is commonly associated with potassium, it would be 
difficult to evaluate accurately the numerator of the general age equation 
as applied to the K*°->Ca* transmutation. It is conceivable that certain 
potassic minerals might originally have been quite free of calcium. These 
minerals would serve as valuable indices of age, but accurate isotopic 
analyses of the quantitatively separated calcium would be even more 
essential in these studies than in lead age measurements. 

Another possible age method based on the decay of potassium has been 
suggested by Evans (1940, p. 282), although no direct confirmation of 
this idea has yet appeared. It appears that the high energy, low intensity 
gamma radiation of potassium is not associated with the beta-ray disin- 
tegration of K**. Instead it is tentatively considered to be due to K- 
electron capture by a definite fraction of the K* nuclei. The decay 
product of this change would be A*°, the most abundant (99.63 per cent) 
isotope of argon. 

Postulating this to be the process, only about 7X 10-° ce of A*® (re- 
ferred to O° C and 760 mm Hg) would be formed in 100 million years. 
While it would be possible to determine such quantities of argon fairly ac- 
curately this method would be complicated by the likely presence of much 
larger amounts of ordinary argon, of nonradiogenic origin. Shepherd 
(1938, p. 311) has found argon to be a very common constituent of both 
lavas and plutonic rocks. He reports at least traces, which for his meas- 
urements means more than 10~‘ cc per g, for 25 out of 29 rock samples. 
These concentrations are several thousand times greater than could have 
been produced by the potassium present, which ranges from 0.11 to 4.46 
per cent. Thus, while both potassium age methods are worthy of further 
investigation, it does not appear likely that either of them will be capable 
of extensive application. 


CONCLUSIONS 


From a critical review of the fulfillment of the basic requirements by 
the various radioactive age methods, it is concluded that: (1) the re- 
tentivity of rocks for helium is the major source of uncertainty in the 
present application of the helium method, (2) analytical errors, earlier 
misconceptions, and the lack of attention to the basic requirements have 
necessitated the discarding of much of the age research of the past 3 
decades, (3) isotopic abundance measurements are essential in the de- 
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termination of ages by the lead method, (4) more specific geological 
dating is needed in the formulation of a useful lead time scale, (5) more 
direct comparisons between the lead and helium age ratios are urgently 
needed, (6) although other methods have been discovered, the lead and 
helium methods still remain the outstandingly promising radioactive 
methods of quantitative age determination, and (7) recent developments 
in the application of the helium method to selected common rock minerals 
promises a method of dating igneous rocks that should be of wide applica- 
tion as an aid to geologic correlations. 


GEOLOGIC NOTES ON ROCK SPECIMENS REFERRED TO IN TABLES 


4A, Chelmsford granite collected by N. B. Keevil from H. E. Fletcher Co. quarry, 
N. Chelmsford, Mass. Age: possibly Mississippian (Emerson, 1917). 

7A Rockport granite, Lat. 42° 40’ N., Long. 70° 35’ W., collected by N. B. Keevil 
from Upper Pigeon Hill Quarry, Cape Ann, Mass. Intrudes Salem gabbro 
diorite, which in turn cuts pre-Cambrian rocks and is possibly of Devonian age. 
It is cut by diabase dikes interbedded with Triassic sedimentary beds at Deer- 
field and Mt. Holyoke, Mass. Age: probably pre-Triassic and post-Cambrian. 

8A Rockport xenolith, fine-grained and cognate, from the same large hand speci- 
men as 7A. 

11A Granite collected by N. B. Keevil from quarry in Franklin, Me., (see Dale, 
1923, p. 219). Age: Middle Paleozoic. 

12A: Quincy granite collected by N. B. Keevil from J. 8S. Swingle quarry, Quincy, 
Mass. Age: possibly Mississippian, probably late Paleozoic. 


15A Fitchburg granite collected by N. B. Keevil from Rollstone Hill quarry, Fitch- 
burg, Mass. Age: probably Carboniferous. 

17A Granite collected by N. B. Keevil from Barryfield Quarries, 2 miles east of 
Kingston, Ont. Age: post-Grenville and pre-Keweenawan. 

23A Rhyolite collected by A. D. Howard 100 feet north of Tower Creek Bridge. 
Yellowstone Park. 

28A Obsidian collected by C. Max Bauer from Obsidian Cliff, Yellowstone Park. 

31A Lebanon granite collected by E. P. Kaiser (1938, p. 107) near Hanover, N. H. 

33A — sample collected by A. C. Lane, 1937, from Ruggles Mine, Grafton, 

1B__Porphyritic basalt collected by A. D. Howard from southeast base of Crescent 
Hill, Yellowstone Park. 

3B Crescent Hill Basalt collected by A. D. Howard 4 mile west of crossing of 
Geode Creek, Yellowstone Park. 


4B Basalt from the well-known Giant’s Causeway in Ireland, submitted as a large 
hand specimen for intercheck measurements by Arthur Holmes (1937a, p. 155) ; 
also see 5I and 71. 


6B Basalt from Lat. 40° 40’ N., Long. 40° 45’ W., collected from a flow by A.C 


Hawkins, 1936, 50 feet from the surface in a working quarry at Oldwick, N. J. 
The top of the flow has been removed by erosion, and the bottom is not visible. 
Similar flows in the area of slightly different periods of eruptivity have been 
identified as being of Triassic age from fossils in adjoining sedimentary beds 
(red shales and sandstones). Age: probably Triassic. 


7B Basalt collected by A. D. Howard % mile east of crossing of Geode Creek, 
Yellowstone Park. 
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Basalt collected by A. D. Howard from Overhanging Cliff, Yellowstone Park, 


Diabase from Lat. 43° 15’ N., Long. 79° 2’ W., collected by P. Price 3 feet from 
chilled contact of a north-south dike at the 17th level (1975 feet) of the Horne 
Mine, Norando, Quebec. In a previous publication (Evans, et al., 1939, p. 946) 
this rock was indicated as post-Cobalt in age. Subsequently, a personal com- 
munication from H. C. Cooke informs us that this dating is very uncertain, 
and the most that can be said is that this rock is post-Timiskaming (Cooke, 
et al., 1931, p. 143). 

Norite. Lat. 45° 1’ N., Long. 109° 25’ W., collected July 1937 by A. C. Lane 
from the metamorphosed Stillwater complex of Quad Creek, Mont., 1 mile 
north of state line. The Stillwater complex is bordered by a Cambrian con- 
glomerate at an erosional contact and is cut by granite and by basic dikes. 
Age: probably early pre-Cambrian. 

Trap rock, Lat. 44° 5614’ N., Long. 109° 2944’ W., elev. 9800 feet, collected from 
a small working quarry by A. C. Lane, July 1937, 2 feet from the margin of 
a 100-foot dike at Long Lake, Wyoming. The dikes of this series are intrusive 
into the early pre-Cambrian complex (Stillwater series), and appear to be 
older than the middle Cambrian rocks of Beartooth Butte, but the exact rela- 
tions are obliterated by overburden for a considerable distance. Age: possibly 
pre-Middle Cambrian and post early pre-Cambrian. 

Trap rock, Lat. 44° 58’ N., Long. 109° 29’ W., elev. 10,600 feet, collected by A.C. 
Lane, July 1937, as part of Beartooth Range suite, from center of 25-foot dike, 
dip 45 degrees, strike north-northwest. 

Monchiquite dike, Lat. 55° 45’ N., Long. 5° 3’ W., submitted by Arthur Holmes 
in 1936 from Riasg Buidhe on the Island of Colonsay, Scotland. Geologic re- 
lations nowhere provide evidence of age beyond the fact that it cuts formations 
regarded as Torridonian. Age: post-Devonian. 

Monchiquite dike, same as 26B but from Kilchattan instead of Riasg Buidhe. 


— collected by K. G. Bell (1939, p. 10) from Dalhousie Mt. in Nova 
cotia. 

Tower Falls basalt collected by A. C. Lane, 1937, in Yellowstone Park at about 
1.4 of Yellowstone Valley Conference (1937, p. 16), 134 feet vertically and 30 
feet horizontally from the outcrop surface. Age: Recent. 

Trap, Lat. 44° 5644’ N., Long. 109° 2914’ W., elev. 9800 feet, collected by A. C. 
Lane, as part of Beartooth Range suite, from altered glassy margin of Long 
Lake Dike by new road from Red Lodge to Yellowstone, near Long Lake, 
Mont., at about 35 of Yellowstone Valley Conference (1937, p. 11). Age: 
possibly latest pre-Cambrian. 

Trap rock from Gogebic Range, Mich., collected by L. M. Scofield, Jan. 16, 
1936. “Normal” from the 26th level of the Brotherton mine and cuts the upper 
phase of the Palms (Huronian). 

Granodiorite from La Motte, Quebec, collection of R. A. Daly. Other measure- 
ments by Keevil (1938a, p. 145). 

Tholeiite from Cleveland Dyke, Bolam County, Durham, Eng., submitted as 
a small hand specimen for intercheck measurements by Arthur Holmes (1937a, 
p. 118). 

Hornblende granodiorite, Lat. 48° 14’ N., Long. 77° 54’ W., collected by W. C. 
Gussow in 1936 from a specimen freshly blasted from a road out 10 feet from 
the surface. The granodiorite is intrusive into a sedimentary complex of 
Keewatin or possibly earlier age and is thought to be older than the diabase 
dikes in the area and the Cobalt series to the southwest. Age: probably pre- 
Cobalt and post-Keewatin. 

Quartz dolerite from Whin Sill, Scordale Beck, Westmoreland, England, sub- 
mitted as a small hand specimen for intercheck measurements by Arthur 
Holmes (1937a, p. 118). 
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4M Hornblende gneiss obtained from collection of R. A. Daly, designation 1737. 
Metamorphosed pre-Cambrian, known as Shuswap gneiss from British Co- 
lumbia. 

5M Feldspar augen gneiss submitted by R. A. Daly, designation 1848. 

6M Granodiorite from contact phase of Castle Peak stock, Long. 120° 50’ W., col- 
lected by R. A. Daly during survey of 49th parallel, designation 1442. 


7M Schist collected by A. C. Lane, 1937, as part of Beartooth Range suite. Speci- 
men is from a hornblende schist inclusion in the orthogneiss; some premeta- 
morphic material in it. Location about 34 of Yellowstone Valley Conference 
(1937, p. 11). 
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ABSTRACT 


In the helium method of age determination the age of the material to be tested is 
given by the ratio of the content of helium to the content of radioactive elements in 
the material. In an unaltered igneous rock, where presumably the minerals are of 
the same age, it would be expected that the distribution of helium in the various 
minerals should correspond with the distribution of the radioactive elements. In 
each mineral the ratio of helium to radioactivity should be the same. This was 
found to be not the case. 

Pyroxene and feldspar were separated from six samples of Triassic diabase. Age 
measurements on each yielded mean age ratios of 103 million years for the pyroxene 
samples and 36 million years for the feldspar samples. Six samples of late Triassic 
magnetite from ore occurrences in West Virginia, Pennsylvania, and Nova Scotia all 
gave age ratios closely grouped about a mean value of 134 million years. The con- 
stant age ratio, despite widely varying contents of radium and thorium, suggests that 
the magnetite has retained most of its helium. 

Additional evidence on the high helium retentivity of magnetite is shown by further 
work on magnetite specimens of various geological ages. These age ratios appear to 
show a sequence and spacing compatible with geological knowledge, and the results 
are in fair agreement with corresponding age ratios given by the lead method. 
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INTRODUCTION 


The present paper reports the progress made to date in a fundamental 
investigation of the retention of helium by rocks and rock minerals. This 
question constitutes the major uncertainty in the application of the helium 
method of geologic age determination. At the beginning of this re- 
search, in the fall of 1939, the following features characterized the status 
of the helium method: 

(1) Higher helium age ratios were observed for rocks of greater geologic age, as a 


general rule, but there were many exceptions. 
(2) Acidic rocks gave lower helium age ratios than basic rocks of the same geologic 


age. 
(3) In general, the time scale based on helium age ratios was considerably less than 
that given by the lead age ratios of radioactive minerals. 
(4) Anomalously high helium age ratios were observed for certain rock specimens. 
Preliminary helium age measurements on late pre-Cambrian rocks 
from the Sudbury district, Ontario, had suggested that this region might 
offer an unusually good opportunity for helium retention studies. Ac- 


TaBLe 1—Age measurements on Sudbury rocks 


10° years 
29B CO 3.29 0.096 1050 
Norite....... 2.68 0.55 151 
29A Micropegmatite.................+. 2.97 1.01 91 
9A ee 16.48 1.34 379 
10A; Creighton granite.................. 9.7 1.64 182 
15B; 11.3 0.26 1320 


*Specimens not analyzed for thorium. An average thorium content would reduce these age ratios 
by about one half. 


cordingly, a suite of type rocks were obtained and their helium and 
radium contents measured. The results, together with the calculated 
maximum age ratios, are summarized in Table 1. It was immediately 
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apparent from the extreme randomicity of the results that these rocks 
involved too many uncertainties to be of value for a fundamental investi- 
gation. Instead of continuing the work on this region, it was decided 
that the problem would best be furthered by restriction to rocks which 
involved as few uncertainties (or variables) as possible. For this pur- 
pose specimens of igneous rock were sought which fulfilled as many of 
the following prerequisites as possible: 


(1) No apparent metamorphism since time of formation. 

(2) Conditions of origin not in doubt, preferably with crystallization under low 
pressure to reduce possible inclusion of primary helium. 

(3) Bi-mineralic in composition and of sufficiently large grain size to facilitate 
mineral separation. 

(4) Opportunity for lateral and vertical range in sampling. 

(5) Geologic age known. 

(6) Sufficient radioactivity for accurate analysis. 

(7) Preferably glaciated surface to minimize effect of deep weathering. 


The Triassic Palisade diabase and Watchung flows of New Jersey 
appeared to be satisfactory material, and large, 15 to 25 pound, speci- 
mens were collected from quarries and road cuts over a lateral range of 
60 miles and up to 150 feet below the surface. Collection localities and 
petrographic notes for these specimens are included in a later section. 


HELIUM LOSS 

REGIONAL SCALE 
Before beginning detailed work on helium retention, a study was made 
of the helium content in various parts of the Triassic rock masses to dis- 
cover if there were any general relationship between helium content and 
depth of specimen below the surface, grain size, or proximity to the 
boundaries. Although the actual position of each specimen, with respect 
to the top or bottom of the flow or sill, was not closely known, it was 
assumed that a sufficient number of different positions were represented 
in the sampling to reveal any major variations existing. Of course, un- 
even distributions of radioactivity would cause similar uneven distribu- 
tions of helium, and thus other minor variations would be masked. Hence, 
only large-scale differences were sought. Table 2 presents the results of 
these analyses. It will be noted that a remarkably uniform distribution 
of helium in each rock mass is indicated. From these results, it appears 
that the helium retention within the rock mass is independent of the loca- 

tion with respect to the surface or boundary. 


CRUSHING AND GRINDING 


As indicated by Goodman and Evans (1941, p. 517), conflicting observa- 
tions have been made on the possible loss of helium during the preparation 
of the sample for measurement. Accordingly, a number of tests were 
made as a check on this possible source of error in the analytical technique. 
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Helium analyses were made on the same rock sample ground to sizes 
varying from 144 inch to — 150 mesh, with the results shown in Table 3. 
Within the experimental uncertainty, there is no evidence of loss of 
helium from any of the samples over the lengths of time indicated. 


Taste 2—Distribution of helium in rock masses 


Sample : Depth below Helium in 
No. surface in feet 10-5 
1st Watchung 
9B; 35 0.59 
2nd Watchung 
9B, 4 0.72 
3rd Watchung 
Palisade 
Standard Diabase 


The values in Table 3 are also indicative of the variation in multiple 
helium analyses. The particular specimen studied has the coarsest texture 
of any in the suite, and, as we shall see in the next section, shows the 
greatest difference in helium and radioactive content between its mineral 
components. Therefore, this specimen would be expected to show the 
maximum variation in the suite. The range of helium values observed 
for the samples of 4 inch material is only slightly greater than for the 
— 48 + 100 mesh samples, which suggests that the sampling uncertainty 
is less than the analytical uncertainty for these quantities of helium. In 
both cases the individual measurements are in agreement within about 
10 per cent of the mean value. 


AGE MEASUREMENTS 
LATERAL VARIATION 


Measurements of helium, radium, and thorium made on rock specimens 
from the Palisade sill over a lateral distance of 45 miles show a remarkably 
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uniform distribution. Age ratios calculated from these measurements * 
are given in Table 4 and are somewhat higher than the value obtained 
by Keevil and Goodman (Evans, et al., 1939, p. 942, Oldwick basalt) on 
one other specimen of New Jersey trap rock. While these recent de- 
terminations are definitely not in agreement with each other within the 


Taste 3.—Tests for helium loss on grinding 
Staten Island diabase 
(16B;) 


. P Sample weight Time* Helium in 

Bise of material (in grams) (in days) 10-5 cc/g 
5 68 1.03 
—48+100 3 0.92 
4 4 0.84 
—48+100 mesh................... 6 16 0.99 
7 17 0.80 


* Elapsed time between grinding and analysis. 


analytical accuracy of the measurements, there is no evidence of a 
systematic variation with lateral distribution, such as might be expected 
if any significant migration of helium occurred over a regional scale. 


SEPARATE MINERALS 


Radioactivity —Before considering the age ratios obtained for sep- 
arate minerals, mention should be made of the distribution of the radio- 
active elements in rocks. Waters? (1909, p. 677; 1910, p. 903) appears 
to have been the first to measure the radioactivity of mineral separates 
from common rocks. Although the gold leaf electroscope used by Waters 
was not sufficiently sensitive to detect any activity in the rock as a whole, 
his conclusions concerning certain mineral constituents have since been 
amply confirmed. The major activity was found to be concentrated in 
minute inclusions of zircon, allanite, and sphene, with somewhat less 
in magnetite, apatite, and rutile. The dimensions of the individual min- 
erals, sufficiently active to be measured by Waters’ instrument, were 
less than the ranges of most of the alpha rays. Hence, it is not surprising 
that these minerals were found to contain only a fraction of the helium 


tThe method of calculation is discussed by Goodman and Evans (1941, p. 507). 

2In view of recent developments, it is particularly interesting to note that Waters, in the first paper, 
stated the following: ‘‘The object of the investigation was to find in what minerals the radioactivity 
of igneous rocks is most concentrated, and, if possible to separate them with a view to their use as 
material for finding the age of the rocks in which they occur.” 
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content expected. It is regrettable that for this reason these interesting 
initial investigations were discontinued. 

Other workers, notably Strutt (1909, p. 298), Fletcher (1911, p. 102), 
and Holmes and Paneth (1936, p. 385), have made extensive measure- 
ments of the helium, radium, and thorium contents of zircons. The 


TasLe 4—Distribution of helium, radium, and thorium in the Palisade sill 


Sample i Helium in Radium in Thoriumin Age ratio in 
No. Location 10-5 cce/g 10-2 g/g 10 g/g 10° years 
16Bi Staten Island........... 0.91 ‘ 2.14 62 
16B2 North Bergen........... 0.96 0.19 2.37 72 
16Bs 0.93 0.16 2.34 75 
1.02 0.14 2.08 92 
16B; ee 1.01 0.1 2.60 76 


radioactive content varies widely; accessory zircon has in general high 
activity, while larger crystals in zirconium-rich rocks are often only feebly 
radioactive. 

Evans and Williams (1935, p. 441) have summarized the limited num- 
ber of measurements that had been made up to 1935 on other igneous rock 
minerals. They also report radium measurements on mineral separates 
of two Lassen dacites. Keevil (1938, p. 406) has subsequently studied the 
distribution of helium, radium, and thorium in the Quincy granite. 

It may be concluded in general that, although there are individual ex- 
ceptions, the radioactivity of biotite and other ferromagnesian minerals 
is higher than that of quartz and the feldspars. As is well-known to 
petrographers, these minerals usually contain most of the zircon, allanite, 
and sphene which may account for the higher radioactivity of the 
biotite and other ferromagnesions. There is considerable evidence that 
the radioactive elements also occur as individual atoms within the crystal 
structure of the essential rock minerals. 


Distribution of helium—As the radioactive elements undergo spon- 
taneous disintegration, the alpha particles (helium nuclei) are emitted 
with sufficient energy to penetrate the surrounding mineral matter a 
distance of 40 micra or less. After expending their energy, the alpha 
particles become normal helium atoms and appear to lodge themselves 
in the crystal structure of the mineral they have penetrated. This 
process may be pictured as an involuntary solid solution of the helium by 
the mineral. Unless the helium atoms escape from the solid structure and 
accumulate in open spaces, nowhere will radiogenic helium be present as 
gaseous inclusions. Since most rock minerals have dimensions much 
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greater than 40 micra, the major portion of the helium produced within 
a particular mineral is, initially at least, contained within the boundaries 
of that mineral. It might be expected that the firmness with which the 
helium atoms are held by a crystal structure would differ according to 
the atomic arrangement in the mineral. 


Taste 5.—Measurements on separated minerals 


FELDSPAR PYROXENE 

Sample Helium in Radiumin Thoriumin |} Heliumin Radiumin Thorium in 

No. 10-5 g/g 10-* g/g 10-5 10-" g/g 10 g/g 
Se 0.25 0.23 2.24 1.57 0.23 2.48 
0.31 0.18 1.76 1.46 0.18 2.51 
0.57 0.13 1.44 1.29 0.19 2.72 
Ear 0.35 0.12 1.82 1.37 0.14 2.36 
Sea 0.30 0.11 1.81 1.46 0.18 3.49 
eee 0.38 0.14 2.10 1.52 0.25 2.22 
Average...... 0.36 0.15 1.86 1.45 0.19 2.63 


If a particular mineral permits migration and loss of helium atoms 
from its structure, a lower age ratio would be obtained for this mineral 
than for an adjoining one with a higher helium retentivity. The observed 
age ratio of the rock containing these minerals would then be the 
weighted mean of the age ratios of the individual minerals and would 
depend upon the proportion and helium retentivity of all the minerals 
composing the rock. The true age of the rock would not be shown by 
direct analysis of the whole rock but could be found only from measure- 
ments on isolated minerals that have retained all their radiogenic helium. 


Age measurements.—The present research has shown that the age ratio 
in feldspar is much lower than in pyroxene separated from the same 
rock. Six specimens of Triassic diabase were ground sufficiently to unlock 
most of the grains, and the minerals were separated in an electromagnetic 
separator. The final samples consisted of pure feldspar separates and 
almost pure pyroxene concentrates. The mineral compositions of the 
pyroxene concentrates were calculated from Rosiwal analyses of the 
original specimens. 

The helium, radium, and thorium contents in the feldspar and in the 
pyroxene of each specimen of diabase are given in Table 5. The corre- 
sponding age ratios calculated from these value are compared in Table 6. 
The grouping of the age ratios for the two minerals clearly shows that 
there is a definite disparity. 
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The variation in the age ratios for the pyroxene is probably greater 
than would be the case if it could be separated from the feldspar more 
completely. Because of the extreme difference in helium age ratio, small 
quantities of feldspar would significantly affect the results obtained for 
the pyroxene fractions. 


TasLe 6.—Age ratios of separated minerals 


Age ratios of minerals in 10* years Age ratios of rock in 10® years 
Sample 
No. 

Feldspar Pyroxene Calculated* Observedt 
Aer 18 105 56 62 
28 109 82 72 
66 90 82 75 
__. 34 90 69 76 
35 102 92 98 
Average....... 36 103 79 79 


* Weighted mean from observed age ratios for separate minerals and Rosiwal analysis of their pro- 
portions. 
+ Observed on rock as a whole; see Table 4. 


A comparison of age ratios for the rock as a whole is also included in 
Table 6. The mean value, calculated from the age ratios for the sep- 
arate minerals and the Rosiwal analyses, in each case is in good agree- 
ment with that observed by direct measurement of the unseparated rock 
specimen. This correspondence furnishes additional evidence that there 
is no appreciable loss of helium upon crushing and grinding the rock 


sample. 
DISCUSSION OF RESULTS 


It seems probable that the age ratios shown by the feldspar, at least, 
are too low. If so, there are two reasonable ways in which such 
anomalous age ratios could have been caused; the feldspar has either lost 
helium or gained radioactivity during its geological life. 

In thin section the feldspar is seen to be fresh with all evidence against 
its having undergone growth or recrystallization since the formation of 
the diabase. Bombardment of alpha particles out of the more radioactive 
into the less radioactive minerals cannot explain the anomaly because 
this effect would make the helium content of feldspar too high rather 
than too low. For the same reason, the loss of minute radioactive grains 
during grinding and separation does not afford an explanation. 


te 


AGE MEASUREMENTS 553 


There is a possibility that the whole rock has become uniformly con- 
taminated late in its life by solutions containing uranium and thorium. 
At present this explanation does not seem likely, for unless there has been 
selective contamination of the feldspar no single quantity of uranium 
or thorium subtracted from the values observed for each mineral will 
bring the age ratios into agreement and at the same time keep them 
within reasonable limits for the Triassic age. It is also quite unlikely 
that contamination could be so uniform for the entire rock body. 

The most plausible explanation, at present, is that a portion of the 
helium in some way has escaped from the feldspar. Keevil (1940, p. 311) 
has presented a theoretical consideration of the probability of migration 
of helium atoms held in the crystal structures of common rock minerals. 
He concludes that an appreciable helium loss over geologic periods 
of time is highly improbable. However, since helium leakage appears 
as the most likely explanation of the age ratios observed, it seems prob- 
able that Keevil’s theoretical treatment fails to include all the factors 
involved. 

Assuming that helium does escape from minerals, and in varying de- 
grees for different minerals, it is possible that the fractional loss may be 
a fairly specific quantity for each mineral, even with widely varying 
helium concentrations. This state of affairs would afford an explanation 
of the grouping of the age ratios about the mean values observed for the 
two minerals in Table 6. 

TRIASSIC MAGNETITES 

Although the foregoing measurements indicated helium loss from 
feldspar, the evidence was inconclusive concerning the degree of reten- 
tion of helium by the pyroxene. It seemed probable that information 
concerning this question could be obtained from age measurements on 
other Triassic minerals, since, if higher age ratios were observed for 
these minerals than for the pyroxene, helium loss from the latter would be 
suggested. For this purpose four specimens of magnetite from deposits 
associated with the Triassic intrusives in West Virginia and Pennsylvania, 
together with two specimens associated with Triassic flows in Nova Scotia, 
were obtained. The results of helium, radium, and thorium measure- 
ments and the calculated age ratios for these six magnetites are given in 
Table 7. The age ratios range in value from 126 to 145 million years 
with a mean of 134 million years. The six values are in good agreement 
within the analytical uncertainty of the individual measurements. Com- 
paring these age ratios with those for the pyroxenes in Table 6, it is 
seen that the magnetite values are all larger than even the largest value 
of 119 million years observed for 16B,, and the mean value for the 
magnetites is about 30 per cent greater than the mean value for the 
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pyroxenes. Thus it appears likely that the pyroxene also has lost part 
of its helium but that the fractional loss from pyroxene is smaller than 


that from feldspar. 
In addition to the numerous interchecking measurements discussed in 


the papers by Evans and Goodman, the agreement in the magnetite 


Taste 7.—Age ratios of Triassic magnetite deposits 


Sample Helium* in Radium*in Thorium*in Age ratio in 
No. 10*cc/g  10-%g/g 10-7 g/g years 
1 22.4 31.6 2.44 126 
2 Goose Creek, Va........... 15.7 14.8 2.25 137 
3 Boyertown, Pa......... 5% 17.6 27.8 1.32 133 
4 French Creek, Pa - 8.6 8.5 1.07 145 
5 Gerrish Mt., N.S 7.8 4.9 15.4 128 
6 11.4 18.0 8.2 135 


*The units used are smaller than in the preceding tables in order to make the variation between 
samples more apparent. 


age ratios also affords evidence that the radium and thorium determina- 
tions are accurate on an absolute basis. The proportions of radium and 
thorium in these specimens varied over such a wide range that the age 
ratio for No. 5, the Gerrish Mt., N. S., specimen, is essentially a He/Th 
ratio while the age ratio for No. 6, the Lakeville, N. S., specimen, is essen- 
tially a He/Ra ratio, with intermediate variations for the other specimens. 
If any systematic error occurs in either the radium or thorium measure- 
ments, it should be reflected in a lack of agreement in the corresponding 
age ratios for these specimens. 


MAGNETITES OF VARIOUS AGES 


While the age ratios for Triassic minerals indicated that magnetite has 
a higher retentivity for helium than either feldspar or pyroxene, still it 
could not be concluded that magnetite retains all its helium. Accordingly, 
a new program of research was initiated to investigate what may be con- 
sidered the specific retentivity of common minerals. As mentioned above, 
there is a possibility that this factor may be fairly constant for indi- 
vidual mineral species. As indicated by Goodman and Evans, no reliable 
helium time scale existed at the time the present researches were begun. 
Accordingly, as a tentative measure, it was decided to consider magnetite 
as having complete helium retentivity and to use it as a standard against 
which other minerals could be compared. For example, considering 
magnetite as unity, the above results indicate that the specific retentivity 
of feldspar is about 0.3 and that of pyroxene about 0.8. By testing other 
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mineral types in addition to further measurements on those already 
studied, it is hoped that some information may be gained on the absolute 
retentivity of helium in magnetites as well as the possible use of other 
minerals for age determinations. 


TaBLeE 8.—Age ratios of various magnetite deposits 


Sample Heliumin Radiumin Thoriumin Ageratioin 
No; Specimen and geologic age 10-* ce/g 10-" g/g 10+ g/g pena 
Mid-Tertiary 
7 Great Western Mine, Utah.... 0.85 0.62 4.0 26 
8 Black Magnetic Mine, Utah... 0.98 0.67 1.4 37 
Early Tertiary 
Mémé, T. N., Haiti.......... 0.72 0.33 0.31 60 
10 Fierro, New Mexico.......... 0.40 0.12 0.9 59 
Early Cretaceous 
11 Prince of Wales I., Alaska..... 3.53 1.23 0.64 81 
12 Lynn Valley, Br. Columbia... . 1.26 0.42 0.44 88 
Late Jurassic 
Texada I., B. C. 
13 M. LT. apecimen.........% 1.66 0.34 0.50 130 
14 Smithsonian Inst. specimen. . 2.35 0.51 0.22 131 
Late Triassic-Early Jurassic 
1 2.24 0.32 0.24 126 
2 VO... 1.57 0.15 0.22 137 
3 Boyertown, PS... . 1.76 0.28 0.13 133 
4 French Creek, Pa............ 0.86 0.085 0.11 145 
5 Gerrishn 0.78 0.049 1.54 128 
6 1.14 0.18 0.08 135 
Carboniferous 
15 Magnitnaya, Urals........... 4.73 0.44 0.27 300 
Devonian 
16 0.81 0.048 0.33 316 
Blagodat, Urals 
17 1.48 0.099 0.23 366 
18 Smithsonian Inst. specimen. . 2.60 0.158 0.42 396 
Pre-Cambrian 
19 64.9 0.300 16.2 1100 
20 26.6 0.600 0.51 1100 


Nore: The age ratios for samples 19 and 20 include a correction for the decay of uranium and thorium. 
(See Figure 6 of Goodman and Evans, p. 507.) The exact agreement is merely a coincidence. 


The first step in this program has been the establishment of mag- 
netite—helium age points in the geological time scale, against which the 
age ratios of other minerals may be compared. The results of this 
work to date are given in Table 8. The age ratios obtained range from 
26 million years for a Mid-Tertiary magnetite to 1100 million years for 
two pre-Cambrian magnetites. Except for the two Late Jurassic speci- 
mens, for which the geologic evidence is somewhat uncertain, the age 
ratios are in the proper geologic sequence. Some evidence that magnetite 
retains most, if not all, of its helium is given by the fact that the magnetite- 
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helium time scale is in fairly good agreement with the post-Proterozoic 
points given by lead ratios and isotope ratios for radioactive minerals 
(Goodman and Evans, 1941, Fig. 10). 


CONCLUSIONS 


It is concluded from the investigations to date that: 

(1) Age determinations by the helium method made directly on rock 
samples are likely to be in large error owing to loss of helium from some 
or all of the mineral constituents. 

(2) Mineral species differ in their ability to retain helium. Feldspar 
loses a large proportion of its helium; pyroxene retains more of its helium 
than feldspar but less than magnetite. 

(3) From the general agreement between the helium age ratios shown 
by magnetite and the existing time scale indicated by the lead method, 
magnetite appears to retain most, if not all, of its helium. 


GEOLOGIC NOTES ON ROCK SPECIMENS REFERRED TO IN TABLES 


29B Sudbury gabbro. 10 feet below original surface in road cut. Lot 3. Concession 
III, Graham Twp., Ontario. 


20B: Sudbury norite. 15 feet below original surface in railway cut. Lot 1. Con- 
cession IV, Wisner Twp., Ontario. 


29A Sudbury micropegmatite. 6 feet below original surface in railway cut. Lot 1. 
Concession III, Wisner Twp., Ontario. 


9A Murray granite, Sudbury District. 5 feet below surface in railway cut. Lot 10, 
Range V. McKim Twp., Ontario. 


10A; Creighton granite, Sudbury district. 5 feet below surface in road cut. Lot 9. 
Concession I, Snider Twp., Ontario. 


26A Late aplite phase of the Creighton granite, Sudbury district. 4 feet below 
surface. Lot 1. Concession III, Snider Twp., Ontario. 


15B; Olivine diabase, Sudbury district. 2 feet below surface. Lot 10. Concession V, 
McKim Twp., Ontario. 


9B: Second Watchung basalt. Murray Hill, N. J. 25 feet below surface in quarry. 

9B. First Watchung basalt. Summit, N. J. 30 to 50 feet below surface in quarry. 

9B; First Watchung basalt. Scotch Plains, N. J. 35 feet below surface in quarry. 

9B, First Watchung basalt. Same quarry as 9Bs, 70 feet below surface. 

9B; First Watchung basalt. Chimney Rock Quarry, north of Somerville, N. J. 50 
to 100 feet below surface. 

9B. Second Watchung basalt. Livingstone, N. J. 5 feet below surface in road cut. 

9B: — Watchung basalt. Preakness Mtn., N. J. 4 feet below surface in road 
cut. 

9B, Third Watchung basalt. Riker Hill, N. J. 4 feet below surface in road cut. 

9B. Third Watchung basalt. Mountain View, N. J. 7 feet below surface in road 
cut. 

9B: Third Watchung basalt. Packanack Mtn., N. J. 15 feet below the surface in 
old quarry. 
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16B, Palisade diabase. Staten Island. 8 feet below surface in quarry. 


16B: 


16Bs 


16B, 


16Bs 


18B 


Rosiwal analysis (per cent by volume): Feldspar 62 

Pyroxene 36 

Magnetite 1 
Coarse-grained, fresh diabase. Feldspar entirely fresh, in labradorite laths with 
small amount of more acid intergrowth. The pyroxene occurs in large anhedral 
crystals, interstitial to the feldspar and commonly poikilitically including 
feldspar laths. Most of the pyroxene is fresh; where altered it has been re- 
placed by a strongly pleochroic green biotite with a little chlorite and magnetite. 
Its composition is in a range between augite and pigeonite with an optic angle 
between 40 and 50 degrees. 


Palisade diabase. River Road, North Bergen, N. J. 20 feet below the surface 
in road cut. 
Rosiwal analysis (per cent by volume): Feldspar 44 

Pyroxene 53 

Magnetite 3 
Fine-grained diabase. All minerals are fresh. The pyroxene occurs in small 
rounded grains between the feldspar laths. 


Palisade diabase. River Road, Guttenberg, N. J. 150 feet below surface in 


quarry. 
Rosiwal analysis (per cent by volume): Feldspar 49 

Pyroxene 48 

Magnetite 2-3 
The rock has a fine, holocrystalline, diabasic texture with occasional large 
crystals of pyroxene. A small amount of acid feldspar occurs in graphic inter- 
growth between the labradorite laths. The pyroxene forms rounded and irregu- 
lar grains. All minerals are fresh. There are two pyroxenes present, one with a 
composition in the range of pigeonite with an optic angle of about 15 degrees, 
and the other between pigeonite and augite. 


Palisade diabase. Edgewater, N. J. 70 feet below surface in road cut. 
Rosiwal analysis (per cent by volume): Feldspar 47 
Pyroxene 50 
Magnetite 3 
Fresh, fine-grained diabase. Feldspar laths show a minor amount of strain. 
The pyroxene is mostly pigeonite with optic angles close to 10 degrees. 


Palisade diabase. Kingston, N. J. 20 feet below surface in quarry. 
Rosiwal analysis (per cent by volume): Feldspar 54 

Pyroxene 45 

Magnetite less than 1 

Biotite less than 1 
Medium-grained, fresh diabase. Feldspar mostly labradorite laths with some 
more acid intergrowths occurring interstitially. Pale-green pyroxene toward 
pigeonite in composition. 
Standard Triassic diabase, Virginia. A number of rock samples have been pre- 
pared under the supervision of the National Bureau of Standards, for intercheck 
of radium and thorium measurements in various laboratories working on the 
radioactivity of terrestrial materials. One of these samples, the Triassic diabase, 
was used in the helium-age investigation. 


1. Magnetite ore. Cornwall, Pa. Specimen of magnetite and maghemite with small 
amount of greenish gangue. 


2. Magnetite in basic pegmatite. Goose Creek, Loudoun County, Va. Specimen of 
pegmatite with coarse, bladed augite crystals and basic feldspar. Considerable 
amount of coarse magnetite separated from it by electromagnet. 

Age: Palisade. The pegmatite occurs in a diabase dike. 


3. Magnetite ore. Boyertown, Pa. Specimen of fine-grained magnetite with small 
amount of greenish gangue. 
Age: Palisade. Ore deposit associated with late Triassic diabase. 
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4. Magnetite ore. French Creek, Pa. Specimen of medium-grained magnetite with 
chloritic gangue and small amount of pyrite. Magnetite separated by electro- 
magnet. 

Age: Palisade. Deposit associated with Triassic intrusives. 


5. Magnetite. Gerrish Mountain, Nova Scotia. The specimen was taken from a 
magnetite-quartz veinlet in a Triassic flow. 
Age: Palisade. 


6. Magnetite. Nova Scotia. 2 miles northwest of Lakeville. Associated with a 


Triassic flow. 
Age: Palisade. 


7. Magnetite ore. Great Western Mine, Iron Springs, Iron Co., Utah. Smithsonian 
Institution Collection No. 42174. 
Age: Probably Early Miocene. Deposits associated with post-Eocene andesite 
which is similar to Miocene lavas of Wasatch Mtns. (Leith and Harder, 1908, 
p. 46.) 

8. Magnetite ore. Black Magnetic Mine, Irontown, Utah. Smithsonian Institution 


Collection No. 42177. 
Age: Probably Early Miocene. See Specimen No. 7, Great Western Mine. 


9. Magnetite and calcite. Germain Tunnel, Mémé T.N., Haiti. Specimen consists 
of large masses of foliated magnetite in a gangue of coarsely crystalline calcite. 
The magnetite is apparently pseudomorphic after hematite. 

Age: Tertiary. Deposits associated with intrusions cutting early Tertiary beds. 
(Geological Survey, Republic of Haiti, p. 443.) 


10. Magnetite ore. Fierro, New Mexico. Specimen is unweathered but appears to 
have suffered from shear stress. 

Age: Early Tertiary(?). Deposit associated with porphyry belonging to a group 
of intrusive rocks placed as latest Cretaceous or, more probably, earliest Tertiary 
by Lindgren, Graton, and Gordon (1910, p. 35). 


11. Magnetite in copper ore. Copper Mtn., Prince of Wales I., Alaska. Chalcopyrite 
and magnetite with gangue. Magnetite separated by electromagnet. Smithsonian 
Institution Collection No. 91562. 

Age: Early Cretaceous(?). 


12. Magnetite ore. Lynn Valley, Vancouver, British Columbia. From Copper Queen 
Claim. 
Age: Probably Early Cretaceous or late Jurassic. 


13. Magnetite. From the contact of the iron and copper deposits on Texada Island, 
British Columbia. The magnetite was separated from a greenish gangue by elec- 
tromagnet. M. I. T. Collection No. 5943. 

Age: Probably Nevadan. 

14. Magnetite. Texada Island, British Columbia. Same place and same age as 
Sample 13, but different collector. Smithsonian Institution Collection No. 42200. 


15. Magnetite and hematite. Magnitnaya, Urals. Massive ore, somewhat oxidized. 
Magnetic, but with reddish color and streak. 

Age: Probably Upper Carboniferous or Permian. Associated with granitic intru- 
sions cutting Lower Carboniferous strata (Zavaritsky, 1937, p. 92). 

16. Magnetite ore. Vyssokaya, Urals. Magnetite separated from greenish gangue by 
electromagnet. 

Age: Devonian(?). Same age as Specimen No. 17, Blagodat, Urals. 

17. Magnetite ore. Blagodat, Urals. Small amount of greenish gangue separated 
from magnetite. 

Age: Devonian(?). Deposit associated with syenite intrusion believed by Derwies 
(1937, p. 39) to be Middle Devonian. 

18. Magnetite ore. Blagodat, Urals. Smithsonian Institution Collection No. 103318. 
Age: See Specimen 17. 
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19. Magnetite ore. Kiruna, Sweden. Fine-grained magnetite, separated from a 
gangue of apatite about 5 per cent by weight. 
Age: Pre-Cambrian. 


20. Magnetite ore. Carufel Copper Prospect, Hart Twp., Ontario. Fine-grained 
magnetite with a considerable amount of greenish gangue. 
Age: The ore occurrence is believed to be of a contact metamorphic type asso- 
ciated with intrusive diabase similar to the dikes thought to be of Keweenawan 
age in the Sudbury district. 
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ABSTRACT 


In Part I observational data and hypotheses pertinent to the cooling problem of the 
earth are reviewed. The possible importance of heat transport by slow thermal 
convection currents is indicated. In Part II the theory of the cooling of a solid 
radioactive earth is presented by a simple method which identifies this problem 
with an equivalent one in the cooling of a nonradioactive earth. The facts as to 
slowness of cooling for temperature perturbations of long wave length are reviewed. 
The importance of the large heat capacity of the earth is emphasized by means of 
simple examples. If the earth is cooling the rate must be very small, but, if heating 
occurs, higher rates are possible. Thermal considerations indicate that solidification of 
the mantle probably occurred from the bottom up. Because the steady-state tempera- 
tures in a radioactive earth are significant in determining whether the earth is heating 
or cooling at depth, such temperature distributions have been computed for a dozen or 
more types of radioactive earths. It is probable that thermal equilibrium in the earth 
as a whole is far from achieved. The earth may be heating everywhere—cooling every- 
where—or heating at some depths and cooling at others. The transient temperatures 
are expressed, and the transient heat flows computed for a number of types of suitable 
radioactive earths which are amenable to convenient computation. The case of heat 
flow in layered earths is especially examined. The observed surface heat flows fix 
definite upper limits to the total heat generation which may exist within about 100 
or 200 kilometers of the surface. It is important to improve the reliability and 
scope of present information concerning (1) surface heat flows, (2) the structure 
and composition of the crust, and (3) the association of heat generation with crustal 
layers. Until better knowledge of all three of these factors in the same region is 
available, estimates of the amount of radioactivity existing in and below the crust 
will involve large uncertainties. Even granting excellent thermal information about 
the crust, it is possible to deduce by thermal theory only very limited information 
about regions below a depth of 200 or 300 km. Significant concentrations of radio- 
activity may exist all the way to the center without violating known criteria. At 
depth the problem of the cooling of the earth is nonunique to such a degree that 
little information about the deep interior may be deduced from purely thermal data. 
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INTRODUCTION 


Since Lord Kelvin’s famous controversial estimate of the age of the 
earth in 1862, much attention has been given by geologists, physicists, 
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and mathematicians to the many aspects of the thermal problem of the 
earth. Jeffreys has published an excellent survey of modern develop- 
ments in this subject in his book, The Earth (1929). The contributions 
of Ingersoll and Zobel (1913), Adams (1924), and Holmes (1933) are 
of much interest. Recently Spicer (1937) and Van Ostrand (1940) have 
computed and published tables relating to the cooling of homogeneous 
nonradioactive and radioactive earths. Lowan (19383) has contributed 
to the theory of the cooling of general types of nonhomogeneous radio- 
active earths, and Gutenberg (1939) has briefly summarized the present 
status of the subject in a chapter in his recent book. The limitations 
of space prevent the author’s making adequate acknowledgment, or 
even reference, to the copious literature in this field which has developed 
since Kelvin’s stimulating work. The appended bibliography is typical 
rather than comprehensive. 


PART I. REVIEW OF HYPOTHESES AND OBSERVATIONAL DATA 
INITIAL STATE OF THE EARTH 


No accepted theory of origin.—It is usually postulated that the earth 
had a fiery beginning—more specifically that it was formed by the 
collision or near collision of a passing star (or two stars) with the sun. 
For many years attempts have been made to prove that such a cataclysm 
could actually account for the major features of the present solar system. 
However, these efforts have never been altogether successful, and recent 
work by Spitzer (1939) has revealed seemingly fatal objections to all 
solar collision theories of the origin of our planetary system. Irrespective 
of theories of origin (and without complete justification), it is generally 
held that the earth was originally hot and entirely in the molten state. If 
so, it may be established (see for example Jeffreys, 1929, p. 79) that a crust 
would cool to the point of solidification within a few tens of thousands of 
years, and that soon thereafter would have cooled sufficiently for an 


ocean to condense. 


Constancy of surface temperature——As noted by Kelvin, Jeffreys, and 
others the surface temperature would quickly reach an equilibrium essen- 
tially established by a balance between solar radiation and that re- 
radiated to space—the heat conducted to the surface from below was 
always entirely too small to affect appreciably the surface temperature.* 
Assuming constancy of the sun’s heat this mean temperature would not 
differ appreciably from that observed today. Furthermore it is well 


1 Today solar radiation supplies an average of about 4.2 x 10-8 cal./em®sec. to the surface. Conduc- 
tion from below, probably about 1.2 x 10-° cal./cm*sec. 
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established from independent geological evidence.that the mean surface 
temperature of the earth has remained remarkably constant? as far 
back in time as the record can be traced—namely, for more than half 
the earth’s lifetime. This essentially constant surface temperature is 
conveniently adopted as the zero of reference. 


POSTULATE OF A RIGID MANTLE 


Three major divisions in the earth are recognized by geophysicists: 
The heavy core (mean density, 10.7) has a radius of about 3400 kms., 
or slightly over half that of the earth. It appears to be fluid and it is 
sometimes assumed to be composed of nickel-iron. Overlying the core 
is the thick lighter mantle (density, 3.4 to 5.7) generally held to have the 
approximate composition of ultrabasic rocks, but becoming richer in iron 
toward the core. Above the mantle is the thin crust consisting of several 
layers of lighter rocks whose total thickness is perhaps 20 to 60 kms. 
The crust appears to be subject to regional variations, and its structure 
differs especially between continental areas and the basin of the Pacific 
Ocean. There are minor structural features of interest in both the crust 
and the mantle, but the two major discontinuities in the earth occur 
between the crust and the mantle and between the mantle and the core. 
The seismological evidence supports conclusions based upon the tidal 
yielding of the solid earth in assigning to the mantle a coefficient of rigidity 
somewhat greater than that of steel. But the stresses involved in these 
seismic and tidal phenomena are of very short period relative to geologic 
time intervals, and it is not yet clearly established what effective 
rigidity the mantle may possess toward stress differences prolonged in 
geological time. However, it will be assumed in the thermal discussions 
which follow that the mantle, or at least the major part of it, is rigid, 
and has been rigid since its formation—“rigid” here meaning that convec- 
tion velocities have been small in the sense of a paragraph which follows, 
so that conductive heat transport predominates. It is also assumed that 
the rigid mantle was formed so rapidly that its initial temperature was 
essentially independent of the distribution of radioactivity. In the ab- 
sence today of any accepted theory of the origin of the earth it is perhaps 
futile to inquire further into the genesis of the primeval rigid mantle and 
fluid core here postulated. Nevertheless, because of the prevalence of the 
assumption that the earth was originally thoroughly molten, it is perhaps 
of interest to recall here a process by which the postulated earth might 
have been derived from a previously molten state. 


2From the present point of view, variations in surface temperature of several tens of degrees are 
insignificant, since we are dealing with internal temperatures measured in thousands of degrees. 
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SOLIDIFICATION OF THE MANTLE 


The different possible modes of solidification of an originally molten 
earth were well discussed by Hopkins (1839) over 100 years ago. In 1890 
Kelvin outlined simple thermodynamical principles which are basic to 
modern discussions of the solidification process, such as those of Adams 
(1924) and Jeffreys (1929), but Kelvin lacked requisite thermal data 
about rocks upon which to develop his projected conclusions, and Adams 
appears to have first discussed upon a quantitative basis the process of 
solidification at shallow depths. 

A thoroughly molten initial stage of the earth is here assumed to have 
existed, with temperatures sufficiently high to permit a free and rapid 
convection. In this stage, differentiation by density is presumed to have 
occurred to such an extent that the present heavy, molten core separated 
from the lighter material of the mantle, with the formation of a major 
convection system in the mantle and possibly another one in the core.® 

Cooling in the mantle was essentially by convection, and accordingly 
the temperature distribution there was almost the adiabatic distribution. 
As the mantle rapidly cooled by convection solidification commenced 
at the depth at which the curve representing the adiabatic temperature 
first intersected, or became tangent to, the curve representing the depth- 
freezing point temperature. The adiabatic temperature gradient in the 
crustal magma is determined primarily by the ratio of the volume co- 
efficient of thermal expansion to the specific heat in accordance with 
the equation (Adams, 1924): 


The value of the gravitational acceleration g remains nearly constant 
throughout the mantle. Near the surface the absolute temperature V 
may be taken as about 1200°, the specific heat, Cp, as .25 cal./deg.gr., or 
nearly 1X10" ergs/deg.gr., and the thermal expansion, a, as 4X<10-°. 
With these values, the adiabatic gradient is about 4.6 X 10-*° C./cm. The 
behavior of the adiabatic gradient at greater depths is very uncertain of 
estimate. Bridgman (1931, p. 136) states that the thermal expansion for 
liquids (not, however, silicates) at 12,000 atmospheres is generally be- 
tween 14 and \& its value at atmospheric pressure. Similarly, the value 
of C, in liquids generally first decreases slightly with pressure, then 
increases until at 12,000 atmospheres it has nearly recovered its original 
value; as a function of temperature the specific heat increases for nearly 
all rock minerals from a value of .2+ .02 at 100° C. to about .3+ .03 at 
1000°C. Because of the decrease of the thermal expansion and the in- 


3It is immaterial for present purposes whether convection existed in the core, or not 
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crease of the specific heat, it is probable that the adiabatic gradient 
decreases at depth. 

The form of the depth-freezing curve is also difficult of estimation. 
The early work of Barus (1892) on basalt leads to estimates, based 
upon the Clausius-Clapeyron relation for the change of freezing point 
with pressure, of a freezing point gradient about tenfold greater than 
the adiabatic, or 3 X 10-°°C./em. This, if true, means that freezing 
first commenced at depth. Barus’ work, however, was done on anhydrous 
samples which are certainly not representative of crustal magmas. 
Furthermore, the simple form of the Clapeyron relation does not apply 
for multiple component systems. The recent work of Goranson (1938) 
on the equilibrium of silicate-water systems under the temperatures of 
800°-1200°C. and pressures up to 4000 bars (corresponding to depths 
of somewhat over 10 kms. in the earth) well indicates the complexity 
of equilibria relationships, even in the simple systems studied. In his 
experiments, the freezing point temperatures decreased instead of in- 
creased with increasing pressure throughout the observed pressure range, 
0-4000 bars, but, as Goranson points out, the behavior at higher pres- 
sures of the parameters in the generalized Clapeyron equation which 
applies is such as to reverse the trend and lead to eventual increase of 
freezing point with pressure. Bridgman (1931, p. 206) in commenting 
upon the effect of pressure upon melting point says, 


“we have to take the present experimental evidence as making it exceedingly 
probable that the melting curves of all substances rise indefinitely, at a continuously 
decreasing rate, until pressures are reached so high as to introduce entirely new 


phenomena .. .” 

There exists, unfortunately, little experimental evidence on this subject 
about silicates or magmas. Knowledge of equilibria in silicate systems 
is still confined to relatively simple systems, at relatively low pressures 
and temperatures. 

Jeffreys (1932) revised his suggestion (1929) that solidification in the 
mantle commenced at a depth of several hundred kilometers, and sup- 
ported the hypothesis that the process commenced at the outer boundary 
of the core itself. Such an assumption follows naturally if one accepts 
the hypothesis that the mantle is today predominantly solidified and 
rigid. Otherwise the cooling of the mantle to its present degree of rigidity 
from the postulated molten state would present difficulties. This subject 
will be further discussed in Part II. If solidification in the mantle pro- 
ceeded from the bottom up, the magma would freeze at a relatively rapid 
rate, leaving the core essentially at its original temperature, and insulated 
by a rapidly thickening shell of silicates. During the process the surface 
of the earth would remain molten and would dissipate large amounts of 
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heat by radiation. Both gravitational differentiation and differentiation 
by selective freezing would be effective during the process. It is esti- 
mated that by such a process the complete mantle with its initial tempera- 
tures essentially at the freezing point would have formed in less than 
20,000 years. The heat loss thenceforth would be by conduction and at 
a rate vanishingly small relative to that of the earlier convective stage. 
As a consequence of the rapid freezing the initial temperatures in the 
mantle would in no way be expressive of the equilibrium temperatures 
finally to be established by the radioactivity trapped in the rocks. At 
depth in the mantle the subsequent course of the temperature might be 
either increasing or decreasing as controlled by the concentration of the 
radioactive heat sources. Evidently, much freedom is permitted for 
further hypotheses concerning the subsequent thermal and mechanical 
history of such a mantle. 

The outlined process of solidification involves selective freezing and 
possibly also some segregation by density. It is, therefore, probably not 
a reversible process, and any subsequent remelting of the solidified mantle 
would require higher temperatures—how much higher is, of course, un- 
known. Thus remelting of the mantle by radioactive heat would require 
heating through the necessary temperature increment, plus supplying the 
latent heat of fusion, which for rocks under ordinary conditions is ap- 
proximately 100 cal./gr. 

The heavier core was well insulated thermally by an enormous layer 
of relatively poorly conducting material. If originally liquid, it no doubt 
has remained liquid to this day and has probably maintained a nearly 
uniform temperature by convection or conduction or both. Its present 
temperature probably differs from its original temperature essentially by 
the amount of the radioactive heat which may have been produced in it. 

As the mantle formed, its augmenting thickness would increasingly 
favor a tendency for the initiation of slow convection currents in the 
growing shell. Their possible contribution to the processes of heat transfer 
in the earth is important and will be considered in a following paragraph. 
However, the existence of such currents is still a question for research. If 
they exist, present theories for the cooling of the earth, which have all 
been based upon the hypothesis of a solid earth, will need radical revision. 


RADIOACTIVITY OF ROCKS 


It has long been recognized that the radioactive heat being generated 
in rocks is large enough to be a major influence in the earth’s thermal 
history. The problem of estimating the amount of radioactivity which 
may exist at depths of several hundred kilometers is still little under- 
stood, but the approximate determination of the average radioactive con- 
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tent of the bulk of the crust seems within the scope of present possibilities. 
A beginning in this problem is just now being made. During the last 
several years an extensive and informative program for the study of 
radioactivity of rocks has been undertaken by Evans and his associates. 


TABLE 1.—Radioactive heat generation 


In Rocks 


(After Evans and Goodman, 1941) 


Heat generated 
’ No. of ‘ 
Rock type areas 

eampled | “ened | 10-6 in 10-0 
cal /g. yr. cal /ec. sec. 

43 2.7 5.04.6 | 4.34.5 

Intermediate igneous ................ 7 3.0 2.04.4 | 1.94.4 

1 2.98 82.13 | .76+.12 
Columbia River basalt .............. 1 2.90 1.34.08 | 1.2+.08 


In METEORITES 
(After Paneth and Koeck, 1931, and Evans, 1938) 


Heat generated 


Type of meteorite No. o pay | 
in 10-6 in 10-8 
cal/g. yr. cal /cc. sec. 
Stony meteorites several dozen 5t 85 .52 


* Includes the Deccan trap and Columbia River basalt samples. 
+ High densities are assigned to correspond with the densities appropriate at depth in the earth, in 


accordance with Washington's Model. 


A summary of the results of Evans and Goodman (1941) in the form of 


mean values for the heat generation in the major rock types appears in 
Table 1. These averages are based chiefly upon their own work, with 
minor modifications due to their inclusion of selected determinations of 
others. Historically, the average amounts of radioactivity reported have 
tended to diminish as the subject has developed and techniques have 


Z improved. It is important to note that the present values of Evans and 
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Goodman are only about one-third those assigned to corresponding rock 
types 10 or 15 years ago. 

Radium determinations have also been made on three or four dozen 
different stony meteorites and upon a like number of iron meteorites. 
The values for the radioactive heat generation in meteorites included 
in the table are mainly dependent upon the work of Paneth and Koech 
(1931) with some weighting of earlier work as summarized by Evans 
(1938). 

Of the rock samples measured those which are probably most repre- 
sentative of the bulk of the material lying at depths of tens of kilometers 
are the Deccan traps and Columbia River basalts. Undoubtedly this basic 
basaltic material as now measured for radioactivity has lost a large part 
of its volatile constituents and probably a part of its uranium and 
thorium as it flowed into the lower-pressure and low-temperature environ- 
ment (Hurley, personal communication). Therefore, we may expect the 
mean value of about 1 & 10-1 for these basalts to be somewhat exceeded 
by the source magmas for these rocks. 

The average value of the heat generation in the random group of 43 
acid rock samples in Table 1 was not expected to be applicable to the 
granitic layer of the crust. The composition as well as the thickness 
of this layer appear to be subject to important regional variations. (This 
inference receives additional support from recent detailed seismic studies 
of the crust in New England and in Wisconsin, by means of quarry 
blasts.) It seems probable that the group of granitic samples selected 
for laboratory determinations tends to represent a more highly differen- 
tiated type of rock than that representing the bulk of the granitic layer 
of seismology. Thus, if the seven granites representing the deeply eroded 
pre-Cambrian Canadian Shield be selected from the 43 acid rock speci- 
mens in the table it is found that their mean value of heat generation 
is only 14 that of the general mean—that is to say, it is quite comparable 
with the lower values listed for basic rocks. This suggests an average 
value for the granitic layer considerably lower than the value 4.3 « 10-** 
listed. The measurements of Evans and Goodman indicate that the 
selected samples of both basic and acid rocks which are possibly most 
representative of the large bulk of crustal materials are notably lower 
in radioactivity than the listed general averages. Such preferred samples 
indicate that the heat generation in the granitic layer may be as low as 
1 or 2X 10-** cal./cc.sec., and in the basaltic layer 1X 10-78. Much 
remains to be done in establishing the average composition of the crustal 
layers and in determining the probable mean radioactivity of both the 
granitic and basaltic layers. 
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RADIOACTIVITY AT DEPTH 


It is generally held that the radioactivity must be essentially restricted 

to the upper few tens of kilometers of the crust and must diminish in 
concentration very rapidly with depth. Thus Jeffreys (1929, p. 144) 
writes, “No satisfactory alternative to this hypothesis is known, while 
several lines of evidence support it.” Also, in Holmes’ words, (1933, p. 
171), “We must conclude with Strutt (Lord Rayleigh) that the radioac- 
tive elements are practically confined to a radioactive layer which cor- 
responds with a crust a few tens of kilometers in thickness.” And (1915, 
p. 64) “. . . the view that the earth is getting hotter can not be ac- 
cepted.” Bowen (1928), however, takes the opposite point of view. “The 
plain facts are that there is nothing that can be called definite knowledge 
to preclude the possibility that the earth is heating up.” He also states 
in regard to the hypothesis of a rapid dimunition of radioactive content 
with depth, 
“This result is not acceptable. Below the granite layer there must, in any rational 
scheme of things, be intermediate and basic rocks whose radioactive content is far 
from negligible. Below these again the ultrabasic rocks must be regarded as having 
the lower but still important quantity of radioactive substances usually found in 
ultrabasic rocks.” 

There is geochemical support to the inference that radioactivity has 
been concentrated in the upper layers. However, the degree to which 
upward concentration of radioactivity has actually been accomplished 
may not well be estimated. No purely thermal criteria appear capable 
of assigning any but wide tolerances for the amount of radioactivity 
that may exist at depth. Indeed, the heat generation at depth may 
be comparable with that observed in stony or iron meteorites, or even 
with that inferred for deep crustal basalts, without leading to excessive 
heat flow at the surface. Accordingly, it is believed that no cogent 
thermal reasons exist for rejecting the hypothesis of a heating earth, or 
for concluding that an extreme degree of concentration of radium has 
occurred at the surface. Whether the tendency for heat generation to 
weaken or melt the mantle can be used to determine approximate upper 
limits for the amount of radioactivity at depth is a question still full 
of speculative possibilities. 


#R. W. Goranson (personal communication) writes, ‘‘We know that thorium and uranium will rob 
iron of its oxygen. Thus they must lie in the oxide layer. We know, furthermore, that they tend 
to concentrate in the oxygen-rich rather than oxygen-poor oxides and in fact are found most abun- 
dantly in the very latest magmatic concentrates. 

“It is to be expected then that the radioactive material of the earth will be concentrated in the 
upper layers of the earth. The evidence obtained from radioactive measurements on rocks corroborate 
such a picture. This latter evidence has little per se importance because of crustal contamination 
of the magma and thus of the rocks available to us.” 
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AGE OF THE EARTH 

Radioactive determinations of the age of old crustal rocks indicate 
that the time which has elapsed since solidification of the oldest rocks 
is about 114 to 2 billion years. With the mode of formation of the crust 
which has been described, this also represents the interval elapsed since 


TasLe 2—Thermal conductivity of some igneous rocks 
(After Birch and Clark, 1940) 


| Conductivity, in 10-% cal/em. sec. deg. C 
: No. of 
Rock type as roma | At 50°C At 200°C 
sented 

| Range Mean Range | Mean 
Granite........ 2.67 4 | 5.6-7.8 6.7 5.1-6.5 6.0 
Anorthosite.... 2.73 3 | 4.0-4.5 4.2 4.3-4.7 4.4 
Diabase....... 2.98 3 5.0-5.5 5.2 5.0-5.4 a 
Gabbro........] 2.98 3 | 47-54 | 49 [4.751 | 4.9 
Pyroxenite..... 3.23 2 | 9.2-10.0 9.6 7.8-8.7 8.3 
3.3 3 | 10.-11.4 10.5 7.5-8.8 8.1 


the mantle was essentially at its freezing point temperature throughout. 
We shall adopt, then, as the initial time for the thermal problem one 
about 114 to 2 billion years ago, when the mantle is assumed to have 
solidified. 
THERMAL CONDUCTIVITY AND DIFFUSIVITY OF ROCKS 

Birch and Clark (1940) have recently published the results of their 
excellent measurements concerning the thermal conductivity of rocks and 
minerals in the temperature range 0°- 300° C. A digest of their values 
for certain common igneous rocks is furnished in Table 2. At the higher 
part of the observed temperature range, the conductivities tend to con- 
verge, and at 200° C. all samples, igneous or sedimentary, measured by 
Birch and Clark have values between .004 and .009 cal./cm.sec.deg. 
Thus, the possible range for the thermal conductivity appears to be 
narrowly limited, if we may judge from surface results. But methods 
for extrapolation to conditions of far higher temperature and pressure 
are not yet known, and it is uncertain what the conductivity may be at 
depth. 

The thermal diffusivity is obtained by dividing the conductivity by 
the density and the specific heat. At a given temperature the specific 
heats of nearly all igneous rocks are nearly the same; for example, 
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0.20 + 0.02 at 100° C.; but the values all increase with temperature and, 
at 1000° C., lie in the range 0.3 + 0.03. Values of the density of the com- 
mon igneous rocks of Table 2 are listed in column 2 of the table. The range 
of density of these rocks is small. At depth (Bullen, 1940) the density 
of the earth is reasonably well known. It varies between 3.4 to 5.7 in 
the mantle; in the core, the mean density is between 10.6 and 10.8. The 
uncertainties in the values of the diffusivity at depth reside, therefore, 
chiefly in uncertainties regarding the conductivity and specific heat. 
Relative to the thermal conductivity, the diffusivity appears to decrease 


with depth. 
OBSERVED HEAT FLUX IN CONTINENTAL AREAS 


An evaluation of the heat flow to the surface from below requires both 
the measurement of the geothermal temperature gradient and the thermal 
conductivity of the rock. Although, as it has been noted, the conduc- 
tivities of different igneous rock types are much alike, the sedimentary 
rocks vary as to conductivity throughout a wider range. Most of the 
numerous measurements of geothermal gradients have been made in 
sedimentary rocks, under conditions where the conductivity of the rock 
was not well known. The few exceptions are the excellent recent measure- 
ments by Bullard (1939) in South Africa, Benfield (1939) in England, 
and the comprehensive studies by Hotchkiss and Ingersoll (1934) in 
the deep mine of the Calumet and Hecla Copper Company in the 
Michigan Peninsula.® 

Table 3 lists the heat flow results of Ingersoll (1934), Bullard (1939), 
and Benfield (1940). 

By “equilibrium heat flow” is meant that which would occur in a crust 
undisturbed by recent surface-temperature changes, such as those oc- 
casioned by the retreat of recent glaciation. The heat flow represents 
both the heat liberation by shallow and (to a minor extent) deep radio- 
activity and that escaping from the earth’s original store of heat. 

Much more exhaustive heat-flow studies, with careful attention to 
correlations with the other pertinent geophysical factors, are needed. 
Continental, insular, and, insofar as possible, suboceanic areas need to be 
investigated and compared.* It is commonly believed that the local 
heat flux is predominantly occasioned by relatively shallow distributions 


5 Ingersoll measured the diffusivity, rather than the conductivity, of the Calumet and Hecla 
trap rock, obtaining a value of .0075. With a density of 2.8, and a specific heat of .22, this 
corresponds to a conductivity of .0046. Birch and Clark’s measurements upon similar igneous rocks 
indicate, however, a somewhat higher value, .0052; and this value is adopted, pending results of 
measurements upon the rock itself, which it is expected will become available soon. The equilibrium 
geothermal gradient was 1.96 10 °C./em. 

* Suboceanic areas appear to differ structurally from continental areas only at shallow depths. Below 
the crust, results of seismology indicate that the earth has a high degree of symmetry. If this be 
true, it is to be expected that the radioactivity in sub ic areas will differ from that in continental 
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of radioactivity. It is important to learn what degree of correlation 
exists, in areas where the crustal structure has been well established by 
seismic methods, between the observed heat flow and the carefully esti- 
mated regional crustal heat generation. 


3.—Observed heat flow in crust 


Equilibri heat fi 

Michigan 

1.0 1.0 
South Africa 

0.97 

1.28 

Britain 

1.03 1.42+0.14 

1.20X107% 


* This location about 500 miles southwest of the others. 


The author has reviewed a considerable amount of data concerning 
observed geothermal gradients and has assigned, insofar as possible, 
rock conductivities suitable to the geological areas involved. These data 
include bore-hole gradients in the Canadian Shield (communicated by the 
courtesy of L. C. Graton), results in 6 additional North American areas, 
1 South African district, and 18 European borings. The mean heat flux 
so computed was 1.15 X 10~-°. The method is not ideal, and the result 
is of uncertain accuracy; but the sample taken was fairly large, and 
the result tends to support the view that the mean heat flow is not far 


regions only at shallow depth. Accordingly, the difference in heat flow between continental and 
Oceanic areas is probably not a fundamental factor in the cooling of the earth as a whole. Since the 
observational data for the problem are almost entirely obtained in continental areas, the thermal 
problem here being considered is really equivalent to the cooling of an earth of purely continental 
type. But an earth of oceanic type would differ only by the absence of 10 or 20 kms. of overlying 
granitic material in the outer crust. This is not a major change. 
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from 1.2 X 10~-* cal./em*sec. The subject is little explored, however, 
and regional values in the range 0.6 x 10-* to 2.5 & 10~® seem entirely 
possible. It would be most desirable to reduce this range of uncertainty 
through further field studies. 


COMPARISON OF TRANSPORT OF HEAT BY CONDUCTION AND CONVECTION 
The magnitude of the heat which could be transported by even ex- 
ceedingly slow convection currents in the mantle of the earth (if these 
exist) is great compared to that which can be transported by conduction. 
In the following, the question of the relative effectiveness of these two 
processes in the interior of the earth is examined, and the order of magni- 
tude of the convection velocities which would be thermally significant 
is established. Consider a sphere of uniform conductivity, in which con- 
vection is occurring only in meridional planes—that is, with velocity com- 
ponents % and %, the longitudinal component » being zero. Transport 
of heat under these circumstances is expressed by the equation 


aT eT , 2eT aT , aT 
Here 7 denotes the temperature; = time, p the density, and c the 


specific heat. At the center, r = 0,5 must vanish. Thus ar =? ar 
namely twice the mean value of sas on the interval (0,7). 

Now the familiar types of temperature distributions which have been 
postulated in the earth agree in letting the gradient dT/dr decrease grad- 
ually in absolute value from the surface to the center, where it vanishes— 


that is, a is uniform as to sign. Under such circumstances the average 
3 ST, has a magnitude comparable to the nee value a and 
twice this average value (which is the same as 2 dT/dr) will nowhere 
be greatly less than the local value. This inference may be easily verified 
by examination of the temperature distributions in the earth which have 
been proposed in the literature, such as those of Adams (1924), Jeffreys 
(1929), and Gutenberg (1939). In the earth the —_ temperature 
gradient a in general exceeds the transverse gradient 2 ce whereas by 
virtue of the condition of continuity of the flow the maximum values of 
the components of velocity », and v are necessarily of approximately 


the same magnitude. Thus the convective term peo, 27 represents the 


value 


greater part of the convective transport of heat, whereas the term o2 ar 
represents the order of magnitude of the conductive transport. It is clear, 


| 
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therefore, that the convective transport will predominate whenever pct, 
greatly exceeds “ ; namely when 


where k is the thermal diffusivity, o/pc. If surface measurements offer 
any guidance, the diffusivity at depth will not be large—perhaps of 
the order .01. Thus when r = 4000 kms. = 4 X 10 ems., the convective 
transport will predominate if v, greatly exceeds 5X 10-™ cm./sec. 
This velocity corresponds to a displacement of only some 16 meters per 
million years and only 32 kms. in 2 billion years—the lifetime of the earth. 

A solution by Pekeris (1935) for the thermal convection currents in an 
earth of high kinematic viscosity (3 10"'c.g.s.) yields convection cur- 
rents of velocity 3 X 10-* to 3 & 10-7 cm./sec.—namely, 60 to 6000 fold 
the limiting value above. Whether convection exists now in the mantle 
of the earth, or has ever been effective since an early stage when the crust 
solidified, are still open questions. If convection has existed, its influence 
upon heat transport very probably has far exceeded that of conduction. 
But the fact that the heat transport occasioned by vanishingly small 
convection velocities may far exceed that due to the conductivity process 
does not, of course, imply that such convective transfers are large. Both 
types may be, and in fact probably are, exceedingly small. 

The problem of a viscous earth cooling primarily by convection is 
obviously of exceedingly great complexity. The unknown mechanical 
properties of the magma, the unknown driving forces that may have 
operated, and the analytical difficulties of the problem are all obstacles 
of the first magnitude. With the remark, therefore, that convective cool- 
ing might easily predominate, we shall leave this subject and consider 
the cooling problem under the alternative hypothesis that the mantle 
of the earth has remained essentially rigid. 


PART II. THEORY OF THE COOLING OF A RADIOACTIVE EARTH 
EQUIVALENT NONRADIOACTIVE EARTH 


In the following simple procedure for solving the problem of the cool- 
ing * of a radioactive earth, the nonhomogeneous partial differential equa- 
tion is reduced to one of homogeneous type by considering the radioactive 
earth as thermally equivalent to a nonradioactive one with suitable 
altered initial temperature.* In essence, the method consists merely in 
adopting the final steady state temperature as a new zero of reference 


7™For brevity, the word ‘‘cooling,” in the present paper, will be used to denote temperature change 
of either sign—i.e., it may mean either “cooling” or “heating.” 

8 The author has recently noted that essentially this device was used by Jeffreys (1921, p. 122-149; 
1916, p. 575-591) in connection with the cooling of a flat earth at shallow depth. 
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and solving for the departure from this temperature. In the notation of 
equation (1) the heat-flow equation for a symmetrical solid earth in which 
the heat generation is H(r) cal./ce.sec. is 


aT _1a aT 

Here, H, p, c, and o may all be functions of r. The boundary conditions 
are as follows: 


(a) The surface temperature, at r= a is uniformly zero. JT =—0, 


(b) At the center of the earth, r = 0, the gradient ar vanishes. 
(c) The initial temperature distribution, T,(r), is assigned. 


When the radioactive heat term H is essentially constant in time, as 
it is in the earth, the ultimate steady-state temperature, T,, satisfies 


aT, 
o=52 (on Ze) (6) 
an equation which is readily solved for 7, by integration, namely, 
1 


Clearly, T,, must satisfy the same boundary conditions at r= 0 and 
r = aas does T(r,t). Denote by t/(r,t) the difference, T (r,t) — T,, (r), -— 


and form the difference between Equations (4) and (5). It results that 
t(r,t)satisfies the simple heat-flow equation with the radioactive term 


absent, namely, 
er 1420 ( ort Or (8) 


Par ar 
At the earth’s surface 7(r,t) also vanishes, as does = at r = 0; and 
at the initial time, t= 0, 1(0,r) = T,—T,. Thus, t(r,t) is precisely 
the temperature distribution in a nonradioactive earth, cooling to zero 
from an initial temperature which is the difference between the actual 
initial temperature and the equilibrium temperature which ultimately 
would be attained. The cooling problem for the radioactive earth is 
therefore essentially equivalent to the simple problem of the cooling of 
a nonradioactive sphere from prescribed initial temperature conditions. 
No distinction need be drawn, for the present, between that part of the 
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equivalent initial temperature which is the true initial temperature and 
that part which, in effect, expresses the presence of radioactive heat. It 
is worthy of remark that the initial equivalent temperature t(r,0) 
involves all three of the major uncertainties in the problem of the cooling 
earth—namely, this initial condition is the difference between the true 
initial temperature distribution, 7, at the time of solidification, and 
the ultimate distribution, 7, which in turn is determined only if both 
the heat generation and the conductivity be known at all depths. Fa- 
miliar results in the classical heat-flow problem of the sphere may now 
be applied in the examination of the present problem. For convenience 
the transient temperature, t(r,t), for the equivalent problem will be 
referred to as the “equivalent” temperature. 


SOLUTION FOR THE EQUIVALENT TEMPERATURE 


General statement.—Separation of the variables in Equation (8) by 
the standard substitution 


leads to the following boundary value problem for R,(r) 
(10¢) 


It is well known (see for example Lovitt, 1924, Chap. VI) that the 
characteristic numbers, 6», for this type of boundary problem are 
positive and that the fundamental functions of the boundary problem, 
Ra(Bn,r), are of oscillatory type and form a complete normalized orthog- 
onal system. 

Similarly, the functions T',,(t) satisfy 


oT. 
namely, 
Thus, the initial temperature, at t = 0 is, by Equation (9) 


where the constants A, are to be determined by the relation 


The above will suffice as a general outline of the mathematical structure 
of the cooling problem for a nonhomogeneous radioactive earth. 
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Amount of cooling, conductivity constant—The major features of the 
problem will be retained, and the essence of the matter much more simply 
arrived at, if it be assumed that the conductivity, o, and the product oc 
in Equation (10a) are constants. In view of the uncertainties about 


TasB_e 4.—Cooling after 2 x 10° years 
In an earth of thermal diffusivity .01, for various fundamental modes, n = 1 to 20. 


n 1 2 3 4 5 6 ? 
.985 .940 .872 .783 682 577 .472 
n 8 9 10 12 15 20 22 

.376 .289 .110 .031 .002 .001 


the physical properties at depth, this assumption is apparently as realistic 
as any. With o and pc constants, the problem possesses Fourier’s (1822) 
well-known solution, 


r(rt) = An exp ( - 
1 a 
where 


Equation (15) states that the equivalent initial temperature t(r,t) is 
An sin nar/a 


expressed as the sum of fundamental modes, of type , each 


of which cools (or heats, if the sign of A, be negative) to the ultimate tem- 


perature zero in accordance with its individual exponential, exp(-= it) 


Several of the fundamental modes—namely those for n = 1, 2, and 10— 


are plotted in Figure 1, each with amplitude A, = =. In the case of 


the earth, the “wave length” of the temperature perturbation correspond- 


ing to n=10 is 2x een = 1274kms. The amount of cooling after 2 


billion years in a solid earth of diffusivity 0.01 which is associated with 
each of the modes n = 1 to 10, inclusive, and n = 12, 15, and 20 
is listed in Table 4. The cooling is here expressed as the ratio of the 
present amplitude of the given mode to its original amplitude. 

It is evident that relatively little cooling has occurred throughout the 
whole lifetime of the earth in the case of the longer perturbations (n = 5 
or less). All perturbations of wave length greater than 1800 kms. (n < 7) 
retain more than half their original amplitude, and the fundamental 
mode, n = 1, has lost only 114 per cent of its original temperature. 


= 
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But when n > 10, corresponding to perturbations of wave length 1250 
kms. or less, the present temperature amplitudes in the assumed model 
are less than 20 per cent of their original values; and original tempera- 


THE FUNDAMENTAL FUNCTIONS 
An sin ane 
An=2> and n=1,2 and 


1000 2000 3000 4000 5000 6000 
RADIUS IN KM. 


Ficure 1—Fundamental functions for n=1, 2, 10 


ture perturbations of wave length 800 kms. (n = 15) or less have now 
almost entirely disappeared. The above statements, of course, depend 
upon the value, 0.01, assumed for the diffusivity. However, since the 
diffusivity enters the exponent only as a first power, whereas the factor 
n is squared, the amount of cooling is decidedly more sensitive to wave 
length than to the diffusivity. 


MODES OF MAXIMUM PRESENT RATE OF COOLING 


It has been noted that temperature perturbations of short wave length 
now retain only a small part of their original temperature and have in 
effect cooled completely, whereas those of long wave length have cooled 
scarcely at all. Evidently, there exist modes of an intermediate wave 
length for which the cooling rate is now a maximum. The modes of 


maximum present cooling rate are, indeed, simply found. Let t,(r,t) 
represent the nth fundamental mode, 


The rate of cooling, a , of this mode is 
Or sin nrr/a es 


a? 


‘ 
1.0 
CORE SURFACE 
8 
6 n=! a=6370 
4 
2 
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The rate, “ , is a maximum with respect to s when s = 1.0, and, if 


t = 6.3 X 10"* seconds (2 billion years) and k = 0.01, this implies that n is 
8 or 9—that is, the mode which would now be cooling most rapidly in 
a solid uniform earth has a half wave length of about 750 km. The 
present rate of change of temperature in this mode is 


(20) 
where t, represents the initial state of the mode, to = An 
temperature change per 100 million years is 


At = 5.84 X X 3.15 X = 1.8 X (21) 


If t) had an initial amplitude of 1000°, the present rate of temperature 
change would be about 18° C. per hundred million years, which is clearly 
very small. 

In a similar way it may be readily shown that the maximum present 
rate of change of the surface thermal gradients occurs for the modes 
n = 5 or 6. (It is here assumed, for purposes of comparing rates of 
change, that all the amplitude factors, A,, are equal; it is also assumed 
that k = 0.01 and t = 63 X 10"*, as before). The maximum rate 
turns out to be small. After 100 million years, the present gradient 
for the mode of maximum rate of change, n = 5 or 6, will have changed 
by only about 2 per cent. 


LARGE THERMAL INERTIA OF EARTH, EXAMPLES 


General considerations.—The fact that the fundamental modes for small 
values of n have suffered so small a percentage change in temperature 
during the earth’s lifetime is, of course, a direct result of the earth’s large 
size and small thermal diffusivity. Indeed, in 1822 in his original solu- 
tion of the problem of the cooling sphere, Fourier noted that, “the final 
cooling is very slow in spheres of great diameter.” However, the magni- 
tude of the heat capacity of the earth and its importance to the present 
problem are not always fully appreciated by geologists. Accordingly, it 
seems well to emphasize this feature by mention of several elementary 
examples which will convey in a simple, vivid way a concept of the 
magniture and effects of the earth’s thermal inertia. 


Heat capacity of earth—If the earth’s mean density is 5.5 and specific 
heat is about .3, then the earth’s thermal capacity is 1.8 & 10? cal./deg. 
The observed rate of heat loss is about 1.2 X 10-® cal./em.? sec. or over 
the entire surface, about 6.1 X 10” cal./sec. This loss rate, maintained 
throughout the lifetime of 2 billion years, would suffice to cool the entire 
earth by only some 200° C. 


ae 
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To make available the large supply of original heat, it is only necessary 
that adequate means, in the form of suitable initial thermal gradients and 
values of the conductivity at depth, be provided for transporting heat 
to the surface. It is certain that a nonradioactive earth model could be 
designed, if desired, which would be compatible with all the thermal 
conditions on the problem, as we now know them—namely, the condi- 
tions (a) of zero surface temperature (b) of agreement with the known 
surface heat-flux and gradients at the present elapsed age of the earth, 
and (c) the looser, more vague condition of “reasonableness” for the initial 
internal temperature distribution, and for the values of the thermal con- 
ductivity and diffusivity at depth. However, in view of the known 
existence of large amounts of radioactive heat generation, a proposed 
model for a nonradioactive earth possesses little interest. The point 
which it is desired to emphasize is merely that the heat capacity of the 
earth is large—great enough to supply the total estimated heat losses 
without aid from radioactivity and without notable loss in mean tem- 


perature. 


Slowness of thermal diffusion in earth—On a uniform solid earth 
having the thermal diffusivity observed in surface rocks, let the surface 
temperature be suddenly changed, and maintained thereafter at a new 
value. In such a case it is readily found that the temperature of the 
earth’s center will remain essentially unaffected (that is to say, to an 
extent less than 10 per cent of the imposed surface temperature change) 
for 200 billion years, or about 100-fold the earth’s lifetime. If the dif- 
fusivity at depth be 10-fold the surface value, the diffusion process is 
correspondingly speeded; an interval of 10 times the earth’s lifetime 
must elapse before the 10 per cent change is produced. The size of the 
earth is so great, and the thermal diffusivity of its materials apparently 
so small, that its interior is nearly thermally independent of its exterior. 


Cooling rates by inspection of heat-flow equation—A curve purporting 
to represent present temperature distribution in the solid earth must 
satisfy the following four obvious conditions. At the surface it has 
(1) the value zero, with (2) an initial gradient to correspond with ob- 
served values, (3) at the center, the slope is zero, and (4) the central 
temperature is unknown, but “reasonable.” By “reasonable” is meant 
a temperature measured in thousands of degrees, say between 3000° C. 
and 15,000° C. If an additional assumption (which seems universally to 
have been made) be accepted—namely, that the temperature always in- 
creases with depth, but at a diminishing rate—then the smooth curves 
which represent any temperature distribution of the character described 
must necessarily have small values of the curvature. In curves which 


> 
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have been proposed by Jeffreys, Adams, and Gutenberg, and in arbi- 


trary curves sketched by the writer (drawn for a central temperature of 


5000° C.), the values of the term or yet at depths of 500 km. or 


greater are about —1 to —4X 10-**. Thus, in accordance with the heat- 
flow equation (Equation 8), cooling by conduction occurs at the approxi- 


mate rate 2 = —k (1 to 4 X 10-**). This, with k = .01, corresponds 


to a temperature drop of between 30° C. and 120° C. per billion years. 
At greater depths the cooling so estimated is decidedly less. Such small 
rates of heat loss by conduction would be compensated by a radioactive 
heat generation of 1 to 4X 10-** cal./cc.sec—namely, about 1 per cent 
of that observed in surface basalts or granites. Thus, if such an earth is 
cooling at all at depth, the rate must be exceedingly slow. On the con- 
trary, it is clearly possible for heating to occur at a rate more than ten- 
fold that mentioned. This would mean temperature increases of 300° C. 
to 1200° C. or more per billion years. 


MANTLE PROBABLY SOLIDIFIED FROM”BOTTOM UP 


The above three elementary examples indicate that cooling rates by 
conduction, in the deep interior, must necessarily have always been very 
small. This conclusion is, of course, a matter of common knowledge 
and is in agreement with those which have already been noted in con- 
nection with the cooling of the fundamental modes. Such results show 
that the solid mantle of the earth at depth cannot have cooled significantly 
during the earth’s lifetime. Therefore, if the mantle was initially molten, 
it probably solidified progressively from the bottom up, beginning at 
the boundary of the molten core. The initial solidification of a great 
thickness of the mantle either at the surface or at some intermediate 
depth would have introduced effective thermal insulation. Subsequent 
cooling of the magma below would have been essentially halted thereby; 
and unless by chance this magma had already nearly cooled to its freez- 
ing point at the time when solidification began above, it would have 
remained fluid until today. The hypothesis that the mantle was once 
fluid, and is now almost everywhere solid, therefore implies either (1) 
that freezing began at the bottom, with subsequent rapid cooling of the 
magma above by convection, or (2) that the mantle had cooled nearly 
to its freezing point at all depths below that where solidification first 
occurred and was thereby able to congeal without much additional loss 
of heat. Clearly hypothesis (1) implies a less delicately balanced type 
of condition, and a priori possesses greater probability of truth. Fur- 
thermore the rapid cooling implied by hypothesis (1) would not be 
materially affected by the presence of radioactivity, whereas hypothesis 


te 
we 
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(2) is tenable only if the heat generation at depth be less than a few per 
cent of that in surface rocks, for otherwise no cooling could occur. 


EFFECT OF FLUIDITY OF THE CORE 


The preceding discussions neglected consideration of the probable 
fluidity of the core. Fluidity, aided possibly by superior thermal con- 
ductivity, would tend to preserve uniformity of the core’s temperature. 
It would, in effect, make more available the large thermal capacity of the 
core for maintaining constancy of temperature at the inner boundary of 
the mantle. However, the effects of fluidity in the core probably are in- 
significant because the temperatures deep in the mantle in any case change 
only slowly with time for the modes of long wave length. In the funda- 
mental mode it has been seen that the cooling for a solid earth amounts 
to only 1 or 2 per cent. 


EQUILIBRIA TEMPERATURES AND ASSOCIATED HEAT GENERATION 


General statement—The answer to the important question—Is the 
earth heating or cooling at depth?—clearly depends upon knowledge 
of two things: (1) the initial temperature distribution in the earth, and 
(2) the equilibrium temperature ultimately to be established by the radio- 
active heat. If the former temperature exceeds the latter everywhere, 
the earth is cooling everywhere; if the latter exceeds the former, it is 
heating everywhere; and if the curves representing the two distributions 
intersect, cooling occurs at some depths and heating at others. The rate 
of heating or cooling depends upon the difference between the two curves. 
If by chance the initial state nearly coincided with the final one, the 
internal temperature would, of course, remain essentially unchanging with 
time. 

According to the process for the formation of the mantle outlined in 
Part I, the initial temperatures in the solid mantle are those of the depth 
freezing-point curve of its materials in the physical and chemical environ- 
ment under which the freezing occurred. This curve is certainly not 
satisfactorily known, but the initial temperatures are generally assumed 
to increase from about 800° at the base of the crust to several thousand 
degrees at the base of the mantle. 

The final temperatures depend in a simple way (see Equation 6) 
upon the radioactivity and conductivity and may be computed for a 
range of different hypotheses concerning these quantities. Examples of 
steady-state temperatures for some dozen earth models of different types 
will be discussed. These models are of the following general classes: 

(1) Layered Earths, in which the radioactivity and the conductivity 
is uniform in shells but changes discontinuously at the shell boundaries. 
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(2) Earths of constant conductivity, in which the radioactivity de- 
creases exponentially with depth. 

(3) Earths of constant conductivity in which the radioactivity varies 
in accordance with certain other useful formulae. 


Layered earths—The surface heat flow in the layered earths chosen for 
computation has generally been set approximately equal to that observed 
today in the crust. This clearly is an arbitrary choice—the flux today 
may still be notably below its equilibrium value, and the temperatures 
computed in the examples may accordingly be too small at depth. The 
two outermost layers in the earth models are intended to represent the 
crust. The deeper layers bear no relation to actual earth structure. They 
represent merely assigned transitions in the heat generation and thermal 
conductivity for the purpose of approximating a desired type of variation 
for these quantities. The radioactivity is assumed to decrease rapidly 
in successive layers at depth, and the conductivity is chosen to increase. 
Since the temperature increments depend merely upon the ratios Hi/o, 
(see Equation 22), the type of layering chosen tends to minimize the 
steady-state interior temperatures. 

Let r; denote the outer radius of the shell throughout which the con- 
ductivity is o;, and the heat term H,;. When H;, and 9; are constants, inte- 
gration of Equation (5) yields, for the temperature 7;(r) in the ith layer, 


Cu + Cyr =r Tilr) (22) 


The constants C,,, C4 are determined by the boundary conditions at 
the interfaces of the shells, which are 


OTs 41 


The constants C,; and C.; are conveniently evaluated successively from 
the surface down, starting with temperature zero, and an assigned heat 
flux at the outer surface, r = a. The total heat generation in all layers 
is clearly fixed by the assigned surface flow at r = a. 

In Table 5 are listed the steady-state temperatures produced at the 
depths indicated, in earths of the specified layer thicknesses, heat gen- 
eration, and thermal conductivity. The influence of the thermal con- 
ductivity is separately indicated by including comparison temperatures 
for the case of constant conductivity, s = 0.01. The majority of the 
functions tabulated in Table 5 are shown plotted in Figure 2. In com- 
paring these plots, it should be noted that the temperature and heat 
generation scales in cases (1) and (2) differ from the rest; the scales 
in the cases (3) to (11), inclusive, are the same. 
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Taste 5.—Examples of equilibria temperature distributions 


For various types of radioactive earths 


Case (1) Surrace Heat Fivux, — 1.375 X 10-* sEc. 


Depth Heat Heat Thermal Depth Temp. Temp. 
km. eal. /cc. yr. eal. /ec. sec. conduct. km. °C. = 01 

0-10 25.5 8.1 -0055 0 0 0 
10-32 8.1} xX 10-6 2.6} X 10-8 -0062 10 177 97 
32-6370 0 32 279 160 
Center 279 160 


Case (2) ExponentTiaAL Decrease Rapioactivity with Depru Surrace Heat 


Friox, — 1.2 X 10-* cau./cm.? sEc. 


Convuctiviry —.006 


Depth Heat Heat Temp. 
km. cal./cc. yr. cal. /cc. sec. 
0 14.5 4.6 0 
6.37 11.3 3.6 117 
12.74 8.9 2.8 205 
63.7 13. px 434 
127.4 -116 .037 528 
254 .000 0 532 
Center 0 0 532 
Cases (3), (4) anp (9) Surrace Heat Fivux as In (1) 
Depth Heat Hea Conductivity Temperature Temp. | Depth 
km. cal./ec. yr. cal./cc. sec. 3) (4) (9) (3) (4) (9) |co=0 km, 
0-15 12.1 3.8 -0055 .0055 .0055 0 0 0 0 0 
15-40 4.0 | -0124 .0062 .0062 297 297 297 163 15 
40-200 -178} X 10-*] .057}X 10-3) .024 024 .008 428 559 559 326 40 
200-570 -080 .025 030 03 .01 734 865 1477 1060 200 
570-6370 .061 .020 -15 15 -05 1231 1362 69 2552 570 
1954 2085 5137 | 13,392 | Center 
Case (5) Oceanic Crust. Granttic Laver ABSENT. OcEAN—4.3 KM, DEEP 
(Compare Case (6)). Surrace Heat Frvx, —1.1 X cau./cm.? sEc, 
Depth Heat Heat Thermal Temp. Temp. Depth 
km, cal. /cc. yr. cal./cc. sec. conduct. = 01 km, 
42 “008” 475 204 34:3 
194-564 | ¢ | ‘o50¢ X 108 ‘010 1556 1158 194° 
564 -6366 -006 .002 .05 2057 1659 564 
2271 2729 Center 
Cases (6) anp (7) Surrace Heat Fiux, —1.375 X 10-* cau./cm.* sec. as In (1) 
Depth Heat Heat Conductivity Temperature Temp. | Depth 
km. cal./ce. yr. eal. /cc. sec. (6) (6) (7) o = 01 km, 
0-10 8.64 2.7 -0055 .0055 0 0 0 0 
10-40 2.7 .86 -0062 .0062 225 225 124 10 
40-200 1.31 10* 42 X10-8 -08 .08 700 700 418 40 
200-570 .16 -05 01 .0048 1781 1781 1282 200 
570-6370 ‘ .05 .05 282 2782 1783 570 
2496 2996 Center 


— 


L. B. SLICHTER—COOLING OF THE EARTH 


TasLe 5.—Ezamples of equilibria temperature distributions—Continued 


Case (8) Ocgeanic Crust. Granitic Layer ABSENT. OCEAN 4.3 KM. DEEP. 


(Compare Case (9)). 


Surrace Heat —.805 X 10-* cax./cm.? 


Depth Heat Heat Thermal Temp. Temp. Depth 
km. eal./ce. yr. cal. /cc. sec. conduct. 

29 057 “008. 282 162 29'3 
189.3-559.3 ‘og ¢ % 10° ¢ % 10" “010 1180 896 | 189.3 
559 .3—Center -061 .020. -050 2672 2388 559.3 
4840 13 ,228 Center 


Case (10) Surrace Heat Fivux, —2.75 X 10-* cau./cM.? 


(Compare Case (1)). 


Depth Heat Heat Thermal Temp. Temp. Depth 
km. cal. /cc. yr cal./cc. sec. conduct. =. = 01 km. 
10-40 ye 36 “0062 427 235 10 
40-200 2.4¢ % 10° 10" "0080 1187 706 40 
200-—Center 0. 2229 1540 120 
1666 200 
2387 1666 Center 
Case (11) Unirorm Raproactiviry, —.179 X 10-* cax./cc. yr. = .057 X 107" caL./cc. BEC. 
Unirorm Conpuctivity, —.006 
Surrace Heat Fiux, — 1.2 X (as (2)). 
2 
TEMPERATURE = 63,700 € = =) 
Teme. Temp r/a 
0 0 12,100 9 
15.9 318 
31.8 637 47 ,800 5 
47.8 955 61,200 
64.3 1274 63,700 .0 
128.5 2550 
194.0 
327.0 6370 


The cases (1) to (11) are arranged in order of increasing radioactivity 
at depth, and so involve progressively higher final internal temperatures. 
Except as will be noted, the surface heat flux is always either 1.2 x 10-° 
or 1.375 & 10-® cal./em*sec. Cases (5) and (8) are the same as cases 
(6) and (9), respectively, except that the “granitic” layer has been 
omitted, to obtain an earth of “oceanic” type. The surface heat flow 
is correspondingly reduced from 1.375 X 10~-* to 1.1 10~* cal./em.’sec. 
in case (5) and to 0.80 * 10~* in case (8). In case (1), no heat gen- 
eration is assumed below a depth of 32 kms. Case (10) is like case (1) 
except for the addition of heat generation at depth, in particular in 
the zone 40 km. to 200 km. where it is assumed that the heat generation 
is 0.76 X 10-'* cal./em.’see., or midway between that of the Deccan- 
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Columbia River basalt average, and that for stony meteorites. The 
equilibrium surface heat flux is thereby raised to 2.75 X 10-° cal./em2 
sec., and the central temperatures from 280° C. to 2390° C. A further 
illustration of the high ultimate temperatures which may be produced 
by a minor content of radioactive heat is seen in case (11). A uniform 
radioactivity of about 6 per cent that observed in the surface rocks 
lowest in radioactivity leads to ultimate central temperatures of over 
60,000° C. 

The exponential distribution of radioactivity, of which case (2) is 
an example, is considered in greater detail in the following paragraph. 


Exponential distributions —The exponentially decreasing type of heat 
distribution is often assumed (Ingersoll and Zobel, 1913; Lowan, 1933; 
Van Ostrand, 1940). If the heat generated is 


(23) 


then, by Equation (6), the steady-state temperature in an earth of con- 
ductivity o is: 


2 2 
r-% 1-2 -(1 - exp Ba (1 wt + exp (— pa) ]....(24) 
where y = r/a. 
At the center, 
T.= ) exp (— 8a) | (25) 
The equilibrium heat flux at the surface is 
oT dH, 2 


Usually the arbitrary constant, 8, is chosen large enough to render 
negligible the exponential term exp (— a). 
In case (2) of Figure 2a the plotted temperature corresponds to the 


values 


H, = 4.6 X 10-® cal./cc. sec. 


Ba= 242, 6? = 1.45 X 10°" 


The corresponding surface heat flow is 1.2 X 10~-® cal./em.’sec.—i. e., 
about that observed at present, but the radioactivity in the example 
decreases exceedingly rapidly from a surface value of 4.6 & 10~**, cor- 
responding approximately to the Evans-Goodman value in granite, to 
less than 10 per cent of this value at a depth of only 125 km. The ulti- 
mate temperature is accordingly nearly constant at depths below 100 
km., at the low value 530° C. Evidently such an earth would be cooling 


fi 
we 
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at depth. With other, and possibly more tenable, values for the constants, 
a high ultimate internal temperature is readily obtained. For example, 
let Ho. = 10-**, corresponding to the heat generation in Deccan traps 
and Columbia River basalts; let o— 0.01. Ba=19. In this case the 
ultimate central temperature of 10,000° is 99 per cent attained at a depth 
of 1400 km. The equilibrium heat flux is 3 & 10-°, and the radioac- 
tivity drops to 10 per cent of its surface value in a depth of 700 km. 
Because of the high ultimate internal temperature, this example prob- 
ably represents a heating earth. 

The transient temperatures for the exponential type of distribution 
have been expressed and computed by Lowan (1933) and Van Ostrand 
(1940). Lowan’s expression for the temperature rise due to the genera- 


tion of radioactive heat is 
n? 


sin [1 — exp 


2 H.a*® 
n (Pa + 


Lowan computed the transient temperatures from Equation (27) for 
the values H, = 10-'* ecal./ce.sec., k = 0.005 and 0.007, o = 0.006, 
Ba = 448, B = 6.66 10-7, and t = 1.6 X 10° years. These values 
correspond to an ultimate surface heat flow of 1.5 « 10~°, and a steady 
central temperature of 375°. The surface gradient in this example has 
essentially achieved its ultimate value, but the temperature at a depth 
of 200 km. is still only about 50 per cent of its final value, whereas at 
30 km. depth it is about 114 times the steady value. 

The sum in formula (27) converges very slowly, especially for small 
values of the depth, a - r, and for large values of Ba; in one of Van Ost- 
rand’s examples (Ba = 6,370) the series was summed to 6740 terms, for the 
case a - r= 1 km., to 2540 terms for the depth 50 km., and to 160 terms 
for depth 200 km. The difficulty arises partly ® from the fact that the 
initial coefficients decrease, when fa is large, only as 1/n until an ap- 
proaches Ba. The assumption of an exponential distribution is not an 
especially convenient one, and other, more tractable choices are recom- 


mended. 


Certain convenient types of heat distributions —The following four dis- 
tributions of radioactivity lead to simply computable expressions both 
for the steady-state temperatures and for the transient temperatures. 
Moreover, when taken in combination, with amplitude factors of appro- 


* For an exponential distribution in which the heat generation is concentrated near the surface, it is 
well known that the solution at shallow depths is conveniently computed by adopting the approxi- 
mation of a flat earth. 
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priate size and sign, they may be used to approximate a wide variety of 
suitable heat-generation functions. These functions are the following: 
(1) The fundamental functions, 


sin nar/a 
(28) 
to which correspond the ultimate temperatures 
T= r, sin nar/a , @ Ho (29) 


nar/a 
Here H, denotes the rate of heat generation at the center of the earth, 


and 7’, its ultimate central temperature. 
(2) The distribution 


to which corresponds 


net 7, He 


= T. cos 


The mode for n= 1 is especially useful, and in subsequent reference to 
eases (1) and (2) the value n=1 will be presumed, unless otherwise 
indicated. 

(3) The constant distribution, 


to which corresponds 
T =T, (1 —r/a), T. = (33) 


(4) The heat distribution which is constant throughout a central core 
of radius r, and zero elsewhere. Thus, 


(34) 
To this corresponds the steady-state temperature 
Hy 2 
Tro = 42(1-2)-(£) cr en (35a) 
3 


(4a) From cases (3) and (4) is obtained by subtraction the steady 
temperature distribution due to heat generation uniform throughout a 
superficial shell of thickness h = a - r,. This is a special case of the steady 
temperatures in layered radioactive earths which have been listed in 
Table 5. 

In Figure 8, in solid line, are plotted the first three temperature func- 
tions above, each with the same central temperature, and the correspond- 


H nar Qa 
3 2a nr (30) 
2 
an 
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ing heat generation functions in dotted line, for comparison. The three 
temperature curves are nearly alike, and they obviously are not sensitive 
to the differences in their respective generating functions. 

The steady-state temperatures illustrated in Figures 2 and 3 are in 
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Ficure 3.—Certain steady-state temperatures and associated heat-generation functions 


T, = temperature at center; H, = heat generation at center 


many instances very large at depth. If the assigned amounts of radio- 
activity (which are very small, judged from observational experience) 
be accepted, then the probability is strong that the earth is heating at 
depth. It is also apparent from inspection of these temperatures that 
a curve representing the initial temperature in the mantle might intersect 
at depth, in appropriate cases, the steady-state temperature curves of 
Figures 2 and 3, thus leading to cooling at certain depths and heating at 
others. Whether distributions of radioactivity of the small concentra- 
tions suggested in these examples are compatible with the observable 
thermal data at the surface will be considered hereinafter. 


TRANSIENT TEMPERATURES 


General remarks.—Formulae for the transient temperatures T'(r,t) 
existing in earths being heated from zero by the above four types of 


= 
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radioactive distributions—Equations (28), (80), (82), and (34) re- 


spectively—are readily caleulated.*° (See Equations (15) and (16).) 
They are as follows: 


(1a) T¢,0 = [1 - exp (- (36) 


narja 
| (87) 


_ Ho a’ (—) ""12n sin nxr/a 
(2a) T(r, t) = 


Ho 2 a 


(4a) T (r,t) = 


sin nxh/a+nr(1 


sin nar/a [1 —exp ( )] (39) 


In the formulae (37), (38), and (39), the convergence of the series 
is satisfactorily rapid, by virtue of the factor 1/n? in the typical term. 
Clearly the transient temperatures for earths of constant conductivity 
in which the radioactivity is distributed uniformly in shells are readily 
found by forming the appropriated sums and differences of formulae of 
the type (38) and (39). Similarly, obvious combinations of formulae 
(36), (37), (38), and (39) may be formed to obtain the transient tem- 
peratures for a wide variety of radioactive heat distributions. 


Special case. Simple relation between transient temperatures and heat 
term.—lIf the heat-generation function vanishes at the surface, then its 
corresponding transient temperature function may be written down at 
once, in terms of the coefficients expressing the heat-generation function. 
To see this, introduce the auxiliary heat function, h(r), defined by 


Then both h(r) and T,,, (r) satisfy conditions similar to those (see Equa- 


tions 10a, 10b, 10c) determining the fundamental functions, R,, of the 
cooling problem; namely, h(r) and T, (r) satisfy 


Te atr=0. (41b) 
or or 


10 For computing temperatures at shallow depths, z, it is desirable to use the equality, 


sin ut = (—)®*1 sin a, in formulae 37, 38 and 39. 


a 

I 
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Clearly h(r) (if continuous) and T ,.(r) may be expanded in terms of 
the fundamental functions, R,. Thus, 


From inspection of Equations (10c) and (41) it is evident that 


Thus, when the radioactivity vanishes at the surface, the simple relation 
(44) holds between the amplitudes A, of the elementary modes of the 
steady-state temperature distribution and those, C,, expressing the 
distribution of radioactive heat. More specifically, if the initial tempera- 
ture is uniformly zero, then the transient temperature is 


T(r,t) = > Re — exp(—B*, ») (45) 


In case both o and pe are constants, formula (45) reduces to one similar 
to Equation (36), namely, 


where the heat function has the expansion 
nr 


n=1 
SURFACE HEAT FLOW 


Equations (36), (37), (38), and (39) express transient temperatures 
in earths being heated from zero by the respective heat-generation func- 
tions expressed in Equations (28), (30), (82), and (34). To take into 
account the effect of the earth’s initial temperature, T,(7,0), the addition 
of a transient term, 7',(r,t) is needed. This term represents the tempera- 
ture in a nonradioactive earth cooling to zero from the initial state, 
T,(r,0). The surface heat flows, F;, associated with each of the solutions 
(36), (37), (38), and (39), amended by inclusion of the term T,(r,t), are 
respectively 


F, = [1 exp (-" kt +om.. 


1 
h= 
1 
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64 
F.= 3 


[1 exp (- = it) | +oam..... (37a) 


2a 1 
22S) - (- a 
1 


F,= + nx (1 — h/a) cos 
1 


[1 — exp at) | + om (39a) 


Here the term om, which represents the present surface flux due to the 

original heat is small compared to the term arising from the radioactivity. 

It is estimated that om, is between .2 and .4 X 10~-® cal./cm.*sec., de- 
re pending upon the assumptions as to the conductivity and the initial tem- 
e peratures. A value of .3 will be adopted as a satisfactory approximation 
3 for present purposes. The first quantities in the sums (7.e., those which do 
i not contain the exponential factor) represent the steady-state flux. This is 
best calculated directly from formulae (29), (31), (33) and (35b), re- 
spectively, for steady-state temperatures. The sums containing the 
exponential factors are rapidly convergent and were separately computed 
to 20 or more terms, using a value .01n? for the exponent, which corresponds 
to k = .0065, t = 2 billion years; likewise to k = .0086, t = 11% billion 
years. The results of these computations for distributions H(r) of type 
1’, 2’, and 3’ in Figure 3 are, respectively, 


F, — om, = 9.97 X 10° H, = 2.35 X10" (49) 
: : Here H, represents the value of the heat generation at the surface, H, 


that at the center of the earth. Formulae (49) and (50) are nearly alike. 
For values of the surface heat, H,, of order 10-**, they yield somewhat 
higher values than are now observed for the surface heat flow and indicate 
that the radioactivity must diminish, rather than remain constant or 
increase, with depth. However, the uncertainties in present knowledge 
of the surface heat flux are such that large amounts of radio-activ- 
ity may exist at great depth without noticeable effect upon the sur- 
face heat flow. For example, if in Equation (48) H, is 10—* cal./ce.sec., 
as in Deccan or Columbia River basalts, the present flux due to this dis- 
tribution is only 9.7 X 10-*—namely only about 8 per cent of that 
being observed. An additional flux of this magnitude would, of course, 
escape notice, since it is within the present observational error in such 
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measurements. If this heat were present, the rise in central temperature 
during the earth’s lifetime would be (in accordance with Equation 51), 
with o = .01, 
(61) 
o x 


The hypothesis of so great a temperature rise does not seem acceptable. 


TaBLe 6.—Surface heat flow 


Produced by a Uniform Radioactive Layer of Thickness h km. 
Heat Generation = H cal./cce. sec. 
Diffusivity = .0065 
Elapsed Time = 2 X 10° years 
Initial Temperature, zero. 


Surface heat flow 
in H X 10’ cal. /cm.? sec. 
a km. 
Present Steady-State 
0 0 0 0 
.02 40.5 .376 .403 
.05 101.4 .862 .998 
.10 202.8 1.59 1.96 
.20 405.5 2.03 3.80 
.50 1014 2.14 8.61 
T 6370 2.22 21.2 


PRESENT SURFACE HEAT FLOW DUE TO LAYERED DISTRIBUTIONS 


The present surface heat flow due to superficial layers of thickness h 
in which the heat generation rate is H may be obtained by differencing 
Equations (50) and (39a). In Table 6 are listed the present flows, and 
also, for comparison, the ultimate flows arising from layers of a number 
of different assigned thicknesses h. The results are also shown graphi- 
cally in Figure 4. The significant conclusion from these calculations is 
that the present flux diverges rapidly from the steady-state value when 
the layer thickness exceeds 200 km. When the thickness is 400 km., the 
present flux is about 14 the ultimate value, and further increase in the 
layer thickness produces almost no increment in the surface flux. By the 
aid of Table 6 or the curve in Figure 4, the present surface heat flow in 
earth models having heat generation distributed in uniform layers may 
be quickly computed. For example, consider the following model, in 
which the successive layers correspond respectively as to radioactive heat 


596 L. B. SLICHTER—COOLING OF THE EARTH 


with (1) granite, (2) Deccan and Columbia River basalt (3) stony 
meteorites, and (4) iron meteorites. 


Layer Heat generation 
(km.) in 10-8 cal/ec.sec. 
0-10 2 
10-40 1 
40-2900 5 
2900-6370 25 


SURFACE HEAT FLOW / 
Hx 10" Cal em?~ sec. 
‘Steady state Flow 
/ 


/ 
20 
/ Present Flow 
/ 
/ 
/ 
/ 
/ 

/ 


100 200 300 400 500 600 700 800 900 1000 
SHELL THICKNESS, h in KM. 


Ficure 4—Surface heat flow produced by a uniform radioactive layer of thickness 
hin km. 
Heat generation =H cal/cc. sec.; diffusivity = .0065; elapsed time 10° years; initial 
temperature, zero. 


The model closely resembles, in respect to heat generation, the earth pro- 
posed by Washington (1925; 1929). The present surface heat flow pro- 
duced by this model is 1.38 & 10~° cal./em.’see. If .3X 10-* be added to 
take account of the flux due to the earth’s original heat, the total present 
flux is 1.68 X 10-*. To see the effect of the depth to which the radioactiv- 
ity extends, let it be successively assumed that the heat generation in the 
above model is zero at all depths below 400, 200, and 100 km. respectively. 
Then the total present surface flows (including in each case an amount 
.3 X 10-* due to original heat) for the several cases are as shown in 
Table 7. 


| 

is 
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Both Tables 6 and 7, and especially the curve in Figure 4, are of use 
in furnishing quantitative illustration of several features: 

(1) An impracticably high degree of accuracy is needed to distinguish 
between the effects of different distributions of heat sources at depths 
below 200 km. An appreciation of the difficulty of estimating the amount 
of heat generation at depth in the presence of shallow heat sources is 


TaBLE 7—Effect of downward extent of heat generation upon present surface flux 


Example No. 1 2 3 4 
Downward extent of heat generation To Center 400km. 200km. 100 km. 
Present surface flux, in 10° cal./cm.’sec. 1.68 1.58 1.37 1.00 


furnished by the following example. If the heat generation in the model 
be assumed uniform throughout, then the zone within 400 km. of the 
surface will account for 90 per cent of the present heat flow of radio- 
active origin. The remainder of the earth, comprising 82 per cent of 
the volume, will contribute only 10 per cent of the flux. In fact, it is 
clear that the present heat flows produced by deep and shallow heat 
sources uniformly distributed throughout thin spherical shells are in the 
ratio of the slopes of the curve in Figure 4 at depth and at the surface. 
At depths greater than 200 km. the slope diminishes to a very small value. 
Indeed, at 400 km. the slope is only 4 or 5 per cent of its initial value, 
which means that the present surface heat flow due to radioactivity at 
depth 400 km. is only 4 or 5 per cent of that occasioned by an equal, 
shallow concentration. 

(2) At lesser depths than 100 km. or 200 km. the slope is nearly uni- 
form, being 80 per cent of its surface value at depth 100 km. and 60 per 
cent at depth 200 km. Throughout this zone, therefore, the heat genera- 
tion is in fair equilibrium with the heat loss at the surface. It follows that 
the manner of distribution in depth of the radioactivity above the 100 or 
200 km. level is without significant effect upon surface heat flows, and 
only the total amount present above this level counts. The observed 
heat flows fix definite upper limits for the total radioactivity which may” 
exist within 100 or 200 km. of the surface. For example, an observed 
flux of 1.2 < 10-® cal./em.*sec. corresponds to an average heat genera- 
tion in the first 100 km. of depth of 1.4 10-"* cal./cc. sec. or less; or to 
an average generation of .75<10-—** or less above the 200 km. level. 

Careful studies of the three pertinent quantities: (1) the surface heat 
flux, (2) the local thickness and character of the layers composing the 
crust, and (3) the mean radioactive heat generation associated with the 
crust (all in the same regions) are necessary in order to ascertain the 
degree to which known shallow distributions of radioactivity are cor- 
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related with observed heat losses. After estimation of the shallow genera- 
tion of heat and allowance for the original heat, the residual flux may be 
associated in a general nonunique manner with heat generation at greater 
depths. In no area are satisfactory observations upon all three of the 
above factors now available. Any of the values for the heat flux in the 
range 1.00 < 10-* to 1.68 X 10—° listed in Table 7 might be in approxi- 
mate agreement with such partial observations as now exist, and any of 
the four hypotheses concerning the distribution of radioactivity is there- 
fore tenable. 

The rise in internal temperatures accompanying deep-seated distribu- 
tions of radioactivity is, however, an important feature and appears to set 
certain vaguely understood limits upon the rate of heat generation that 
may exist at depth. In the model, in which the radioactivity assigned to 
the core is that in iron meteorites, the temperature rise in the core during 2 
billion years is (if pc — 4.0) about 400°C. In the mantle, the rise 
would be greater, possibly about 1500° C. If the suggested mode of origin 
of the mantle be adopted, then the rise permissible without producing 
remelting is probably measured in hundreds and not thousands of degrees. 
It then follows that the maximum permissible amounts of heat generation 
in the mantle at depth are of the order of 10—** cal./ce.sec., or about 10 
per cent of that observed in surface plateau basalts. This degree of con- 
centration of the radioactivity toward the surface is by no means implaus- 
ible and is in accord with geochemical considerations requiring the con- 
centration of the radioactive elements near the surface. 


SUMMARY AND CONCLUSIONS 


The problem of the cooling of a radioactive earth has been reviewed. 
It is probable that the mantle solidified from the bottom up, beginning 
at the boundary of the core. More exhaustive studies of the surface heat 
flow, of the structure and composition of the crust, and of the mean 
heat generation in crustal layers, all in the same regions, are desirable. 
The surface heat flows fix definite upper limits for the amount of radio- 
activity which may exist within 100 or 200 kilometers of the surface. But 
thermal observations and theory seem incapable of furnishing information 
about the thermal state of the earth at depths greater than 200 or 300 
kilometers. Little correlation exists between the radioactivity at depth and 
observed surface heat flux. Heat generation of the order 10-"* cal./ce.sec., 
or 10 per cent of that observed in plateau basalts, may exist at all depths, 
to the center, without producing effects contrary to the facts of observa- 
tion. It is unknown whether the earth is heating or cooling at depth. 
No need is found for a drastic rate of decrease of radioactive substance 


‘ 
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with depth. The oft expressed idea that the indicated amounts of radio- 
activity in the earth are embarrassingly large and greatly exceed the 
amounts to be expected from observed surface heat flows is fallacious. 
Hoskins (1839) summarized the thermal problem in these words, 


“investigation of the earth’s refrigeration . . . still leaves us in a state of perfect 
uncertainty as to the actual condition of its central parts ... from want of the 
experimental determination of values which it must ever be found extremely 
difficult, if not impossible, to attain with accuracy.” 
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